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Abs  tract 

Space  facilities  of  the  next  decade  will  require 
applications  of  highly  advanced  materials  which  have 
properties  exhibiting  excellent  fracture  toughness  and  high 
strength-to-weight  ratios.  The  recently  introduced 
thermoplastic  composite  material,  graphite  polyetherether- 
ketone  (Gr/PEEK)  APC-2,  promises  lower  costs,  lower  part 
weight,  and  higher  operating  temperatures.  This  new  class 
of  organic  material  has  fracture  toughness  properties 
superior  to  those  of  graphite  epoxy.  This  thesis  examines 
the  failure  characteristics  of  Gr/PESK  through  an 
experimental  investigation  and  through  the  application  of  a 
fully  nonlinear  p.ly-by  ply  finite  element  technique. 

Laminates  investigated  were  [0  3  and  [90  3  unidirectional 

i  «s  1  <s 

lay-ups  and  [±453  angle-ply  lay-ups.  An  experimental 
investigation  was  also  done  on  [ 0 / 45 / 90 / - 4 5  3  and 
'0/±45/90  3^^  quas 1  -  isotropic  lay-ups  to  study  failure 

The  experimental  investigation  of  Gr/PEEK  APC-2 
involved  the  testing  of  34  tension  and  compression  coupons 
to  derive  basic  material  properties  for  use  with  the  finite 
element  program.  To  investigate  failure  characteristics 
experimentally,  78  tensile  specimens  were  manufactured  with 
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0.4-inch  diameter  holes.  A  portion  of  these  unidirectional, 
shear,  and  quas i - i so trop i c  specimens  where  tested  at  room 
temperature  to  ultimate  strength.  The  remaining  specimens 
were  tested  to  percentages  of  average  failure  stress  to 
i.nvestigate  the  initiation  and  progression  of  laminate 
failure.  To  supplement  the  test  results,  two  post-failure 
analyses  were  conducted.  Videotapes  were  made  of  each  test 
to  ultimate  strength.  Also  stereo  x-rays  were  taken  of  each 
test  specimen  subjected  to  a  percentage  of  failure  load. 

For  comparison  to  experimental  results,  a  fully 
nonlinear  progr ess i ve-p ly- f ai lure  finite  element  program 
was  employed.  The  models  of  stress-strain  responses  and 
growth  of  failure  were  found  to  closely  approximate  the 
results  of  the  experimentation. 

This  investigation  provided  futher  data  on  the 
application  of  tensile  loads  to  Gr/PEEK  containing  circular 
discontinuities.  This  study  also  proved  that  a  nonlinear 
finite  element  program  can  closely  approximate  progressive 
ply  failure  in  a  Gr/PEEX  laminate.  The  excessive 
nonlinearity  of  this  material  proved  the  need  for  using 
nonlinear  techniques  when  analyzing  laminates  of  Gr/PEEX. 
Furthermore,  this  research  reinforced  the  proposition  that 
the  thermoplastic  matrix  does  produce  a  reliable  composite 
that  should  be  considered  for  use  on  aircraft,  spacecraft, 
and  space  facilities. 
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A  STUDY  OF  FAILURE  CHARACTERISTICS  IN 
THERMOPLASTIC  COMPOSITE  MATERIAL 

I .  Introduction 

Man  has  studied  materials  and  their  failure  mechanisms 
since  he  began  building  facilities,  seeking  transportation, 
and  manufacturing  goods.  Knowing  how  a  material  fails  under 
given  loading  conditions  leads  to  improvements  in  material 
capabilities.  The  study  over  time  of  earth,  wood,  metals, 
and  plastics  has  led  to  more  advanced  materials;  and  often 
materials  have  been  combined  to  produce  optimum  effects. 

Over  the  past  few  decades,  this  idea  of  combination  has  led 
researchers  to  develop  advanced  composite  materials. 

In  the  next  decade,  space  facilities  and  spacecraft 
will  require  application  of  highly  advanced  materials  that 
have  properties  which  exhibit  high  fracture  toughness,  high 
s treng th - to- we i gh t  ratio,  and  low  thermal  expansion.  The 
recently  introduced  polyether-based  thermoplastic  composite 
material,  graphite  polyetherether-ketone  (Gr/PEEK)  promises 
lower  costs,  lower  part  weight,  and  higher  operating 
temperatures.  This  new  class  of  organic  material  has  a 
fracture  toughness  much  higher  than  that  of  graphite  epoxy 
(Gr/Ep)  [1].  In  general,  thermoplastics  promise  lower 
production  costs  and  improved  maintainability  over  similar 
thermosets  such  as  Gr/Ep  [21. 


1 


By  studying  the  stress-strain  response  to  ultimate 
strength  and  the  initiation  and  progression  of  failure  in 
this  new  composite  material,  this  thesis  contributes  to  the 
vigorous  research  effort  being  conducted  by  many  groups  to 
evaluate  Gr/PEEK  thermoplastic  composite  material  for 
practical  applications. 

A .  Purpose 

The  purpose  of  this  thesis  is  to  investigate  the 
initiation  and  progression  of  failure  in  Gr/PEEK  containing 
a  circular  discontinuity  and  subjected  to  tensile  loading. 
With  this  purpose  in  mind,  the  objectives  of  this  study 
wer e  : 

(1)  To  determine  basic  material  properties  for  the 
Gr/PEEK  at  room  temperature  (RT) , 

(2)  To  determine  the  experimental  s tress - s trai n 

response  of  [0  ]  and  [90  ]  unidirectional  laminates, 

1  <S  1  <3 

[±45]  angle-ply  laminates,  and  [ 0/ + 4 5/ 90 / - 45  ]  and 
[0/±45/90]^^  quas 1 - isotropic  laminates  of  Gr/PEEK  specimens 
containing  a  0.4-inch  diameter  circular  discontinuity  and 
loaded  to  ultimate  strength  at  RT , 

(3)  To  experimentally  investigate  the  initiation  and 
progression  of  failure  in  similar  specimens  containing  a 
hole  and  loaded  to  percentages  of  average  failure  stress, 

(4)  To  analytically  model  the  failure  process  of  the 
unidirectional  and  angle-ply  laminates  using  a  fully 
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nonlinear,  progress i ve -p 1 y- f ai 1 ure  finite  element  program, 
and 

(5)  To  compare  experimental  results  of  the 
quas 1 - isotropic  specimens  and  to  investigate  their  potential 
for  delaraination  by  comparison  to  Qr/Ep. 

These  objectives  were  accomplished;  and  the  report  of 
all  theories,  procedures,  results,  and  conclusions  is 
contained  in  this  thesis. 

B .  Background  and  Overv i ew 

A  composite  material,  as  defined  in  this  thesis,  is  a 
material  consisting  of  reinforcements,  such  as  fibers,  in  a 
continuous  matrix,  such  as  a  polymer.  These  two 
constituents  are  combined  on  a  macroscopic  scale  to  form  a 
useful  material. 

Composites  have  a  long  history  of  use.  Prehistoric 
civilizations  used  straw  in  mud  bricks;  plywood  existed  in 
early  Egyptian  construction;  and  medieval  knights  used 
combinations  of  various  metals  in  their  swords  [31.  In  the 
last  century,  reinforced  concrete  began  to  give  engineers 
greater  applications  in  construction. 

Conventional  composite  materials,  or  'reinforced- 
piastic'  (RP) ,  was  first  commercialized  as  fiberglass  in 
1939.  Fiberglass  still  accounts  for  90%  of  the  RP  market 
[4].  But  the  emphasis  on  research  in  the  past  few  decades 
has  been  on  advanced  composite  materials. 
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The  adjective  'advanced'  is  used  to  distinguish  between 
composites  of  new  ultrahigh  strength  and  stiffness  fibers 
and  composites  of  familiar  fibers  such  as  glass.  The  matrix 
of  these  new  composites  can  be  either  a  plastic,  such  as 
epoxy,  or  a  metal,  such  as  aluminium.  These  advanced 
composites  have  two  major  advantages  among  many  others: 

(1)  improved  strength-to-weight  and  s 1 1 f f ness  -  to- we i gh t 
ratios  and  (2)  the  ability  to  tailor  a  composite  to  meet 
design  requirements. 

Composite  components  can  be  made  that  have  the  same 
strength  and  stiffness  as  high  strength  steel  components  but 
are  70%  lighter.  Some  advanced  composites  are  up  to  three 
times  stronger  than  aluminum  but  are  60%  lighter.  These 
advanced  composite  materials  are  ideal  for  we i gh t -sens i t i ve 
structural  applications  such  as  space  structures  and  space 
vehicles.  Some  examples  of  current  use  include  the  motor 
casings  on  the  space  shuttle  solid-fuel  rockets  [4]  and  the 
wings  of  the  X-29  f orward -swept -wi ng  developmental  aircraft. 
The  F-15  and  F-i6  have  proven  the  value  of  composites  in 
their  use  as  horizontal  and  vertical  stabilizers. 

Tailoring  a  composite  to  meet  strength  and  stiffness 
requirements  in  particular  directions  leads  to  economy  of 
manufacturing.  Because  fibrous  composites  are  anisotropic, 
i.e. .  material  properties  are  different  in  all  directions, 
components  can  be  designed  for  strength  and  stiffness  in  a 
particular  direction.  There  is  little  waste.  Part 
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consolidation  and  reduction  of  pos  t- f  abr  i  cat  i  on  machinin;^ 
provide  up  to  50%  cost  savings  over  analogous  machined-metal 
assemblies  [3].  Conventional  isotropic  materials  must  be 
designed  for  the  worst  case,  resulting  in  excess  strength, 
stiffness,  and  weight  in  all  other  directions  [3].  These 
advantages  of  composite  materials  promise  to  lead  to  new 
aircraft  and  spacecraft  designs  that  significantly  depart 
from  past  efforts  based  on  conventional  isotropic  materials. 

In  the  composite  industry,  the  dominant  reinforcement 
of  composite  materials  is  continuous  fibers  (as  opposed  to 
whiskers,  filaments,  or  spheres  of  reinforcing  material)  . 
Throughout  this  thesis,  reference  to  composites  will  imply 
cont 1 nuous - f iber-re inf orced  composite  materials.  Also 
throughout  this  thesis,  reference  to  composites  will  imply 
the  use  of  a  polymer-based  matrices,  such  as  epoxy  or  PEEK. 

In  applications  of  composite  materials,  the 
conventional  problems  of  design  and  analysis  arise.  As  in 
any  structure,  composite  components,  or  laminates,  have  to 
have  holes  to  serve  various  purposes.  In  practice, 
composite  components  are  connected  through  adhesively  bonded 
joints,  through  bolted  joints,  or  through  a  combination  of 
each.  Whether  caused  by  bolts,  or  in  hostile  environments 
by  bullets  or  even  meteorites,  holes  in  laminates  lead  to 
analysis  and  design  problems  in  composites  that  are  unique 
from  conventional  isotropic  metals.  Loads  are  not  easily 
redistributed  around  a  hole  in  the  non-homogeneous  material, 
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and  discontinuous  fibers  cause  shear  lag  effects  which  lead 
to  design  difficulties  [3]. 

The  numerous  difficulties  in  determining  the  tensile 
strength  of  notched  composite  laminates  involve  the 
consideration  of  an  infinite  number  of  parameters  which 
include : 

>1)  The  infinite  number  of  possible  ply  lay-ups,  or 
stacking  sequences, 

(2)  The  great  variety  of  fibers,  matrices,  and 
fiber/matrix  combinations, 

(3)  The  uncertain  and  inconsistency  of  failure 
cr 1  ter la , 

(4)  The  infinite  variety  of  geometries  (thicknesses, 
notch  size  and  shape,  width-to-hole  size,  etc.)  ,  and 

(5)  The  great  variety  of  environmental  conditions 
influencing  the  behavior  of  composites. 

The  techniques  available  for  the  prediction  of  strength 
reduction  or  ultimate  strength  in  a  material  with  a  hole  are 
almost  as  numerous.  In  general,  the  techniques  include 
finite  element  methods  [5,  6.  7,  8]  and  fracture  mechanics 
approaches  [3,  9,  10,  11,  121. 

Fracture  is  caused  by  higher  stresses  around  holes, 
cracks,  or  flaws  in  the  material.  Fracture  has  a  new 
dimension  in  composite  materials  because  of  the  presence  of 
t.Jo  or  more  constituents.  In  analyzing  the  failure  of 
composites,  one  can  consider  fracture  of  the  individual 
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constituents  or  the  separation  of  the  interface  between  the 
constituents  [3].  In  analyzing  the  effects  of  a  hole,  the 
classical  approaches  to  analysis  that  are  used  with  metals 
and  other  isotropic  materials  do  not  work  in  general  with  a 
composite.  The  differences  exist  because  the  stress 
concentration  factors  (SCF’s)  for  composites  have  not  been 
characterized  successfully  [13].  Isotropic  SCF’s  are  based 
on  the  conditions  at  the  hole  boundary,  while  the  strength 
of  composite  laminates  seems  related  to  the  in-plane  elastic 
stress  region  adjacent  to  the  hole  boundary.  The  measure 
of  strength  for  composites  containing  circular  disconti¬ 
nuities  must  be  based  on  a  more  complete  description  than 
simply  stress  concentration  factors  [13]. 

Early  attempts  at  predicting  the  strength  of  composites 
containing  notches  involved  the  application  of  a  linear 
fracture  mechanics  concept.  These  attempts  evolved  into  the 
two  -  par ame ter  fracture  theories  that  became  popular  in  the 
early  1970's  [8].  Some  particular  fracture  mechanics  models 
include  the  linear  elastic  fracture  mechanics  (LEFM) 
approach  [9],  the  point  stress  criterion  [9],  the  average 
stress  criterion  [9],  a  'three  parameter  model'  [13],  and  a 
f  1  r  s  t -p  1  y  -  f  ai  1  ur  e  (FPF)  model  [11].  But,  from  further 
review,  it  was  seen  that  the  fracture  mechanics  approaches 
confine  themselves  to  particular  ranges  of  parameters,  such 
as  geometry,  notch  size,  etc.  In  general,  these  theories 
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required  the  use  of  crack  length  corrections  to  reflect  the 

presence  of  damage  zones. 

Awerbuch  [10]  reviews  several  fracture  models  for 

predicting  notch  strength  of  composite  laminates,  and  notes 

that  emphasis  is  placed  on  semi -empi r i cal  techniques  which 

are  simple  to  use.  But  as  References  [5-8]  point  out, 

with  the  development  of  accurate  non-destructive 
evaluation  methods,  it  became  apparent  that  the 
damage  zone  corrections  (of  fracture  theories) 
were  arbitrary  and  subject  to  question.  As  a 
result,  attempts  were  made  to  develop  analysis 
methods  that  model  the  actual  damage  accumulation 
process  [ 8  ]  . 

Thus,  to  obtain  a  satisfactory  design  incorporating  the 
use  of  fiber-reinforced  composite  laminates,  one  must  know 
the  entire  stress -strain  behavior  to  ultimate  strength  of 
the  component.  Although  experimentation  and  repeated 
testing  will  yield  this  knowledge,  it  is  both  expensive  and 
inconvenient  considering  the  large  variety  of  parameters 
listed  above.  As  an  alternative,  another  technique  would  be 
to  analytically  relate  the  properties  of  the  individual 
plies,  or  laminas ,  to  the  laminate.  With  this  technique, 
one  would  have  to  know  the  material  biaxial  properties  of 
the  laminas  for  various  load  combinations,  which  can  easily 
be  obtained  by  experimentation  [5]. 

To  obtain  these  biaxial  properties  of  individual  lamina 
under  simple  loading  conditions  (uniaxial  tension,  uniaxial 
compression,  and  pure  shear) ,  one  may  estimate  these 
properties  from  the  properties  of  the  constituent  materials. 
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The  'Kule  zi  Mixtures'  is  a  well  known  estimation  scheme 
[3].  Fan  and  Knapp  [14]  provide  an  alternative  to  the  rule 
of  mixtures  in  the  estimate  of  composite  properties. 

Starting  with  the  Second  Law  of  Thermodynamics,  and  leading 
the  development  through  a  heat  balance  equation,  they  use 
variational  calculus  to  derive  an  ultimate  strength  equation 
that,  incorporates  nonlinearity  present  in  the  material. 
However,  this  approach  is  advocated  for  advanced  design 
purposes ;  and  verification  testing  is  still  necessary  [14]. 
These  estimation  procedures  were  not  effective  for  this 
study  so  mechanical  properties  of  the  laminas  were  obtained 
by  exper 1  men  t . 

Even  if  the  behavior  of  the  lamina  can  be  accurately 
predicted,  the  problem  of  determining  the  performance  of 
laminas  in  a  laminate  still  remains.  The  degradation  of 
strength  of  individual  laminas  may  not  precipitate  the 
failure  of  the  laminate.  For  example,  in  the  case  of  a 
cross-ply  laminate  (a  composite  laminate  with  alternating 

o  o 

layers  of  0  and  90  orientations  of  plies)  subjected  to 

normal  stresses  (o'  .o'  )  acting  in  the  direction  of  the 

12 

fibers,  the  degradation,  of  strength  of  laminas  in  the 
transverse  direction  does  not  result  in  the  failure  of  the 
laminate;  the  laminate  is  capable  of  sustaining  stresses 
higher  than  those  causing  transverse  degradation  of  the 
individual  laminas  [5].  In  general,  the  behavior  of  the 
la.'nina  under  complex  stress  states  is  predicted  by  using  a 


criterion  to  relate  the  parameters  of  strength  obtained  irom 
the  simple  tests  to  the  states  of  stress  in  the  laminate 

■  c  T 

This  lamina- laminate  response  has  been  the  subject  of 
many  theoretical  and  experimental  studies  [5].  Noyes  and 
Jones  modified  Hill’s  failure  criterion  [3,  9]  to  determine 

the  onset  of  failure  in  a  ply  or  plies.  Chiu  modified  this 
approach  to  reduce  the  affected  stiffnesses  and  strengths 
gradually,  instead  of  instantaneously  [5J. 

These  early  attempts  assumed  linearity  of  s tr ess - s t r a i n 
response,  even  though  tests  indicated  significant  nonlinear 
behavior  of  the  composite  lamina  under  transverse  and  shear 
loadings.  Petit  and  Waddoups  relaxed  this  linearity 
assumption.  They  used  properties  of  unidirectional  laminae 
in  conjunction  with  a  piece-wise  linear  approach  to 
determine  response  of  the  laminate.  I.n  their  method,  the 
loading  IS  proportional  and  is  applied  in  small  increments 
to  compute  the  current  stress-strain  states.  Biaxial 
strains  are  used  to  obtain  tangent  moduli  from  s tr es s - s tr ai n 
data  generated  under  simple  load  conditions  [6].  However, 
this  method  does  not  consider  that  strains  in  the 
3  tress -s  tram  curves  are  simple,  while  the  lamina  strains  in 
the  la.minate  are  not.  (This  aspect  will  be  discussed  in 
.more  detain  in  Section  II  F.)  Additionally,  stress-strain 
curves  of  laminates  obtained  by  this  method  are  very  much 
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dependent  upon  the  size  of  the  increment  used  in  the  applied 
loads  [ 5  ]  . 

Hahn  and  Tsai  [15]  proposed  a  method  which  allows  only 
for  the  nonlinearity  of  the  shear  s t r ess - s tr a i n  response  by 
the  e  q  ua  t,  i  o  n 

ir=So'-rS  ^  '  1  ) 

a  <s<s  ts  otstsa  a 

where  o'  and  £  are  shear  stress  and  shear  strain, 

<s  o 

resoec 1 1 ve 1 y  ,  and  S  and  S  are  constants.  This 

<S<S  <3<S<S<S 

approach  does  not  distinguish  between  tensile  and 
compressive  behaviors  of  the  lamina,  and  it  does  not  deal 
with  the  problem  of  failure  [5,  6], 

In  the  composite  industry,  at  least  30  failure  theories 
exist  for  laminated  composites  (as  of  late  1986)  [16].  Some 

o:  these  theories  are  applied  directly  to  the  laminate, 
while  the  others  are  applied  to  the  individual  laminas.  In 
addition,  at  least  12  theories  exist  of  post-failure 
behavior  of  laminated  composites  [15].  Nahas  concludes  that 
some  failure  theories  are  special  cases  of  the  maximum 
stress  theory.  He  also  found  that  some  theories  are  based 
on  assumptions  that  are  not  always  true  for  composite 
materials.  Some  of  the  more  general  theories  have  been  used 
only  for  special  cases,  which  reduce  them  to  the  Tsai-Wu 
theory,  For  theories  such  as  Sandhu's  strain  energy  failure 
criterion,  theory-experiment  correlation  coefficients  are 
included  which  make  calculated  and  test  results  agree 
rlosely.  Although  such  coefficients  must  be  determined  for 
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each  material  [16j,  once  they  are  obtained,  the  finite 
element  method  using  Sandhu ' s  criteria  may  be  used  for  a 
variety  of  geometries.  Rowlands  [17]  remarks  in  his  survey 

in  the  Handbook  o  f  Compos i tes  that 

whereas  most  of  the  strength  predictions  assume  a 
linear  stress-strain  relationship  in  the  associated 
lamination  analysis,  those  by  Sandhu  and  Petit- 
Waddoups  provide  for  nonlinear  lamina  material 
response  [17], 

The  composite  strength  analysis  proposed  bv  Sandhu 
differs  from  the  more  conventional  formulations  in  the 
following  respects  [17]; 

(1)  The  nonlinear  lamina  stress-strain  responses  to 
failure  are  represented  analytically  by  cubic-spline 
functions.  The  tangent  moduli  of  these  functions  are 
employed  to  evaluate  lamina  and  laminate  stiffnesses  and 
compliances  during  load  increments. 

(2)  In  keeping  with  the  nonlinear  provision,  ply 
degradation  is  based  on  an  energy  to  failure  criterion. 

(3)  Equivalent  strain  increments  are  defined.  (See 
Section  II  F.) 

Thus  considering  the  above  justifications,  this  study 
uses  Sandhu's  fully  nonlinear  progr ess i ve - p 1 y - f a i 1 ur e  finite 
element  .method  to  analyze  the  damage  accumulation  process  in 
jr^PEEK  tension  specimens  containing  a  circular  disconti¬ 
nuity.  This  study  uses  the  technique  of  Sandhu  for 
predicting  the  response  to  failure.  Previous  works  have 
describe'!  this  technique  in  detail  [5,  6],  and  others  have 
used  this  technique  for  predicting  behavior  in  Gr/Ep  [5.  7, 


3,  18,  19],  In  applying  this  technique  to  predict  the 

behavior  of  Gr/PEEK,  no  previous  work  could  be  found. 

As  a  general  description  of  this  analytical  technique, 
cubic  spline  interpolation  functions  are  used  to  represent 
basic  s tress -s trai n  curves  obtained  from  simple  tests  [5j: 

1.  Longitudinal  tension  and  compression. 

2.  Transverse  tension  and  compression. 

3 .  Shear . 

Cubic  spline  functions  yield  smooth  composite  stress-strain 
curves  (from  sets  of  three  data  points)  from  which  the 
computer  program  can  determine  accurate  moduli  of  elasticity 
over  the  entire  range  of  the  curves  [6]. 

This  functional  form  of  stress-strain  curves  provides 
an  accurate  representation  of  stresses,  strains,  and  moduli 
over  the  entire  range  of  curves.  This  representation  is 
used  in  conjunction  with  an  incremental  constitutive  law. 
This  law,  which  relates  the  increments  of  stresses  and 
strains,  is  used  to  generate  the  response  of  Gr/PEEK 
laminates  subjected  to  tensile  loads.  The  basic  data  and 
constitutive  relations  are  then  used  to  determine  the 
ultimate  1 oad- carry i ng  capacity  of  the  laminates  under 
incremental  loading  by  a  plywise  app  1 1 c;<.t i on  of  Sandhu’s 
failure  criterion  using  strain  energy  under  longitudinal, 
transverse,  and  shear  loadings  as  independent  parameters 
[5.61.  The  theory  behind  this  technique  is  discussed  in 
detail  in  Chapter  II.  and  the  procedure  used  in  this  thesis 
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for  accomplishing  the  analysis  is  described  in  Chapter  III. 
The  results  obtained  by  this  technique  are  compared  with 
experimental  results  of  this  study. 

The  laminates  studied  through  both  analysis  and 

experimentation  were  the  [0  ]  and  [90  ]  unidirectional 

1  <s  ta 

lay-ups  and  the  angle-ply  lay-up.  The  two 

quas 1  -  1  so tr op  1 c  lay-ups,  [ 0/ +45/90/ -45  ]  and  [0/i45/90] 

2S  2S 

were  only  studied  experimentally  for  comparison  to  similar 
Gr^^Ep  laminates.  All  test  specimens  contained  16  layers  and 
were  symmetric  about  their  mid-plane.  Also,  they  all 
contained  a  0.4'  diameter  hole  at  their  centers.  Thus  the 
laminate  geometry  represented  a  finite-width  plate 
containing  a  stress  concentrator  subjected  to  uniaxial 
tension.  Furthermore,  bending  was  not  induced  in  the 
specimens  by  tension  because  of  the  symmetry  of  the  lay-ups. 

Basic  properties  were  first  derived  with  the  experi¬ 
mentation  process,  and  then  the  five  types  of  laminates  with 
holes  were  tested  to  their  ultimate  loads  and  to  percentages 
of  ultimate  load.  Testing  to  ultimate  yielded  the  laminate 
stress-strain  responses  for  the  entire  range  of  strength. 
Testing  to  percentages  of  failure  load  yielded  partially 
failed  specimens  which  could  be  examined, by  stereo  x-ray. 

In  addition  to  the  x-rays,  one  other  medium  allowed  for 
a  post-failure  analysis.  A  video  camera  filmed  each  test  to 
ultimate  strength.  Through  repeated  viewing  of  the  tape  and 
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through  stills,  or  photographs,  taken  oft  the  tape,  one 
could  deduce  further  results  from  the  experimentation. 

A  description  of  the  procedures  used  in  the 
experimentation  phase  of  this  thesis  is  presented  in 
Chapter  IV.  Results  of  both  experimentation  and  analysis 
are  presented  and  discussed  in  Chapter  V.  Conclusions  on 
the  entire  study  are  reported  in  Chapter  VI. 

C .  Material  Choice 

A  poiyaromatic  resin,  such  as  po 1  ye  there therke tone 
(PEEK) .  combined  with  a  continuous  fiber  is  known  as  an 
aromatic  polymer  composite  (AFC)  [13].  The  Gr/PEEK  chosen 
for  this  study  is  refered  to  as  APC-2.  an  improved  laminate 
over  APC-1.  Whereas  APC-1  contained  Courtauld  XAS  graphite 
fibers,  APC-2  contains  Hercules  AS-4  graphite  fibers.  The 
AS-4  fibers  make  a  better  bond  with  the  PEEK  matrix  [13]. 
Gr/PEEK  can  be  converted  into  a  range  of  component  shapes 
and  sizes  by  a  full  spectrum  of  fabrication  technologies 
[ 2 .  4  ,  13  ]  . 

Demuts  and  Sharpe  [1]  found  APC-2  to  be  257.  tougher 
than  graphi te/epoxy  (Qr/Ep)  in  static  post-impact  residual 
compression  strength.  They  further  found  that  the  fatigue 
life  of  thermoplastic  laminates  is  twice  that  of  comparable 
thermosets.  They  conclude  that  thermoplastic  matrix 
composites  have  the  potential  for  more  efficient  wing  skins 
at  lower  fabrication  and  maintenance  costs  as  compared  to 
the  more  brittle  thermosets.  Carpenter  [20]  reiterates 
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these  findings,  praising  the  potential  of  thermoplastics  for 
lowering  production  costs  and  improving  maintainability. 

Ramey  [13]  conducted  a  comparison  study  of  strength 
reduction  in  APC-2,  APC-l,  and  Gr/Ep  using  tension  and 
compression  specimens  with  circular  holes  at  their  centers. 
He  found  APC-2  to  be  superior  in  its  strength  properties 
over  the  others.  His  comparisons  included  microphotographs 
of  each  failed  specimen,  and  an  analysis  using  a  fracture 
mechanics  approach.  He  found  that  the  failure  surfaces  o: 
APC-2  were  combinations  of  different  failure  modes,  and  he 
concluded  that  there  was  ‘no  general  failure  prediction 
model’  that  covered  a  range  of  temperatures  or  a  broad  range 
of  hole  sizes.  He  also  suggested  that  a  continuum  mechanics 
approach  would  be  advisable  in  analyzing  the  behavior  of 
APC-2  [13].  Thus,  the  logic  is  reinforced  of  using  the 
pr 0 gr ess  1 ve - p 1 y - f ai 1 ur e  finite  element  analysis  in  this 
study.  The  strain  energy  failure  criterion  used  in  this 
analysis  exhibits  the  characteristics  of  a  continuum 
mechanics  approach. 

D .  Progress  1  on  o  f  Study 

This  thesis  began  with  a  simple  abstract  and  test  plan, 
and  evolved  into  producing  further  information  about  Gr/PEEK 
and  Its  nonlinear  behavior.  The  planning  of  this  thesis 
required  a  large  logistic  effort  on  a  local  level.  Several 
organizations  at  Wr i gh t - Pat terson  AFB  aided  in  the 
completion  of  this  thesis.  The  process  of  seeking  the  help 
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of  these  organ i za t i ons  and  communicating  to  them  a  research 
problem  had  not  been  documented.  For  completeness,  the 
logistics  efforts  of  this  thesis  are  included  in  Appendix  A 
The  author  hopes  that  this  appendix  may  help  streamline  the 
research  efforts  of  future  AFIT  students. 


T  h  e  o  r  V 


:  I . 

This  chapter  discusses  theories  used  to  describe 
stresses,  strains,  and  failure  in  f i ber / po i ymer  composite 
.materials.  First  considered  is  the  behavior  of  the 
romrosite  on  a  mi cr omechan i ca 1  level,  and  then  the 
discussion  turns  to  one  of  macromechanics.  Mi cromechan i cs 
examines  the  interaction  of  the  constituent  materials, 
defining  the  behavior  of  the  heteroi^eneous  composite 
material.  But  in  an  engineering  approach  to  composite 
materials,  one  chooses  to  view  the  entire  ma t e r i a 1 
macroscop  1  cal  1 y .  Macromechanics  examines  the  effects  of  the 
constituent  materials  as  averaged  apparent  properties, 
defining  the  behavior  of  the  homogeneous  composite  material. 

Following  this  review  of  mechanics  of  composite 
materials,  this  chapter  presents  a  brief  formulation  of 
.inear  finite  element  theory,  which  is  then  expanded  to 
handle  .mul  1 1  p  1  e  -  1  ay  er  ed  composite  materials.  Further 
development  leads  to  nonlinear  finite  element  methods. 
Finally,  the  failure  criteria  used  to  match  analytical 
results  to  experimental  findings  concludes  this  chapter. 

A .  Micro mechanical  Behav i or  o  f  Compos i te.  Material . 

Although  a  designer  is  generally  interested  in  the 
resultant  properties  of  the  fiber/matrix  combination,  it  is 
often  useful  to  know  the  micromechanical  behavior  of  the 
.material  to  determine  optimum  resultant  properties.  For 
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example,  toughness  or  resistance  to  temperature  and  humidity 
can  be  maximized  in  a  material  by  starting  at  the  molecular 
level  and  designing  outward.  To  better  understand  the 
failure  of  the  composite,  the  specific  composition  of  each 
constituent  must  be  known.  This  section  discusses  the 
nature  of  polymer  matrices  and  organic  fibers. 

Polymers  are  formed  as  a  collection  of  long-chain 
molecules.  Elemental  units,  or  monomers,  linked  together 
farm  the  long-chain  molecules.  The  units  in  a  chain  can 
vary  in  element  and  size,  yielding  an  infinite  number  of 
possible  polymers  [13].  Organic  or  carbon-based  polymers 
have  chains  based  on  a  unit  of  carbon.  The  PEEK  polymer  is 
formed  from  the  carbon-based  repeat  unit  shown  in  Figure  1 

r  ->  ' 

The  polymer  chains  can  take  many  shapes.  In  rubber, 
for  example,  the  carbon-based  units  form  a  spiral,  thus 
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In  other  ma  ter: ala 


yie*din§  a  highly  deformable  ma t e r i a  1 . 
carbon  chain  may  be  a  zig-zag  and  yield  a  material  which  : 
lesa  ieftrmable. 

When  units  in  the  main  chain  are  linked  to  other 
similar  chains  by  side  chains,  a  network  polymer  is  formed 
This  dense  cross-linked  polymer  is  referred  to  as  thermo¬ 
setting  because  it  does  not  soften  upon  heating  once  it  is 
formed.  A  thermosetting  material  re.mains  hard  and 
eventually  d i s i n t e g r a t e s  under  extreme  temperatures  '12!. 
The  matrix  epoxy  of  Gr/Ep  is  formed  from  this  type  :: 
ma  t  e  r  1  a  1  . 

Linear  polymers  with  small  to  no  amount  of  cross-link 
form  VISCOUS  fluids.  These  fluids  exist  as  tang.ed  .masses 
of  molecular  segments  sliding  over  one  another.  The 
viscosity  of  these  fluids  varies  greatly  with  temperature. 
These  thermally  activated  fluids  are  called  ther mo p 1  as 1 1 c 
polymers.  Below  their  glass  transition  temperature,  T^ , 
they  are  typical  glass.  Above  their  T  .  they  range  from 

O 

fluids  to  deformable  solids,  depending  on  their  amount  of 
oross-li.nks  bet  we  en  molecular  chains  [13!.  Polyetherether 
ketone  'PEEK'  is  a  thermoplastic  .material  e.xhibiting  these 
p  r  -  p  e  r  1 1  e  s  . 

» 

The  strain  in  a  thermoplastic,  as  well  as  in  any 
material .  is  a  combination  of  elastic  strain  plus  creep,  o 
t *  me -de penden t  strain.  The  rate  of  chain  straightening 


.  ension  strain)  or  crumbling  i compression  strain!  depends 
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on  viscosity.  When  a  material  is  unloaded,  elastic  strain 
IS  recovered  instantly;  but  strain  due  to  creep  recovers 
gradually  because  of  the  viscous  drag  on  the  molecular 
chains  [13j.  According  to  the  creep  data  provided  in 
Reference  [2],  for  the  duration  of  each  test,  creep  is 
negligible  for  each  type  of  ply  arrangement.  Therefore,  one 
assumption  of  this  study  is  that  creep,  or  v i sco -e 1  as 1 1 c i ty  . 
ran  be  ignored. 

In  the  manufacturing  of  thermoplastics,  one  must  con¬ 
trol  the  crystallinity  of  PEEK  to  optimize  the  performance 
of  the  polymer.  This  is  done  by  controlling  the  rate  at 
which  the  composite  is  cooled  down  from  its  melted  state  at 
"’20  F-750  F  to  Its  T^  at  300  *F  [2].  In  crystalline 
regions,  the  chains  are  parallel;  in  non-crystal  1  me 
regions  they  are  tangled.  A  single  molecular  chain  may 
exist  in  both  crystalline  and  non-cry s ta 1 1 i ne  regions  [13]. 
If  APC-2  IS  cooled  at  rates  less  than  18  F  per  minute, 
cr ys ta 1 1 1 n 1 ty  will  increase  beyond  optimum,  which  will 
result  in  reduction  of  composite  toughness.  If  cooling 
rates  exceed  1 260  "F/minute  ,  the  polymer’s'  ‘spherulitic 
growth'  will  not  achieve  completion  and  the  optimum  level  of 
rrys tal 1 1 n 1 ty  will  not  be  reached.  This  may  cause  reduction 
in  stiffness  and  resistance  to  solvents’  [2].  Reheating  and 
recooling  (annealing)  the  composite  can  correct  any  errors 
in  cooling  so  that  optimum  crystallinity  can  be  achieved. 
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In  APC-2,  the  PEEK  matrix  surrounds  and  chemically 
bonds  to  AS4  carbon  fibers  to  form  the  composite  material. 

As  alluded  to  in  Chapter  I,  this  combination  of  fiber  and 
matrix  form  a  thin  lamina  when  continuous  fibers  are  laid  in 
a  single  direction  or  woven  in  various  directions.  Figure  2 
shows  two  principal  types  of  laminae  C3]. 

The  fibers  are  the  principal  reinforcing  or  load- 
carrying  agent.  Advanced  carbon  fibers,  such  as  AS4  fibers, 
are  produced  from  precursor  fibers  in  three  steps.  The 
first  step  involves  controlled  oxidation  of  carbon  fila¬ 
ments.  Complex  reactions  between  the  aligned  polymer  chains 
of  the  precursor  form  crosslinks,  and  an  extended  carbon 


LAMINA  WITH 
UNIDIRECTIONAL  FIBERS 


LAMINA  WITH 
WOVEN  FIBERS 


(:a.)  cb.) 


Figure  2.  Two  Principal  Types  of  Laminas 
Note:  APC-2  laminates  follow  the  pattern  of  Figure  (a) . 
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network  is  formed  which  is  aligned  along  the  fiber  axis 
:21].  Individual  filaments  range  in  diameter  from  six  to 
nine  micrometers  [2].  In  the  second  step  pyrolysis  under 
an  inert  atmosphere  drives  off  most  non-carbonaceous  atoms. 
Finally,  the  third  step  exposes  the  carbon  filament  to  high 
temperature  to  form  extended  graphite  ribbons  aligned  along 
a  fiber  axis  [21],  In  this  graphi t i zat i on  of  carbon,  a 
layered  hexagonal  crystal  structure  is  formed  under  heat. 
Graphite  fibers  may  also  be  subjected  to  surface  treatment 
to  improve  their  compatibility  with  the  matrix.  [21]. 

With  the  fibers  carrying  the  load,  the  matrix  acts  to 
support  and  protect  the  fibers  and  to  provide  a  way  of 
distributing  the  load  among  the  fibers  and  transmitting  the 
load  between  fibers.  As  shown  in  Figure  3  [3],  this 
distribution  and  transmission  of  load  is  especially 
important  if  a  fiber  breaks.  The  load  is  transferred  to  the 
matrix  and  then  to  the  other  side  of  the  broken  fiber,  as 
well  as  to  adjacent  fibers.  The  shearing  stress  developed 
in  the  matrix  is  the  mechanism  for  the  load  transfer.  This 
shearing  stress  resists  the  pulling  out  of  the  broken  fiber 
[3].  This  description  of  mechanical  action  between  fiber 
and  matrix  accounts  for  much  of  the  explanation  of  failure 
progression  in  composite  materials. 

To  determine  the  stiffness  properties  of  composite 
materials,  many  analytical  approaches  are  possible.  One  can 
use  the  mechanics  of  materials  approach,  or  the  Rule  of 


n  n 


in  which  3 1  c\p  1  e  approximations  to  the  engineering 
constants  are  derived  by  averaging  the  properties  of  each 
constituent  material  over  the  composite  cross-section. 
Another  approach  is  that  of  elasticity.  This  approach 
includes  bounding  techniques,  exact  solutions,  and  concepts 
of  contiguity. 

Kalpin  and  Tsai  developed  an  interpolation  procedure 
that  IS  an  approx i ma te  representation  of  more  co mp 1 i c a t e d 
mi cr cmechan 1 cs  results  13].  According  to  Jones,  Tsai  gives 

expressions  for  E  ,  E  ,  .  and  G  "  that  are  in  good 

1  2  12  12 

agreement  with  experimental  data  for  a  glass  reinforced 
resin  composite,  but  a  'contiguity  factor'  [3]  is  the  key 
to  the  agreement.  Tsai  found  that  the  constituent  material 
properties  have  the  following  effect  on  the  properties  of 
the  composite: 

Young’s  modulus,  E,  of  the  fiber  makes  a 
significant  contribution  to  E^  of  the  lamina, 

(2!  E  of  the  matrix  makes  a  significant  contribution 

to  E  and  G  of  the  lamina,  and 

2  12 

i3)  The  Poisson  ratios  for  both  the  fiber  and  the 

.ma  t  r  1  ;<  have  little  effect  on  E  and  G  and  have  no  effect  on 

2  12 

E  of  the  lamina. 

1 

Although  engineering  properties  could  be  derived  with 
these  mi  cr o mec han 1 ca 1  approaches,  this  study  conducted 
e  ;-:per  1  ,men  1 3  to  find  these  properties.  Not  only  was  accuracy 


i'93ire<i,  but  nowhere  did  these  methods  accurately  and  easil 
take  into  account  the  nonlinearity  of  material  properties. 

3.  Macromechanical  Behav i or  o f  C o mp o s 1 1 e  Ma t e r i a 1 s 

A  lamina  is  the  basic  building  block  in  a  laminated 
f 1 ber - r e 1 n f or ced  composite.  In  this  section,  tho  behavior 
of  a  lamina  and  a  laminate  will  be  considered  from  a 
.macromechanical  approach.  A  laminate  is  two  or  more 
laminas  bonded  together  to  act  as  an  integral  structural 
element.  This  stack  of  laminas  may  have  plies  of  various 
orientations  of  principal  mater lal  directions  as  shown  in 
Figure  4  13].  Note  that  the  laminate  is  shown  under 

uniaxial  loading,  and  the  orientations  of  the  plies  are 
measured  with  respect  to  this  axis  of  loading. 

F 1  be r - r e 1 n f or ced  composites  are  usually  treated  as 
linear  elastic  materials  since  the  fibers  provide  the 
majority  of  the  strength  and  stiffness.  This  thesis 
attempts  to  refine  that  approximation  by  considering 
nonlinearity.  Although  most  theories  are  based  on  linear 
assumptions,  repeated  exper  i  .men  ta  1 1  on  has  shown  that 
f 1 ber - r e 1 n f or ced  composites  exhibit  nonlinearity,  and  the 
e xper  1  .me n t a 1 1  on  in  this  thesis  also  displays  this 
phenomenon.  (.See  Chapter  V.  Results  and  Discussion.!  To 
treat  the  nonlinearity  problem,  the  linear  elastic  theories 
presented  in  this  section  have  been  employed  such  that  if  a 
given  load  increment  is  small  enough,  linear  elastic  theory 
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applies.  As  will  be  shown  in  Sections  F  and  G  of  this 
chapter,  linear  elastic  theory  is  applied  incrementally 
while  the  basic  properties  of  the  material  vary  with  each 
increment . 

A  potential  problem  in  the  construction  of  laminates  is 
the  introduction  of  shearing  stresses  between  layers.  The 
shearing  stresses  arise  due  to  the  tendency  of  each  layer  to 
deform  independently  of  its  neighbors  because  all  may  have 
different  orientations  and  thus  different  properties  in  the 
direction  of  loading.  Such  shearing  stresses  are  largest  at 
the  edges  of  a  laminate  and  may  cause  delamination  there. 

The  transverse  and  normal  stresses  resulting  from  uniaxial 
loading  can  also  cause  delamination  [3]. 

Whitney  discusses  this  subject  of  interlaminar  stresses 
in  Reference  [9],  and  he  discusses  the  influence  of  stacking 
sequence  on  interlaminar  stresses.  The  most  influential  of 
the  interlaminar  stresses  is  the  normal  stress,  i.e. ,  the 
stress  attempting  to  separate  the  plies.  Sandhu  [22,  23] 

describes  the  tendency  for  a  composite  to  delaminate  through 
the  use  of  a  'delamination  moment  coefficient'  (DMC) ,  which 
13  based  on  the  stacking  sequence.  For  Gr/Ep  specimens 
identical  in  geometry  to  those  of  this  thesis,  the  critical 
DMC  IS  approached  for  the  [0/±45/90]  laminate  under  tensile 
loading.  A  ply-by-ply  plane  stress  analysis,  such  as  the 
one  used  in  this  thesis  (see  Chapter  III) ,  does  not  take 


into  consideration  stacking  sequence. 


Thus,  reductions  in 


strength  caused  by  de 1 aminat i ons  would  cause  unwanted 
disagreement  between  analysis  and  experimentation. 

Therefore,  the  [ 0/ »45/90/-45  ]  laminate  was  also  chosen  to 
evaluate  in  this  study  for  comparison  to  the  [0/t45/90] 
laminates . 

In  describing  the  behavior  of  a  single  lamina,  the 
basic  restriction  of  the  macromechanical  approach  is  to 
linear  elastic  behavior.  This  assumption  is  appropriate  in 
this  thesis  because  nonlinear  behavior  is  modeled  as  a 
series  of  increments  which  deduce  the  varying  material 
properties  from  smooth  curves  of  stress-strain  data. 
Summarizing  the  linear  elastic  theory  for  a  lamina,  the 
generalized  Hooke's  law  relating  stresses  to  strains  can  be 

written  in  contracted  notation  as 

* 

cr  -  C  £  i,j  =  l . 6  (2) 

>-1  J 

where  cr  are  the  stress  components.  C  is  the  stiffness 

matrix,  and  are  the  strain  components  [3].  The 

components  of  the  stiffness  matrix,  and  the  terms  of  its 

inverse  the  compliance  matrix,  S  .  will  be  referred  to  as 

'•j 

elastic  constants  (keeping  in  mind  the  incremental 
assumption  of  the  above  paragraph) . 

A  unidirectional  reinforced  lamina^ in  the  macro- 
.mechanical  approach  is  considered  orthotropic,  i.e,  ,  there 
are  three  mutually  orthogonal  planes  of  material  property 
symmetry.  These  planes  of  symmetry  are  described  by  the  set 
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Figure  5.  The  Principal  Material  Axes  for  a 
Unidirectional ly  Reinforced  Lamina 


of  axes  in  shown  above  in  Figure  5  [3]  as  the  principal 
material  directions. 

Taking  the  inverse  of  Hooke’s  law  above,  the 
strain-stress  relations  become 

^=3(7  i.j  =  l . 6  (3) 

'■  VJ  J 

For  an  orthotropic  material  with  body  coordinate  axes  X  and 
Y  aligned  with  principal  material  directions,  Eq .  2  becomes 
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wh  ere  }■  =  2  £  The  ter  ms  j'  (i^j)  represent 

L)  I.J 

engineering  shear  strain,  whereas  ( i^j )  represent  tensor 

shear  strain  [ 3  j  . 

For  the  study  conducted  in  this  thesis,  the  tension 
specimens  representing  finite-width  thin  plates  were 
subjected  to  plane  stress.  For  a  lamina  in  the  1-2  plane,  a 
state  of  plane  stress  is  defined  by  setting 

<o=i3,r=0,r  =0  (5) 

a  23  13 

in  the  three  di me nsional  strain-stress  relation  given  in 
Fq . i4'  above.  Thus,  the  s tr a i n - s t r es s  relations  in  Eq .  i4) 

above  reduce  further  to 
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where  the  compliance  terms  are: 
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s  =  s  =  s  =  s 

la  2a  ai  a2 

The  stain-stress  relations  in  Eq . 
obtain  the  stress-strain  relations: 


=  0 


(6)  are  inverted  to 
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where  the  are  called  the  ‘reduced  stiffnesses,"  are 


defined  as  [ 3 ] : 
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Q  =  G 
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<3  =Q  =0  =0  =0 

1<3  ''2<5  ^<S1  ^ 


The  preceding  s tress -s tra i n  and  strain-stress  relation 
are  the  basis  for  the  stiffness  and  stress  analysis  of  an 
individual  lamina  subjected  to  forces  in  its  own 
p*ane.  Note  that  the  stiffness  and  compliance  terms  are 
based  on  four  general  independent  material  properties: 

E,  .  E  .  u  .  and  G 

i  2  12  12 

Where  subscripts  1  and  2  refer  to  principal  material  axes. 
Mote  aiso  that  these  four  properties  are  expanded  to  seven 
specific  properties  when  considering  tension  and 
rompression.  These  properties  are  designated  as; 

T  C 

^  i  '  ^  I  '  ^  2  •  ^  2  '  -• 
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where  f^he  superscripts 


and  C  represent 


ensicn  and 


comDression  basic  properties.  The  shear  modulus  G  is  not 

12 

affected  by  the  change  of  unidirectional  load. 

Often  the  principal  directions  of  a  lamina  do  not 
coincide  with  the  coordinate  directions  of  the  component, 
the  X  and  '{  directions,  that  is.  axes  1  and  2  are  oriented 
at  an  angle  S  from  the  X  and  Y  axes.  This  rotation  of 
pri.ncipal  material  axes  is  shown  in  Figure  6  3  ]  .  To 

handle  this  problem,  first  a  transformation  relation  is 
needed  between  the  stresses  in  the  principal  material 
directions  and  those  in  the  body  coordinates.  Second  .  a 
.method  is  needed  to  transform  s  tress  -  s  t^a  i  n  relations  from 
one  coordinate  system  to  another. 

Through  this  procedure  of  transformations,  the 
3 t r es 3 - s t r a  1 n  relation  of  Eq .  (8;  can  be  written  that  will 

account  for  any  orientation  of  the  lamina  [3]: 


f  1 

X 

c 

X 

r  ^ 

11 

^12 

®1<5  ' 

€ 

X 

a 

y 

c 

y 

► 

«12 

«22 

^2. 

€ 

V 

Y 

• 

«2a 

V 

* 

L 

where  the  terms  Q  denote  transformed  reduced  stiffnesses, 
as  opocsed  to  the  reduced  stiffnesses  of  Q 

'-j 

*  are  given  as  131: 


The  values  o 


*  2(Q+2Q  )n^m^  +  Qn* 

11  11  12  <S<$  22 


^12  "  ^'^11  *  ^22'“^  ®aa' 


Q  =  Q  n 

22  11 


(  Q  +  2  Q  )  +  Q  m* 

12  aa  22 


^1<S  =  ^®11  ■  ®12  ■  2  ^  ®12'  ®22"  2 

^2<S  =  ‘«11  ■  «12  •  -  "  *  «12-  «  *  2  Q^)n  m" 

'2a<5  =  '®ii  "  «22  ’  2  -  2  nV  *  Q^tn*  * 


wh  ere  the  values  of  Q  we  re  ;?iven  in  Eq.  (9).  Note  that 

'■j 

there  are  still  only  four  general  independent  constants 
which  characterize  an  orthotropic  lamina. 

Alternatively,  the  strains  can  be  expressed  in  terms  o 
the  stresses  in  body  coordinates  by  inversion  of  Eq  .  (10)  , 


S  S  S 

11  12  1<$ 

s  s  s 

12  22  2<S 

s  §  § 

1<S  2(3  <3<3 


where  values  of  the  transformed  reduced  compliance  terms, 
S  ,  are  given  as  [3]: 


.  2  2  ,  <4 

S*S  )n  m  *bn 

12  <3<3  22 


,22 

n  +  m''‘(S*b  -o  )n  m 

12  (3<3  <3<3 


^  22,4 

on  •►f2S»S  ,'nm  *bm 

11  12  (va  22 


-  S  !  n  m 

12  <S<5 
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=  (2  S  -  2  S  -  S  )n  m  -  (2  S  -  2  S  -  S  )  n  m' 

2<S  11  12  <5<S  22  12  <5<3 

S  =2(2  S  *  2  S  -  4  S  -  S  )  -<•  S  (n'*+  m* . 

CXS  11  2  2  12  <S<S  <5<S 


where  the  values  of  S  were  given  in  Eq .  (7) 


■-j 


2 .  Strength  o  f  Composite  Lami nates 

In  describing  strength  of  composite  materials,  deter¬ 
mination  of  the  strength  of  laminate  is  built  upon  the 
strengths  of  its  laminas.  Because  of  the  anisotropic 
and  hetergeneous  nature  of  composite  materials,  failure 
modes  are  often  quite  different  from  those  of  isotropic 
homogeneous  materials.  In  particular,  the  failure  of  one 
layer  in  a  composite  material  does  not  necessarily  mean  the 
entire  laminate  will  fail.  The  laminate  may  be  capable  of 
sustaining  higher  loads  despite  a  significant  change  in 
st.ffness  13].  This  concept  was  introduced  in  Chapter  I  and 
•  s  shown  graphically  in  Figure  1]3]. 

Because  of  the  various  characteristics  of  composite 
materials,  it  is  difficult  to  determine  a  strength  theory  in 
which  all  failure  modes  and  their  interactions  are  taken 
int.3  account.  Strength  theories  presented  in  Reference  [3] 
are  based  on  a  macroscopic  analysis  in  which  the  strengths 
if  each  lamina  must  be  assessed  as  it  relates  to  the  whole 
.aminate.  They  include  maximum  stress,  maximum  strain,  and 
maximum  listortional  energy  approaches.  As  described  in  the 
;  'll  owing  sections,  this  thesis  uses  a  finite  ele  me  n  t 
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approach  to  the  overall  failure  of  the  laminate,  but  its 
inherent  strain  energy  failure  criterion  for  individual 
elements  is  not  unlike  the  other  analytical  approaches 
me  ntioned  above.  See  Section  II  F. 

In  any  method  of  strength  analysis,  the  designer  is 
striving  to  determine  either  (1)  the  maximum  loads  a  given 
laminate  can  withstand  or  (2)  the  laminate  characteristics 
necessarv  tc  withstand  a  given  load.  [3].  The  overall 
procedure  of  determining  laminate  strength  is  straight- 
fcrward,  but  extremely  tedious;  thereby  implying  the 
need  for  computer  techniques.  A  general  procedure  is  shown 
schematically  in  Figure  8  [3],  and  the  computer  application 
used  in  this  thesis  is  a  specific  version  of  this  general 
cr?c®dur9 • 

I.  .  H r  1  e s  :  n  Lami  nates 

For  the  problems  of  holes  in  laminates,  one  of  the 
first  solutions  to  the  difficult  problem  of  stresses  around 
a  .-.ole  was  given  by  Lekhnitskii  in  1936.  1963,  and  1958. 

13;.  jreszczuk  plotted  the  circumferential  stress  around 
the  hole  for  an  isotropic  material  and  several  unidirec- 
t.onal  materials.  See  Figure  9  [3].  Tlie  usual  stress 

c :  .ncentr  at  1  on  for  an  isotropic  material  is  3,  i.e., 

o 

-1'  o ,  =  3  a  t  "9  =  9  0  .  For  composite  materials  such  as 

grash ite  eooxy,  the  stress  concentration  factor  depends  on 
.  a  m  1  a  t  e  construction  . 


•  LAMINAE  PROPERTIES 

•  PROPORTIONS  BETWEEN  APPLIED  LOADS 
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The  circumferential  stress  at  ©  =  0 ”  is  reduced  for 
composites  relative  to  isotropic  materials.  Because  of 
isotropy  of  material  properties,  the  key  factor  in  failure 
of  isotropic  plates  with  holes  is  the  magnitude  of  the 
stress  concentration  factor  from  which  the  maximum  (failure) 
stress  IS  obtained.  However,  for  orthotropic  materials,  a 
combined  stress  failure  criterion  instead  of  a  maximum 
stress  failure  criterion  is  more  accurate  since  failure  of  a 
lamina  is  a  function  of  the  strengths  in  various  directions. 
Thus,  stress  concentration  factors  alone  are  insufficient 
for  failure  prediction  of  orthotropic  (and  anisotropic) 
plates.  Moreoever ,  for  laminated  plates,  the  comparison  of 
stress  states  with  failure  stress  states  must  be  done  on  a 
layer-by-layer  basis  [3]. 

Extending  the  analysis  to  a  laminate.  Figure  10  [3j 
shows  stress  concentrations  around  a  circular  hole  in  a 
cross-ply  and  an  angle-ply  laminate.  The  stress  shown  is  a 
gross  stress  on  the  laminate,  and  the  stresses  in  each  layer 
must  be  found  by  classical  lamination  theory  or  some  other 
method.  Failure  is  determined  by  application  of  failure 
criterion  to  each  layer.  The  interlaminar  stresses  are 
Ignored,  and  thus  predicted  stresses  are  not  accurate  within 
about  one  la.minate  thickness  from  the  edge. 

These  first  sections  have  presented  an  overview  of  some 
theories  and  approaches  used  in  studying  the  mechanics  of 
composite  materials.  The  next  three  sections  describe  the 
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specific  theories  behind  the  analysis  accomplished  in  this 
thesis  . 

E  .  L.  1  near  Finite  Element  Theory 

For  the  convergence  study  described  in  Section  III  C, 
this  study  used  a  linear  finite  element  program  called 
PLSTR ,  written  by  Dr.  R.S.  Sandhu,  that  is  a  s i mp 1 i f i ca i on 
of  his  nonlinear  program  PLSTREN ,  which  is  described  in 
detail  in  References  [5]  and  [6],  Both  programs  are  a 
ply-by-ply  analysis  of  finite  ele me n t  me shes  representing 
the  laminas  of  a  composite  material.  The  program  takes  into 
account  the  number  and  orientation  of  orthotropic  plies  in 
the  model,  and  it  can  model  an  infinite  number  of  geometries 
and  discontinuities.  This  section  describes  the  theory 
behind  the  analysis,  beginning  with  linear  finite  element 
theory  and  extending  the  theory  for  use  with  multiple- 
1  avered  .meshes  . 

PL-PTR  13  based  on  the  constant  strain  triangle  as  shown 
in  Figure  11,  which  implies  a  two-dimensional  elasticity 
problem.  Quadrilateral  elements  in  this  program  are 
assemblages  of  four  constant  strain  triangles,  which  then 
resemble  linear  elements  [24],  The  assumed  displacement 
field  13  linear  in  the  x  and  y  directions: 


4  2 


u  (  X  ,  y  t 
V  ( X  ,  y ) 


0  0  0 


0  0  0  1 


’hus  £  =  a  ,  c  =  a  ,  and  y 

X  2  y  <S  12 


a  --a  L  24  ]  .  The 

3  5 


constant  strain  triangle  assumes  no  variance  of  strain 
through  the  thickness  of  the  element,  i.e. .  no  bending. 

Through  operations  on  the  shape  factors  (See  Reference 
'24]) ,  the  strains  are  expressed  by. 


=  r  B  ](  d  )  =  [  3  ] 


where  (  d  )  is  the  nodal  displacement  vector  and, 


(  y  -  y  ) 


( y  -y  : 

3  1 


.7,-72- 


(  X  -X  ) 

3  2 


:  X  -  X  ) 

1  3 


(X  -X  1 

2  1 


(x-x)  (y-y)  (x-x)  (y-y!  (x-x)  (v-v  ! 

3  2  ’^2'^3  1  3  31  2  1  12 


where  A  is  the  area  of  the  triangle  and  x  and  v  are  the  x 

I  '  I 

and  y  coordinates  of  the  i node  [  1 '3  ]  . 
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Developing  an  equation  for  the  stiffness  of  an  elemen*  . 
the  differential  volume  of  an  element  for  constant  thickness 
IS  dV  =  t*dA,  where  t  is  the  element  thickness.  (Note  chat 
this  thickness  is  that  of  a  ply.)  The  matrix  [B]  is 
constant,  and  the  material  property  matrix,  defined  in 
general  as  [D],  is  assumed  constant  over  the  element. 

The  above  assumptions  yield  the  equation  for  the 
element  stiffness  matrix  [k]. 


[  k  ] 


[B]’^[D]  [B]  t  dA 

A 


C3]  [D]  [3] 


:  5 ) 


For  an  orthotropic  lamina  whose  material  axes  are 
oriented  at  an  angle  S  to  the  X  and  Y  directions,  the 
material  property  matrix  CD]  becomes  CQ],  of  Section  II  3: 
i.e. ,  the  stresses  are  related  to  the  strains  by 


r  1 

r  "1 

X 

X 

a 

£ 

y 

■  =  [Q]  - 

y 

T 

r 

1  xyj 

L  xyJ 

Since  [Q]  is  constant  for  each  increment. 


that 


[k:  =  [Bl'^CQ]  [31  t  A 


it  follows 


117) 


Now  to  account  for  n  number  of  plies  through  the 

» 

thickness  of  the  material,  an  equivalent  element  stiffness 
matrix  is  formed  as 
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where  i  denotes  the  i  layer,  and  here  n  =  16  plies.  Thus, 

loads  are  related  to  displacement  on  an  element  basis  by  the 
finite  element  equation; 


( f )  =  [ k  :  ( d }  (19) 

•q 

Note  that  even  though  multiple  layers  have  been  allowed 
through  the  thickness  of  the  element,  the  layers  have  common 
nodes.  Thus  strain  or  displacement  gradients  are  not 
permitted  through  the  thickness  of  the  laminate  [Cron',. 

Thus.  16  elements  stacked  together  all  deform  in  the  same 
way  under  a  given  load.  Whether  or  not  an  element  fails 
depends  on  its  own  material  properties  and  on  its 
orientation  to  the  given  load. 


? .  Non  1 1  near  F i n i te  Element  Theory 

As  discussed  in  Chapter  I,  a  more  accurate  analysis  of 
composites  must  include  the  effects  of  nonlinearity  in  the 
material.  The  computer  analysis  portion  of  this  thesis  was 
accomplished  using  a  nonlinear  finite  element  technique, 
which  13  based  upon  the  concepts  of  the  previous  two 
sections.  Sandhu  gives  a  complete  development  of  this 
technique  in  References  [5j  and  [6]. 

This  pr o g r es s 1 ve - p 1 y - f a  1 1 ur e  finite  element  computer 
program  (called  PLSTREN)  was  written  by  Sandhu  to  predict 
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'!)  the  damage  initiation  and  accumulation  process 
in  composite  laminates  containing  stress  concentrations  and 
v2)  the  static  strength  of  composite  laminates 
containing  stress  concentrations. 

The  program  combines  a  nonlinear  point  stress  analysis 
procedure  15,  6],  with  a  damage  accumu lation  approach  1251. 
The  program  is  modular  in  structure  to  permit  the  use  or 
varir'us  finite  elements,  ply  failure  criteria,  and 
p 0 s t - p 1 y - : a  1 1 ur e  unloading  models.  The  present  version  or 
the  program  uses  constant  strain  elements,  a  "total  strain 
energy"  failure  criterion  15,  6],  and  two  different 
pos t - p 1 y - f a  1 1 ure  unloading  models  13].  This  section 
develops  an  incremental  stress-strain  law,  discusses  the 
nonlinear  material  basic  properties,  describes  a  biaxial 
stress  state  correction,  and  presents  the  finite  element 
procedure  used  in  PLSTREN, 

To  account  for  nonlinearity  of  material  behavior,  an 
incremental  form  of  the  constitutive  law  of  Eq . (3)  is 
needed  to  define  the  response  of  the  laminate  under  general 
states  of  stress.  In  formulating  this  incremental 
constitutive  relationship  it  is  assumed  that  15 1 

a.  The  increment  of  strain  depends, upon  the  strain 
state  and  the  increment  of  stress;  and 

b.  The  increment  of  strain  is  proportional  to  the 
f  stress. 
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-V': ‘h  zhese  assumptions,  and  with  the  assumption  that 

the  3 t r a  1 n - i 1 s p 1 acemen t  relation  (matrix  lBj)  remains 

li.near.  for  orthotropic  iaminas  under  generalized  plane 

stress,  the  incremental  constitutive  law  can  be  written  as 

d£  =  S  (a-  )  do-  i  1  .  j  =  1  .  2  6)  2 

1.  w  I-  j 

where  da  ,  do-  are  strai.n  and  stress  increments  and  S  is 


a  function  of  the  current  strains,  a  .  Assuming  that  the 
la.T.ina  remains  orthotropic  at  all  load  levels,  Eq. 


expanded  to  become  13] 


dr 

1 
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L  <sJ 
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. 

dc 

L  ‘sJ 
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which  13  the  incremental  version  of  Eq.  (6).  Since  the 
la.mina  is  assumed  to  remain  orthotropic  at  all  load  levels 


each  increment  ;9] 

S  =  S  =  S  =  S  =3 

1<5  Z<S  <51  <32 


1  ar  1 

t y  ,  in 

Eq .  (21' 

d  O'  . 

1 

dr 

1 

normal  stress  and  strain  incre.ments  in  the 

fiber  direction, 

do-  . 

2 

dr  = 

2 

normal  stress  and  strain  increments  in  the 

transverse  direction. 

d  - 
a 

ic  = 

1  2 

shear  stress  i  .-icr  erne  n  ts,. 

rr 

■r1 

dr  = 

1  2 

shear  strain  mere  me  n  t  s  , 

s 

- 

terms  of  the  compliance  .matri.x  representing 

average  values  during  an  increment  of 
stress;  defined  by  Eqs .  , 
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Eq-.;a*:?n  ■  cfin  be  :n  verged  ’-.o  yield  ihe  i  r.  cremea '.a . 


version  o : 


=  C  t  i-  d  £■  ,  '■  1  .  1  =  1  ,  1  ,  6  i 

>- j  j 


wn  ere 


'•  j 


y  h  e  inverse  or  S 


I  r.  t,  h  e  application  of  the  mere  me  r.  t  a  1  constitutive  law 
multidirectional  laminates.  Eq .  1 

written  in  m.atrix  notation  as 


tor  the  K 


wn  ere 


-  , 


Hrt’. 


■:  d  £  :• 


stiffness  ma t r i x  of  t n e  x  P 1 y 

-  •  t 

stress  in  ere  .me  n  t  in  t  n  e  x  pi  y 

relative  to  the  .ma  t  e  r  i  a  -  axes  1  and  1 

.  .  It,  . 

strain  increment  in  tne  x  p.y 
relative  to  the  .material  axes  1  and  1 


Extending  Eq.  (24)  to  the  general  case  of  multidirectional 
and  a  n  g  1  e  -  p  1  y  1  a  .m  i  n  a  t  e  s  ,  the  incremental  constitutive  law 
for  the  k  ply  becomes; 


{ d  O' } 


c  c  ■ 


where  the  terms  of  IQj  are  given  following  Eq  .  .1/0.. 

To  obtain  the  incre  me  ntai  stress-strain  rel.ations  for  a 
1  a .T.  1  n a  .  the  definitions  of  the  S  terms  from  section  II  3 
are  substituted  into  Eq .  vTi;.  yielding 


<1 ) 

12 


n 


■A’h  ere  = 


and  provided  that,  d<y  and  do  ^0  [5] 

1  2 


The  incremental  elastic  constants  in  [Q’  (or  [SI;  are 
obtained  from  the  basic  lamina  stress-strain  data  by  an 
Iterative  procedure  [3],  In  general,  the  finite  element 
equation,  Eq.  (19),  presented  in  the  last  section  is  also 
written  on  an  incremental  basis: 

d(f}  =  Ck(e)  ]  did)  (27) 

•q 

where  i  •(  f  (■  ,  did),  and  [k(£)  ]  are  the  incre.ments  of  load 

•q 

ano  displacement  and  the  current  stiffness  matri,x, 
respectively.  Note  that  the  stiffness  matrix  depends  on  the 
current  level  of  strain.  This  dependence  exists  through  the 
transformed  reduced  stiffnesses,  [Q]  [19].  From  the 


previous  section. 


n 

k  ;-A-y[B]'^[Q][3]t 


and  one  can  see  that  [k  ]  is  calculated  using  the  i  [  Q  ] 

•q 

If  the  analyzed  ma  terial  exhibits  nonlinearity,  the  basic 
orooerties  of  the  laminas' 
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! o 1 1 o  W3  that  [ 4  ^  and  hence  [ k 


Ir 


■■n:  .  .  V  a n y  with  strain  as  we  11  [191. 

Since  material  properties  are  not  linear,  material 
property  curves  must  be  used  for  this  analysis.  For  this 
thesis  these  curves  were  obtained  from  stess-s train 
e  .\' p  e  r  1  men  1 3  on  unidirectional  lamina.  The  tests  and  their 
:  t  r  r  e  3  p  c  nd  1  .n  g  curves  and  material  properties  are  shown  in 
Table  1  .  These  t  e  s  t  s  are  described  in  T  h  a  d  t  e  r  17  of  this 


To  make  the  experimental  curves  usable  bv  the  computer 
code,  they  are  entered  in  tabular  form  and  are  reoresented 


analytically  by  a  piecewise  ciibic  spline  interpolation 


b 


:ur,o'.:on  .3^.  The  'jse  o:  'he  spline  I'lnc'.ion  renders  smooth 
3lress-s^^aln  curves,  which  are  desirable  lor  d  e  t,e  r  mi  r.  a  c  i  :■  n 
0 elastic  moduli  under  incremental  and  iterative 
computations  [3!;.  With  the  stress-strain  curves  thus 
represented  by  si.mpie  polynomials,  the  tangent  moduli  of 
S  ,  E  .  and  G  can  be  readily  obtained  as  functions  of 

Strain  by  differentiating  the  appropriate  cubic  spline 
functions  with  respect  to  strain. 

To  satisfy  the  numerical  nature  of  the  computer,  the 
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-n  the  program. 

In  this  procedure,  Ik  ]  is  initially  calculated  in  the 

first  incre.ment  using  engineering  linear  elastic  constants, 

which  are  equal  to  the  initial  slope  of  the  basic  property 

curves.  For  the  subsequent  increments,  Ik  I  is  calculated 

using  .material  properties  that  correspond  to  the  state  of 

strain  existing  at  the  end  of  the  previous  load  increment. 

Having  calculated  Ik  I  from  the  n*'*^  load  increment,  a  new 

•»q 

increment  of  load  is  applied  and  an  increment  of 
displacement  is  calculated  using  Eq  .  i  2'^ )  .  Thus  Eq  .  iT"' 


wi^l  take  the  form 


d{f}  =  [k(£)  ]  d{d}  (29) 

^•4*1  n  n-*'! 

where  n  denotes  the  n'''^  load  increment  [19  1. 

The  new  displacements,  did)  ,  can  then  be  used  to 

rn-l 

calculate  an  increment  of  strain  usinjg 

d{£}  =  [B]  did}  (30) 

The  incremental  stress  can  then  be  calculated  from 

d(o'}  =  d(£}  (3  1  ) 

where  i  denotes  the  i''*^  layer  as  defined  previously  :19], 
Given  the  current  levels  of  stress  and  strain,  new 
increments  of  stress  and  strain  are  computed  and  added  to 
the  end  of  the  current  level  to  obtain  a  new  level.  A  mean 
level  of  strain  is  then  computed  by  averaging  the  new  level 
of  strain  with  the  previous  level  at  the  end  of  the 
increment.  These  mean  strains  are  then  used  to  deter.mine  a 
new  set  of  material  properties  since  they  are  readily 
obtainable  as  functions  of  strain  through  the  cubic  spline 
functions  [19]. 

Given  the  new  elastic  properties,  [Q],  [Q],  and  hence 

[k  ]  are  recalculated.  The  same  load  increment,  not  an 

•  q 

additional  one,  is  then  reapplied;  and  the  increment  of 
displacement  and  strain  is  recalculated.  This  procedure  is 
continued  for  a  given  load  increment  until  the  change  of  the 
strain  increment  converges  to  less  than  a  specified  value, 
i . e . ,  until 
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0 . 130  1 


(32 


d£  -  d£ 

1  n 

d£ 

n 

where  n  denotes  the  n  reapp 1 icat ion  of  a  given  load 
increment.  When  Eq .  (32)  is  satisfied,  a  new  load  incremen 

IS  applied.  The  repetitive  use  of  the  procedure  outlined 
above  generates  the  predicted  stress-strain  response  of  the 
laminate  under  consideration  Cl9]. 

Graphically,  this  procedure  is  shown  in  Figure  12. 

IIG],  A  one -d 1  mens  1 ona 1  case  of  shear  loading  is  shown  for 
clarity,  but  the  program  recalculates  elastic  properties 
using  all  seven  curves  during  a  load  increment.  Stepping 
through  the  example,  the  initial  strain  exists  due  to  a 
previous  load  increment.  A  new  increment  of  the  load  is 
then  applied  and  the  resulting  strain,  is  calculated 

based  on  at  point  0.  A  new  modulus  is  now  determined, 

corresponding  to  the  average  of  and  point  A.  With 

new  moduli,  the  load  is  then  reapplied  and  a  new  strain,  y^ 
IS  calculated  using  the  modulus  at  point  A.  Again,  a  new 
modulus  correspond  1 ng  to  the  average  of  y^  and  y^.  point  B, 
IS  calculated  and  a  r eapp 1 i cat i on  of  the  load  yields  y^. 
This  procedure  is  repeated  until  [19] 

y  -  y 

1  n 

-  V  0.001  .33 

V 

i 

n 

Mote  that  superposition  principles  do  not  apply  in  the 
iterative  procedure  that  begins  with  Eq .  (29) .  The 
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s  t,  1  t  t  n  e  s  3  ma  ‘.rix  [kl  is  an  assemblage  of  element  stiffnesses 
bui*t  upon  the  (assumed)  independent  basic  material 
prcperties.  Thus  Figure  12  cannot  be  applied  directly  to 

Eq .  ; 29)  . 

During  the  determination  of  the  incremental  elastic 
constants,  direct  use  of  the  lamina  s tress -s trai n  curves  may 
be  iiicorrect.  For  the  off-axis  laminas,  uniaxial  loads 
create  a  biaxial  stress  state  in  the  lamina.  For  plies 

under  biaxial  stress,  the  biaxial  strains  £  and  c  must  be 

1  2 

modified  before  they  can  be  used  to  determine  the  elastic 
constants  from  the  cubic  spline  s tr es s - s tr a i n  curves.  It 

would  be  erroneous  to  use  these  strains  to  determine  E  or 

1 

from  s tress - s trai n  curves  obtained  under  simple  loading 
conditions.  Recalling  Eqs .  (26) ,  d£^  corresponds  to  the 

curve  ON  in  Figure  13  [5]  on  the  plane  OEHG.  The  simple 
s t r es s - s t r a  1 n  curve  OM  lies  on  the  plane  OEDC . 

3 tress - s trai n  data  similar  to  ON  is  not  avialable: 
therefore,  it  is  assumed  that  simple  equivalent  strain 
increments  can  be  computed  from  the  following  expressions 
[5  :  : 

Idc  d£ 

=  — i  =  - i _  (34) 

•q  E  ,  do' 

1  1  -  2 
t  2  - 

dC7 
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d<y 


•q 


d£ 


2  1 


do- 


dcr 


Note  that  the  error  that  would  result  from  usinj?  d^: 

I 

and  d£^  instead  of  their  equivalent  strain  would  be  •t>r 

d£  and  t-  / ,'?  for  d£  .  where  ri  =  dcr  /  da  . 

1  21  2  21 

Specifically,  the  nonlinear  finite  element  analysis  in 
this  thesis  uses  the  foregoinjg  methods  to  accomplish  the 
following  process.  Assuming  that  stresses  are  uniformly 
distributed  through  the  thickness  of  each  ply,  stress 
resultant  increments,  [dNj,  in  the  X,  Y  coordinate  system 
are  given  by 

p 

CdN]  =  y  Ida],  t,  36) 

K  < 

k  =  1 

wh  ere 


da  =  the  (transformed)  stress  increments  in  the  X.  Y 
coordinate  system 

Ih 

t^  =  thickness  of  the  k  ply 
p  =  number  of  plies  in  a  laminate 


Paralleling  the  development  of  Section  II,  B  and  substituting 
the  incremental  constituitive  law.  Eq .  (25) ,  into  Eq .  i 36 <  . 

the  stress  resultants  beco me 
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[dN  j  = 


y  :q:.  Id. 
<  < 


where  '.d£]  are  (transformed)  strain  increments  in  the  X.  Y 
coordinate  system. 

Noting  that  the  strain  increments  [d£'  are  the  same  for 

< 

all  plies  in  the  laminate,  Eq .  (38)  becomes 


iNj  =  lA:  [dsl 


wn  ere 


3  8  ) 


'  ^  > 


Inversion  of  Eq .  v39)  yields 

Cd^: 


:aj‘‘  IdN: 


i  40  i 


where  lA'.  represents  the  average  compliance  properties  of 

iK 

the  laminate  during  the  in-^l)  load  increment.  However, 

,  -1  th 

.Aj  is  not  known  when  the  n*l)  increment  is  applied. 

To  overcome  this  difficulty,  the  elastic  properties  at  the 
end  of  the  n‘^  load  increment  are  used  in  Eq.  .40),  i.e., 

:d£:  =  :a]'^  CdN:  t4i) 

The  strain  increments  id£ I  obtained  from  Ea.  (41'  are 
used  in  Eqs  .  (25)  ,  (34)  ,  and  (35)  to  calculate  Idc-l,  ,  Id^l  , 

K  '< 

These  stress  and  strai.n  increments  are 


I  ,  and  ds 

l  2 


Ih 


added  to  stresses  and  strains  at  the  n  load  increment  to 
obtain  the  current  stresses  and  strains  in  all  olies.  The 


/ 


curren'  stresses  and  strains  are  emp loved  to  determine  the 
average  elastic  properties  of  the  plies  and  a  new  lA',  ^  is 
computed.  This  procedure  is  repeated,  in  the  manner  of 
Figure  12,  until  the  difference  between  two  values  of 
Id^:'.  13  less  than  0.1  oercent.  Eq.  i4i)  then  beco me s 

[dri  =  [A]'*  IdNl  /42 

r."**4  •^♦1 

The  repetitive  use  of  the  procedure  outlined  above  generate^ 
the  3  t  r  es  3  -  3  t  r  a  1  n  response  of  unidirectional,  angle-pl.y,  :r 
mu* t 1 d 1 rect 1 onal  laminates  181. 

■3  .  Fai  lure  Cri  terion 

The  incremental  loading  procedure  outlined  in  the 
preceding  section  cannot  continue  indefinitely.  The 
prediction  of  the  response  of  the  laminate  should  ou*m*nate 
with  the  incremental  failure  of  elements  and  the  eventua* 
failure  of  the  laminate.  A  level  will  be  reached  where  a 
-amina  or  la.minas  can  no  longer  sustain  additional  loads. 
Thus,  the  failure  state  of  a  la.mi.na  under  genera*  stress 
states  13  determined  by  a  criterion  r  e  1  a  1 1 S  the  behavior  ;  : 
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n : n 1 i n  e  a  r  ma  terial  behavior. 


rollowin§  the  development  found  in  References  '5!  and 
.5'.  a  scalar  function,  f.  defining  the  failure  condition  of 
ma  terials  exhibiting  nonlinear  behavior  can  be  wr i 1 1  e  n  as 

fCo-.^r.Kj^l  (4  3' 

where  c  and  S  are  the  stress  and  strain  states  and  K 
describes  the  material  char ac t e r i s 1 1 cs  . 

A  simple  scalar  function  dependent  upon  both  stress  and 
strain  states  is  the  strain  energy  of  the  ma  t  e  r i a  *  .  In  the 
case  of  or  the  tropic  ma  teria.s.  strain  energies  that  are  due 
to  u,ni  axial  tension  and  compression  along  and  transverse  tc 
the  material  axes,  and  due  to  shear  are  independent 
parameters.  To  measure  the  level  of  the  effect  of  both 
stress  and  strain  states  on  the  orthotropic  materials  under 
a  combined  stress  state,  a  si.mple  scalar  function  is  a 
..near  combination  of  functions  of  different  strain 
energ.es.  Thus,  an  exp.icit  form  of  Eq.  43'  can  be  written 

d  3 


€ 


wr.ere  c  are  the  current  strain  components  and  m  are  the 


par ime 

r  s 
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U3ir.  g  n  cnreceat  1  ng  summation  notation 


- 

* 

O' 

6) 


t, 

1 


:  4  5 


’sin 4  results  of  tests  under  simple  load  renditions 

r 


a  dc 


=  1.2,5; 


vwr.  ere  f  are  the  u  1  t  i  ma  t  e  normal  tensi.e  or  oompress-ve 
and  shear  strains.  lombinir.ii?  Eqs  .  46  and  4"'  the 


;  ; .  . 0  w 1 n i  equation  is  obtained 


■  c 


4  8 


’.sirs  to  arcount  for  the  nonlinear  strains  that  occur  pri: 
;  0  fail  ur  e  119'.  in  compos  i  tes  such  as  Gr  -  PEEK  . 

Specialized  for  the  plane  stress  condition  and  usini? 
rontracted  notation,  Eq.  >44)  beco  me  s  15  1 
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*  ' 

1 1  r  0  c  t 

ion. 

‘.  ’Jn  1  0 ad  1  n g  s  c h e me s 

for 

the  comput 

er 

p  r  0  ^  r 

am  are 

d  e  s  c  r 

ibed  subsequently.) 

I  f 

the  fiber 

failure 

mode 

occurs  , 

the 

lamina  is  assumed  to 

be 

comp  1 e  te 1 y 

failed 

and  'he  loads  are  transferred  to  the  ad;acent  lamina.  The 
ocrurrenoe  of  either  lamina  failure  mode  is  determined  by 
us  in  2  the  strain  energy  failure  criterion  '9j. 

Fiber  failure  is  assumed  to  occur  if 


{  I,  'x  } 

L 

1  t.  {  u  } 


1:  0.1 
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and  the  criterion  of  Eq .  k48)  is  satisfied.  Matrix  failure 

occurs  when  Eq .  .43'  is  satisfied  and  inequality  i49;  is 

not.  It  13  further  assumed  that  the  delamination  moment 
coefficient  [22^  is  below  the  critical  value  for 
delamination  to  occur. 

A  final  point  to  present  is  the  modeling  of  the 
pr o g r e 3 3 1 ve - p i y - f a  1 1 ur e  [8].  Once  lamina  failure  occurs  in 
one  or  more  elements,  the  failed  ele me nts  have  to  be 
unloaded.  The  following  two  failed-ply  unloading  schemes 
ire  possible  in  the  program  PLSTREN . 

In  the  case  of  multidirectional  laminates,  it  is 
.•■el  0  h  9 1  the  failed  elementis)  unload  gradually.  This 

-''^i.s-n  able  since  as  areas  of  laminas  fail,  ad;  a  cent 


areas  and  iaminas  pick  up  the  extra  stresses.  Thus,  the 
a.nalvsis  attempts  to  model  interlaminar  stresses.  But 
according  to  Petit  [25], this  action  of  the  laminate  is  "not 
so  easily  postulated."  It  is  not  known  what  proportion  of 
the  load  the  lamina  retains  nor  what  happens  to  the  modulus 
after  degradation.  The  of  unloading  is  done  by  assigning 
negative  values  to  the  affected  moduli,  which  determined  on 
the  basis  of  the  failure  mode.  This  scheme  is  shown  for 
transverse  tension  in  Figure  15.  Note  the  slope  of  the 
negative  tangent  modulus.  This  slope  vanes  for  all  simple 
loads,  and  its  value  in  this  analysis  is  assumed  to  be  the 
negative  of  the  initial  slope  of  the  respective  basic 
property  curves. 

The  incremental  loading  is  continued  on  the  modeled 
laminate  until  the  stresses  in  the  affected  elementts)  are 
reduced  to  zero.  The  moduli  are  then  set  to  nominally  small 
values.  This  failed-ply  unloading  scheme  was  used  in 
Reference  [23]  and  was  found  to  give  excellent  results  for 
3 t r es s - s t ra 1 n  curves  of  multidirectional  laminates.  The 
scheme  attempts  to  model  the  progressive  nature  of  the 
matrix  cracking  process  within  the  elements.  But  while 
reasonable  for  multidirectional  laminates,  this  failed-ply 
unloading  scheme  gives  results  that  are  too  high  for 
unidirectional  laminates. 

In  the  case  of  unidirectional  laminates,  it  should  be 
assumed  that  the  failed  elements  unload  rapidly.  This  idea 
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IS  reasonable  since  failure  in  a  unidirectional  laminate 
occurs  throuigh  the  thickness  and  is  usually  sudden; 
theoretically,  no  interlaminar  stresses  exist.  The  affected 
moduli  of  a  failed  element  are  set  to  nominally  small  values 
(100  psi)  and  the  same  load  increment  is  applied  again. 

This  r eappl 1  cat  1  on  of  the  load  increment  can  cause  failure 
of  additional  elements.  These  failures  are  secondary  in 
that  they  are  caused  by  load  redistribution  without  an 
increase  in  load.  If  secondary  failures  occur,  the  load 
increment  is  repeated  until  no  further  secondary  failures 
occur.  The  load  is  then  increased  to  the  next  higher  level. 
Generally,  this  unloading  scheme  gives  lower  strength 
predictions  than  the  gradual  unloading  scheme  [8j. 

As  a  summary  of  the  development  presented  for  the 
program  PLSTREN,  a  diagram  of  Sandhu’s  program  is  shown  in 
Fi.gure  16.  The  chart  shown  is  not  a  typical  flowchart,  but 
instead  shows  the  complexity  of  a  nonlinear  progressive- 
ply-failure  finite  element  program. 

This  chapter  has  presented  the  theory  behind  the 
analytical  and  experimental  investigations  of  this  thesis. 
The  following  two  chapters  describe  procedures  used  in  the 
analysis  and  experimentation  phases.  Rer.ults  and 
conclusions  are  presented  in  subsequent  chapters. 
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As  discussed  in  Chapter  I,  the  analysis  of  composite 
laminates  to  predict  the  onset,  progression,  and  completion 
o:  failure  has  been  accomplished  by  many  methods.  These 
range  from  numerous  fracture  mechanics  appro-aches  to  severe 
finite  element  techniques.  In  this  thesis,  the  nonlinear 
incremental  finite  element  technique  described  in  Chapter  1 
was  used  on  three  stacking  sequences  of  Gt'PEEK  laminates. 
The  purpose  of  this  effort  was  to  study  the  initiation  and 
progression  of  failure  at  or  near  a  circular  discontinuity 
in  a  material  which  is  subjected  to  tension.  With  a 
particular  specimen  geometry  chosen,  a  suitable  finite 
element  mesh  developed,  and  finally  with  experimentally 
derived  data,  a  nonlinear  analysis  was  accomplished  for  the 
composite  .material  Gr/PEEK. 

-A  .  Specimen  Geometry 

To  analyze  the  failure  characteristics  of  Gr/PEEK,  th 
specimen  geometry  shown  in  Figure  17  was  chosen  for  this 
study.  Note  that  the  diameter  of  the  hole  is  one-third  the 
width  of  the  specimen.  As  mentioned  previously,  this  6'  X 
1.2'  X  0.084'  specimen  represents  a  thin  composite  plate  of 
finite  width.  These  tensile  coupons  were  manufactured  with 
Id  layers  of  laminas  with  the  prescribed  orientations,  or 
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Table  2. 

Modeled  Stacking 

Sequences . 

Stacking 

Sequence 

Symbo 1 

Type 

:  0 

0T 

Unidirectional 

1  <5 

;  90  : 

90T 

Unidirectional 

1  <S 

:  r4  5 : 

SH 

Angle-Ply 

stacking  sequences.  The  stacking  sequences  which  were 
modeled  and  analyzed  are  listed  above  in  Table  2.  Note  that 
"S'  implies  a  symmetric  laminate,  and  that  all  laminates  are 
symmetric  about  their  middle  surfaces. 

One  reason  for  choosing  these  lay-ups  was  to 
investigate  the  reaction  of  unidirectional  and  angle-ply 
laminates  under  tension.  It  is  hoped  that  these  results  can 
be  related  to  a  future  analysis  of  quas i - i so tr op i c  laminates 
under  tension.  Of  course  interlaminar  stresses  would  play  a 
major  role  in  the  reaction  of  the  plies  in  a  quas i - i so trop i c 
laminate.  It  is  hoped  that  by  studying  individual 
or  1 en tat  1 ons .  one  could  draw  conclusions  about  the 
quas 1  -  isotropic  specimens. 

3 .  F in  1 te  Element  Modeling 

An  acceptable  finite  element  model  of  a  structure 
should  accurately  predict  the  displacement  and  stress  fields 


7  2 
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::  'he  s  h  .-”jc  ur  e  ’_;nder  a  ^iven  load.  To  oboa-.n  oorh 

a  racde*.  ^ he  analyst  must  accurately  model  the  sice,  shape, 
loads  or  displacements; ,  and  boundary  conditions  o:  the 
actual  structure.  Furthermore,  he  must  provide  for  adequat 
refinement  of  the  element  mesh  in  areas  where  gradients  in 
the  displacement  field  or  the  stress  field  are  known  or 
expected  to  be  high.  In  describing  the  fi.nite  ele.ment 
mo  deling  carried  out  in  this  thesis,  each  of  these 
requirements  will  be  discussed  in  this  section. 

The  sice  and  shape  of  the  finite  element  mesh  obvious* 
.modeled  the  rectangular  thin  plate  analogy  that  was  shown  i 
Figure  1"^.  Each  of  the  sixteen  layers  of  laminas  were 
specified  in  thickness,  t^  of  Eq .  (36) .  and  in  orientations 

in  the  transformed  reduced  stiffnesses  of  Eq .  (10). 

Because  of  the  specimen  and  laminate  symmetry,  on*.y  a 
quarter  of  the  specimen  could  have  been  modeled  to  minimise 
.analysis  ti.me  and  computer  costs.  However,  in  the  program, 
the  stresses  and  strains  in  an  element  are  output  for  the 
center  of  the  ele me nt.  Therefore,  to  acquire  the  stresses 
and  strains  at  the  transverse  centerline  of  the  specimen, 
half  of  the  specimen  was  modeled  as  shown  in  Figure  13.  In 
this  way,  output  at  the  centerline  is  easily  rendered 
without  the  need  for  extrapolation. 

The  model  of  Figure  IS  was  used  for  all  three  types  of 
specimens.  Mote  the  coordinate  axes  relative  to  the  model, 
the  diagrammed  constraints,  and  the  listed  information  in 
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*  r.  -3  ;  .  ^  r  -3  .  A  1  3  r  r.  o  *  e  r,  h  a  t  h  e  e  1  e  me  n  Z3  are  :  .  r.  e  r 

3  r  : -jr.  a  :  r.  ^  *ne  hole.  Beca’^se  oi  higher  stress 
3  :  r.  c  e  r.  r  a  :  0  n  3  nearer  the  hole,  finer  d  i  s  c  r  e  1 1  z  a  t :  o  n  t :  the 
zpecimen  will,  yieli  more  accurate  results.  Mote  further 
that  no  transition  elements  Isuch  as  triangular  elements  :r 
elements  with  five  or  more  nodes)  are  used  in  this  mesh. 

A  0  0  7  r i 1 n  g  to  S  a  n  d  h  u  in  a  conversation  with  him;,  the 
crtiram  works  best  if  the  elements  are  as  square  as 
possible.  This  has  been  accomplished  with  the  quadri*at- 
era.s  gradually  increasing  in  size  away  from  the  hole  unti. 
thev  match  the  rectangular  grid  of  the  rest  of  the  speci.men. 

To  3  1  mu late  the  tension  field  on  the  specimen  in  the 
test  .machine  (See  Chapter  IV),  the  model  was  subjected  to 
tension  by  constraining  the  left  end  of  the  model  and 
displacing  the  other  end.  Another  method  of  this  program 
applies  'negative  pressure'  on  the  nodes  at  each  end  of  the 
specimen.  This  pressure  is  a  force  per  unit  length  acting 
on  each  end  and  creating  tension  in  the  specimen.  This 
method  was  not  used  because  one  of  the  nodes  on  which 
pressure  would  have  been  applied  was  fixed  in  this  model. 
Thus,  pressure  on  that  node  would  have  been  meaningless. 
Since  Its  counterpart  on  the  other  end  of  the  model  was  not 
fixed,  an  asymmetrical  stress  field  would  have  been  set  up 
in  the  specimen,  and  th4  mo  del  wo  u  1  d  not  have  s  i  ,mu  1  a  t  e  d  the 
actual  structure. 
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A  3  3  h  :  wn  i  r.  F.^ure  15,  '-he  modeled  ooundarv  r  :  r.  i 3 
were  3uor.  *,  har  ‘.he  node  a'-  rc?rdina'-e  ^.0,  3.^5  was  ed 

1  r.  both  '-he  x  and  y  directions  to  prevent  r  1 4 1  d  body  motion 
: the  .model.  The  longitudinal  centerline  of  the  spec:. men 
d - r  t  a  t  e  d  that  the  rest  of  the  nodes  along  the  line  y  =  3.3 
were  free  '.  o  move  in  the  x-direction  but  were  fixed  in  the 


7  -  d  1  r  e  c  1 1  1  n  . 

■At 

the  ends. 

the  nodes  a  1 

eng  the 

line  X  =  3.3 

we  re  :  -  .x  e  d  .  n 

t  h  e 

x-direct 

ion  but  we  r  e 

free  ',  o 

move  in  the 

V- di rect 1  in . 

The 

nodes  al 

0  n  g  the  line 

X  =  5.3 

we  r  e 

prescribed  with  a  give.n  displacement  and  were  free  in  the 
-  d  .  r  ec  1 1  c  n  .  The  rest  of  the  nodes  in  the  model  were  free 

1  .n  both  direct:  0, ns.  The  goal  of  this  model  was  to  simulate 
the  actuaf.  boundary  condition  caused  by  the  tension  grips  o 
the  test  machine. 

r  .  r  o  r. e  r  g  e  n  c  e  S  tudy 

Tt  determine  the  refinement  of  the  finite  element  mesh 
two  competing  goals  were  sought:  (1)  to  predict  behavior 

as  accurately  as  possible  and  (2)  to  save  C.PU  time  and 
expense  as  much  as  possible.  Thus,  a  convergence  study  was 
accomplished  using  the  three  models  of  various  refinement 
shown  in  Figure  19. 

Since  the  comparative  experimentation  had  not  been 
accomplished  at  this  time  in  the  study,  another  comparison 
was  needed  to  determ.ine  which  model  yie.l  ded  the  best 
results.  It  was  decided  that  modeling  the  isotropic 
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behavirr  o:  the  roup on  would  yield  a  good  basis  tor 
o 0 mpan 1 3  on  To  model  isooropic  behavior,  ohe  16  plies  were 
all  given  a  0-degree  orientation.  The  material  properties 
were  input  as  elastic  engineerin.g  constants  such  that 

E  =E  =18.0x10'* 

1  2 

E 

G  =  - - -  ,  V  =  0  .  3 

'  ^  2  (  1  *  io ) 

Using  these  isotropic  properties,  a  stress  concentration 
tactor  iSCF)  could  be  used  for  the  comparison  of  finite 
element  models. 

For  a  thin  isotropic  plate  of  finite  width  with  a  hole, 
the  3CF  was  determined  using  Figure  20  [27],  A  d/w  ratio  of 
1  3  yielded  a  K^,  or  SCF.  ol  3.46.  Thus,  the  stress  at  the 
hole  should  be  3.46  times  the  "far  field'  stress  (the  stress 
at  least  2  diameters  away  from  the  hole  and  near  the  center 
of  the  specimen)  for  an  isotropic  material  modeled  bv  the 
meshes  of  Figure  19. 

In  the  analysis,  each  mesh  was  given  a  1 :  s  p  .  i  :  =•  “>• 
loading  of  0.00  1  inch  and  stresses  were  oalrt.  .i’-'-. 
linear  version  of  the  program  describee  . 

Another  program  was  also  used  on  the  ~  t-  • 
for  greater  acouraev  with  a  .ecs--' 
pr'^gram  used  'reducei  .•-ri*. 
built  upon  the  z  r  z  ir  i~ 
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Table  3.  Results  of  Convergence  Study 


Model 

Stress  at 
y  =  0.2 
(ps  i  ) 

Far  Field 
Stress 
(ps  i  ) 

Stress 
Concen  . 

Accuracy 

% 

Time  of 
Execut ion 
( sec .  ) 

Coarse  1 

6909  . 

2890  . 

2 . 39 

69  .  1 

9  . 

Coarse  2 

6423  . 

2865  . 

2 . 23 

64 . 5 

200  . 

Med  1  urn 

8636  . 

2865  . 

3.01 

87 . 0 

34  . 

Fine 

9591  . 

2855  . 

3 . 36 

97  .  1 

46  . 

y  =  0.2  were  extrapolated  using  a  computer  algorithm  based 
on  Lagrange's  interpolating  polynomial  [29].  Note  that  the 
fine  model  provides  a  highly  accurate  result  with  little 
increase  in  CPU  time  over  the  others.  Also  note  that  the 
second  program  did  not  provide  increased  accuracy  and  that 
the  time  of  execution  was  significantly  higher  than  that  of 
the  others.  Thus,  the  fine  model  was  used  in  this  thesis  to 
carry  out  the  analysis  of  Gr/PEEK  under  tension. 

D .  Non  1 1  near  Anal  vs  is 

For  comparison  purposes,  this  study  could  have  first 
performed  a  linear  analysis  on  all  five  laminates  using 
constant  elastic  engineering  properties.  But  as  past 
studies  have  shown  [5,  6,  7,  8,  19]  and  as  the  experimental 

s tress -s tra 1 n  curves  of  this  study  show  (see  Chapter  V)  ,  a 
nonlinear  analysis  is  essential  for  an  accurate  study  of 
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composite  materials,  especially  Gr/PEEK.  Due  to  these 
reasons,  a  nonlinear  analysis  was  conducted  in  this  thesis. 

The  nonlinear  analysis  began  once  the  basic  properties 
were  determined  experimentally  (see  Chapters  IV  and  V) . 

These  properties  were  entered  into  each  of  the  data  files 
corresponding  to  each  type  of  specimen.  These  types  of 
input  cards  are  included  in  the  guide  to  the  program  in 
Append i x  B . 

Each  run  of  the  program  was  routed  from  a  CYBER 
mainframe  computer  to  a  CRAY  mainframe,  where  the  program 
was  executed.  Although  the  CRAY  performed  the  excutions  at 
a  higher  expense  per  second  than  the  CYBER,  the  time  of 
execution  was  significantly  faster,  thus  saving  overall 
expense  of  this  analysis.  Once  a  run  was  completed,  or  a 
user-imposed  time  limit  was  reached,  the  CRAY  routed  the 
completed  output  files  back  to  the  CYBER.  If  a  time  limit 
had  been  reached,  the  run  was  restarted  from  the  last 
increment  of  loading.  This  cycle  of  restarts  was  continued 
until  the  modeled  specimen  failed. 

For  each  increment  of  loading  during  a  run,  the 
stresses,  strains,  and  strain  energy  of  each  element  and  the 
displacements  of  each  node  were  sent  to  an  output  file. 

As  the  loading  increased,  elements  which  reached  their 
ultimate  strain  energy  failed  (see  Section  II  Q) .  The 
program  output  listed  this  sequence  of  element  failure. 
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For  each  run  of  the  laminates,  tension  stresses  were 
set  up  in  the  model  by  displacing  the  free  end  a  prescribed 
amount.  In  order  to  decide  what  increment  would  be 
appropriate,  for  simplicity  data  was  taken  from  the 
experimental  results  of  this  study.  The  displacement 
applied  was  simply  the  total  displacement  of  the  failed 
specimen  divided  by  ten  and  rounded  up  to  the  nearest  100 
of  an  inch.  Precision  was  not  critical  because  the  program 
adjusts  the  increment  to  suit  its  needs  in  analyzing  the 
failure  process. 

For  each  restart  of  the  program,  two  output  files  from 
the  previous  start  or  restart  were  important.  In  an  easily 
readable  format,  one  output  file  contained  the  displacements 
of  each  node  and  the  stresses,  strains,  and  strain  energy  of 
each  element.  It  also  plainly  stated  if  the  element  had 
failed.  This  output  file  was  routed  from  the  CYBER  to  a 
printer.  And  in  the  literally  thousands  of  pages  of  output, 
one  could  read  the  analytical  progression  of  failure  in  each 
of  the  specimen  lay-ups. 

The  second  output  file  contained  only  the  displacements 
of  each  node  for  each  increment  of  loading.  This  file  was 
routed  to  a  VAX  mainframe  where  it  becar.ie  an  input  file  to  a 
plotting  program.  The  VAX  was  used  because  it  operated  on  a 
baud  rate  that  was  8  times  faster  than  the  CYBER.  Thus,  for 
each  increment  of  loading,  a  plot  could  be  made  of  the 
distorted  model.  With  the  numerous  plots  of  the  distortions 
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caused  by  the  load  increments,  one  could  see  the  gradual 
analytical  deformation  of  each  specimen  under  tension. 

While  this  analysis  phase  was  taking  place, 
experimental  tests  on  the  specimens  were  being  conducted  as 
well.  Experimentation  is  discussed  in  the  next  chapter,  and 
the  results  of  both  analysis  and  experimentation  are 
presented  in  Chapter  V. 


I V .  Experimentation 

The  purpose  of  experimentation  in  this  thesis  was  to 
determine  the  ultimate  strength  and  failure  progression  in 
Gr/PEEK,  as  well  as  to  provide  basic  property  data  for  the 
computer  analysis  of  the  previous  chapter.  The  end  result 
was  then  a  comparison  of  ei^per imental  and  analytical 
results.  Experimentation  was  conducted  through  the 
Structures  Division  of  the  Air  Force  Flight  Dynamics 
Laboratory  (FDD  of  the  Air  Force  Wright  Aeronautical 
Laboratories  (AFWAL)  at  Wr i ght -Pat terson  AFB ,  Ohio.  Prior 
to  testing,  a  test  plan  was  submitted  to  the  Structures 
Division  from  the  Air  Force  Institute  of  Technology  (AFIT). 
The  test  plan  outlined  the  work  to  be  done  and  was  sent  to 
all  organizations  to  be  involved  in  this  project. 

Highlights  of  this  test  plan  are  presented  in  Appendix  A, 
along  with  an  explanation  of  the  logistics  involved  in 
accomplishing  experimentation  through  the  Structures 
Division  of  FDL . 

Specifically,  the  objectives  of  this  three-phase  test 
program  were: 

(1)  To  determine  basic  properties  for  the  thermoplasic 
composite  material  Gr/PEEK  at  room  temperature. 

(2)  To  determine  the  ultimate  tensile  strength  of 
Gr/PEEK  tensile  specimens  containing  a  0.4-inch  diameter 
circular  discontinuity  at  room  temperature. 
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13)  To  investi;?3te  the  progression  of  failure  in 
Gr/PEEK  at  a  circular  discontinuity  at  room  temperature. 

This  chapter  contains  the  description  of  fabricating, 
instrumenting,  and  testing  Gr/PEEK  coupons:  and  it  presents 
the  procedures  followed  in  accomplishing  the  above  test 
objectives.  Also  discussed  are  the  post-failure 
examinations  of  the  material,  which  were  videotapes  of  the 
ultimate  strength  tests  and  stereo  x-rays  of  the  failure 
progression . 


A .  Spec  1  men  Fabr i ca 1 1 on 

In  accordance  with  the  test  plan,  specimens  were 
fabricated  from  panels  of  Gr/PEEK  composite  panels.  These 
panels  were  supplied  by  the  Fiberite  Corporation,  a 
subsidiary  of  Imperial  Chemical  Industries  of  Great  Britain. 
The  panels  were  supplied  in  the  configurations  listed  in 
Table  4.  Note  that  all  lay-ups  contain  16  plies  and  are 
symmetric  about  their  midsurface.  Out  of  these  panels,  114 
specimens  were  fabricated,  instrumented,  and  tested.  Some 
surplus  Gr/PEEK  was  used  for  destructive  material  analysis 
to  ensure  that  material  specifications  were  met.  The 
remaining  material  was  reserved  for  future  testing. 

Each  panel  was  guaranteed  by  the  manufacturer  to  be 
APC-2  with  Hercules  AS4  graphite  fiber  in  a  Victrex  PEEK 
matrix.  The  fiber  content  was  61%  by  volume  and  687.  by 
weight.  A  random  sampling  of  the  panels  showed  that  the 
panels  supplied  met  the  material  specifications  of  APC-2. 
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Table  4 . 

Configurations  of  Panels 

Supp lied 

for  Specimens 

Laminate 

Stacking 

S 1 

ze 

Type 

Sequence 

(in. 

X 

1  n  . 

)  Quantity 

Uni-  [ 0  ] 

directional 

10 

X 

14 

5 

Un  1  - 

directiona 

<S) 

O) 

10 

X 

14 

5 

Angle-Ply 

[145] 

1 6 

X 

16 

3 

Quas  1  - 
isotropic 

[0/+45/90/-45] 

2S 

1 6 

X 

16 

ijuas  1  - 
isotropic 

l0/-*-45/'90  ] 

2  S 

i  6 

X 

i  6 

A  description  of  the  manufacturing  of  APC-2  panels  is 
contained  in  Appendix  B. 

When  the  panels  first  arrived  from  the  manufacturer, 
they  were  subjected  to  a  C-scan  to  determine  flaws  that 
might  have  disqualified  the  panel  for  fabrication  of 
specimens.  The  Non-Destructive  Evaluation  Branch  of  the  Air 
Force  Materials  Laboratory  conducted  the  C-scan  and  found 
all  panels  to  be  free  of  significant  flaws. 

After  the  C-scan,  specimens  were  fabricated  for 
testing.  Panels  were  cut  into  sub-panels;  then  sub-panels 
were  ‘tabbed. '  By  applying  tabs  to  each  end  of  the 
specimen,  the  testing  machine  was  given  an  area  to  grip  so 
that  tensile  and  compressive  forces  were  distributed  evenly 
without  damaging  the  specimen.  The  Beta  Corporation,  an 
in-house  contractor  of  the  Structures  Division,  supplied  the 
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tabbing  material,  which  consisted  of  i/'i6-inch  thick  G-10 
glass  epoxy  (0/90  woven) .  The  tabbed  subpanels  were 
cut  into  specimens  according  to  the  test  plan.  Also,  the 
specimens  for  the  ultimate  strength  and  failure  progression 
studies  had  0.3970-inch  holes  drilled  in  their  centers. 

The  final  tension  specimen  geometries  are  shown  in 
Figure  22.  These  are  ASTM  standard  tension  specimens.  Note 
that  the  widths  of  the  basic  property  specimens  vary 
according  to  lay-up.  These  widths  were  imposed  so  that  the 
ultimate  strength  of  the  specimen  did  not  exceed  the  load 
capacity  of  the  test  machine. 

To  experimentally  derive  compressive  properties,  there 
are  many  compression  specimen  types  from  which  to  choose. 

The  one  chosen  for  this  study  was  the  Rolfes  compression 
specimen  [30],  and  its  specimen  geometry  is  shown  in  Figure 
23.  The  tabbing  procedure  for  the  Rolfes  compression 
specimen  is  outlined  in  Figure  24.  The  end  result  of  this 
tedious  process  is  a  specimen  covered  mostly  with  tabbing 
material.  It  is  through  this  large  area  of  tabbing  material 
that  the  compression  fixture  grips  the  specimen  and  through 
shear  applies  compression  to  the  coupon. 

3 .  Spec  1  men  Instrumentation 

As  fabrication  of  specimens  was  completed,  the  speci¬ 
mens  were  instrumented.  All  specimens  were  instrumented 
with  appropriate  sizes  of  strain  gage  rosettes.  At  all  gage 
locations,  rosettes  were  placed  back-to-back,  i.e.,  they 
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were  adhered  to  each  face  .  Figure  25  shows  the  strain  gage 
rosette  configuration  and  numbering  scheme  for  the  front 
face  of  a  specimen.  The  rosette  on  the  back  face  simply 
mirrored  the  front  gage. 

Locations  of  each  gage  were  a  function  of  the  type  of 
test  conducted  on  the  specimen.  For  basic  property  tests, 
the  strain  gage  was  placed  in  the  middle  of  the  specimen.  A 
gaged  compression  specimen  is  shown  in  Figure  26.  For  the 
specimens  containing  a  hole,  gages  were  placed  both  at  the 
hole  and  at  the  far  field  location  for  some  specimens,  as 
shown  in  Figure  27.  Other  specimens  had  far  field  gages 
only,  as  shown  in  Figure  28.  The  use  of  rosettes  yields  the 
desired  plane  stress  information:  strain  in  the 
longitudinal,  transverse,  and  shear  directions  of  the 
specimen  axis  system. 

C .  Spec  1  men  Testing 

In  all  tests,  a  20-kip  Instron  universal  test  machine 
was  used  to  apply  tension  or  compression  to  specimens  under 
an  ambient  environment.  The  speci.mens  were  loaded  at  a 
constant  crosshead  travel  rate  of  0.05  inches  per  minute. 

The  crosshead  is  the  component  on  the  te'st  machine  which 
provides  displacement  loading.  This  speed  of  displacement 
was  chosen  as  median  speed  between  a  faster  impact  type 
loading  and  a  slower  creep  type  loading.  For  the  duration 
of  ail  tests,  it  was  assumed  that  creep  was  negligible  (see 
■Section  II  A)  .  Failure  loads  were  measured  using  a 
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20-kip  load  cell,  which  consisted  of  permanent  strain  gajges 
which  reacted  to  the  displacement  of  the  crosshead. 
Videotapes  of  all  the  ultimate  strength  and  failure 
progression  tests  were  made  in  an  attempt  to  document  the 
failure  process  and  provide  for  a  post-failure  analysis. 

The  Instron  test  machine  is  shown  in  the  left  half  of 
Figure  29  in  a  tension  application  mode.  The  specimen  is 
shown  in  the  grips  of  the  tension  fixture,  and  the  crosshead 
13  the  horizontal  member  just  below  the  test  fixture.  The 
load  cell  is  in  the  top  of  the  machine,  above  the  test 
fixture.  Also  shown  in  this  figure  are  the  voltmeter,  the 
video  camera  (both  shown  in  greater  detail  in  subsequent 
figures) ,  and  the  data  collection  devices. 

The  voltmeter  was  used  with  the  ultimate  strength  tests 
to  provide  a  visual  measure  of  the  load  in  the  specimen 
during  tests.  It  measured  load  by  measuring  the  resistance 
across  the  load  cell  during  the  test.  The  amount  of 
resistance  per  pound  of  load  varied  with  the  range  of  load 
for  the  test.  In  other  words,  if  the  machine’s  load  range 
for  a  particular  test  was  5000  pounds,  the  voltmeter  read 
1  volt  per  1020  pounds.  For  a  20,000  pound  load  range,  the 
voltmeter  output  1  volt  per  4082  pounds.  The  voltmeter  and 
specimen  are  shown  in  Figure  30,  with  testing  in  progress. 

The  reason  for  the  voltmeter  was  for  use  with  the 
videotape  of  the  tests.  With  the  videotape  filming  both 
the  failure  process  and  the  voltage  reading  (the  load)  ,  one 
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Figure  29.  Testing  Area. 


Figure  30.  Specimen  and  Voltmeter  During  Test. 
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could  perform  a  repeated  post-failure  analysis  on  each  test. 
Not  only  did  visual  occurrences  reveal  themselves  after 
testing,  but  the  failure  in  the  specimen  could  be  heard  by 
the  'tinking'  of  individual  fibers  and  the  g  -ating  sound  and 
sudden  boom  of  the  failed  specimen.  The  video  camera  and 
its  subject  are  shown  in  Figure  31,  with  testing  in 
progress  . 

In  the  background  of  Figure  29.  the  data  collection 
devices  are  shown.  For  each  test,  strain  was  read  from  all 
three  elements  of  each  strain  gage  and  load  was  read  from 
the  load  cell.  For  the  strain,  signals  from  each  leg  of 
the  rosette  were  read  by  a  strain  gage  conditioner,  sent  to 
an  amplifier,  then  to  a  multiplexer  (which  combines  the 
signals  for  the  computer) ,  and  finally  to  a  VAX  mainframe 
computer.  Load  signals  were  sent  in  a  similar  path  and 
combined  with  the  strain  information.  The  end  result  of 
each  test  was  a  great  volume  of  data  that  could  be  output  in 
any  desired  form. 

The  data  collection  devices  were  able  to  record  test 
data  in  varying  rates.  For  all  specimens  except  the 
90-degree  tension  coupons,  the  sampling  rate  was  2  data 
samples  per  second.  The  data  sampling  rate  for  the 
90-degree  tension  test  was  increased  to  4  data  samples  per 
second  since  these  tests  lasted  a  relatively  short  time.  In 
this  way,  data  collected  was  more  reliable,  and  strain  at 
failure  was  more  precise. 
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The  Instron  test  machine  is  shown  in  Figure  32  in  its 
compression  application  mode.  The  specimen  is  in  the  Rolfes 
compression  fixture,  and  the  crosshead  is  applying  load  by 
moving  down.  Here,  the  load  cell  is  under  the  fixture, 
inside  the  rectangular  box.  Load  cell  data  is  fed  through 
the  cable  attached  to  the  box  and  over  to  the  data 
collection  devices.  In  this  figure  also  note  the  plotter  in 
the  upper  left.  Plotted  for  each  test  was  a  continuous  load 
versus  displacement  curve.  A  ball  point  pen  simply  traced 
the  load  as  the  scroll  moved  with  the  displace  of  the 
crosshead . 

The  Rolfes  compression  fixture  is  shown  in  the  test 
machine  in  Figure  33.  One  pair  of  wires  for  each  leg  of 
each  gage  is  extended  from  the  gages  to  the  terminal  strip 
at  the  left  of  the  figure.  From  this  terminal  strip, 
signals  are  sent  to  the  data  collection  devices. 

Further  views  of  the  compression  fixture  are  shown  in 
Figure  34.  Figure  34(a)  shows  the  unassembled  fixture  with 
one  end  of  the  specimen  in  the  grips.  Figure  34(b)  shows 
the  assembled  fixture.  Note  that  the  grip  in  the  front 
extends  beyond  the  cover  plates.  It  is  on  this  grip  that 
the  crosshead  applies  the  displacement  ipad. 

The  entire  testing  program  is  outlined  in  Table  5. 

Note  the  test  symbols  listed.  These  will  be  used  for 
conciseness  in  the  next  chapter.  Also  note  the  quantities 
of  specimens  required  for  each  phase  of  the  testing.  For 
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Figure  32.  Instron  Test  Machine  in  Compression  Mode. 


Figure  33.  Compression  Specimen  Wired  for  Test. 
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determining  basic  properties,  five  obstensibly  identical 
tests  of  each  type  of  specimen  were  assumed  to  provide  a 
good  average  of  material  property  data.  For  the  ultimate 
strength  tests,  four  specimens  with  gages  at  the  hole  and  at 
far  field  were  tested  along  with  three  specimens  having  only 
a  far  field  gage.  It  was  desired  to  provide  as  much  data  as 
possible  for  this  phase.  For  the  remaining  phase,  three 
tests  for  each  type  of  specimen  were  assumed  to  give  a 
statistically  accurate  account  of  results.  In  all  phases, 
extra  specimens  were  manufactured  in  the  event  of  poor 
results.  Occasionally  in  this  type  of  testing,  strain  gages 
fail  or  tab  adhesive  fails  prior  to  failure  of  the  specimen. 

D .  Bas i c  Property  Tests 

To  determine  the  basic  property  data  of  APC-2,  the 
basic  property  specimens  were  gaged  at  their  centers  and 
tested  under  static  loading  at  room  temperature.  The  goal 
of  these  tests  was  to  provide  sets  of  values  from  which  the 

T 

computer  program  could  derive  the  nonlinear  values  of  E  E 

,  E  ^ ,  E  ,  G  ,  V  ^ ,  and  The  tests  were  listed  in 

1  2  2  12  12  12 

Table  1  of  Chapter  II,  along  with  their  corresponding  curves 
and  the  basic  property  derived  from  the  data.  Note  that  the 
[t45]  laminate  was  used  for  determining  the  t  vs. 

4S  12 

,v  ^data.  According  to  ASTM  Standards  [31],  the  tensile 
coupon  of  ±45  laminas  yields  shear  stress-strain  results 
from  the  following  equations: 
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Table  5. 

Summary  of  Testing 

Program 

Basic  Property 

Tests  (Specimens 

with  no  hole) 

Spec  1  men 

Number  of  Number  of 

Type  o  f 

Type.'  Test 

Spec  1  mens 

Gages 

Gage  • 

0^  Tension 

6 

12 

1 

X 

0  To  mp  r  e  s  s 1 o  n 

5 

10 

* 

o 

90  Tension 

9 

18 

9  0  3  c  mp  r  e  s  s 1 0  n 

8 

16 

- 

■3 

r45  Tension 

6 

12 

n 

34  Specimens 

Tested 

U 1 1 : ma t e  Strength  Tests  (Specimens  with  a  0.4"  dia  hole) 

Type  of 
Gage  » 

2 
3 


3 


3 


7  c  n  t  1  n  u  e  d  . 


Specimen 

Test 

Number  of 

Number  of 

Type 

Symbol 

Spec  1  mens 

Gages 

; 

i  <3 

0T 

3 

3 

6 

6 

:90  : 

1  <5 

90T 

3 

6 

4 

8 

-  -45 :  ^ 

4S 

SH 

3 

6 

4 

3 

:  0  '  *45/90/ -45  ] 

2S 

Q1 

3 

6 

4 

8 

;  0  —  4  5  /  9  0  : 

2  S 

Q2 

3 

5 

4 

3 

Totals 

15 

30 

19 

3  3 

34  Specimens  Tested 
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Table  5  (Continued) .  Summary  of  Testing  Program 


Failure  Progression  Tests  (Specimens  with 

0 . 4  • 

d 1  a  hole) 

Specimen  Test  Number  of  Specimen  at 

Percentage 

Type 

Symbo 1 

»  1 

#  2 

»  3 

[0  ] 

1  a 

0T 

3 

3 

O 

[90  ] 

1  <s 

90T 

3 

3 

3 

[145  j 

4S 

SH 

3 

3 

J0/-I-45/90/ 

-45 j  Q1 

2S 

3 

3 

3 

[  0/  -^45/90  ] 

Q2 

2S 

3 

3 

3 

Totals 

15 

15 

1  4 

44  Specimens 

Tested 

* 

Types  of 

Strain  Gages 

Reference 

Number 

Number 

Type  of  Gage 

App 1 1  ed 

R  e  ma  r  k  s 

• 

CEA-03-062UR-350 

190 

■-I 

CEA-03- 125UR-350 

38 

Far 

field  gage. 

3 

WK-03-060WR-350 

40 

stacked  rosette 

at  hole. 

>  50 ; 


where 


T 

12 


P  /  2  b  d 


i  5  1  ) 


=  load  in  the  specimen  longitudinal  axis 
b  =  specimen  width 
d  =  specimen  thickness 

,  y  ,  €  ,  and  c  have  been  defined  previously. 

12  12  X  y  ^ 

o 

Note  that  Eqs .  (50)  and  (51)  are  derived  for  the  t45 

laminate  through  the  use  of  a  coordinate  transformation. 
Taking  Eq .  (51) ,  for  example,  to  get  the  desired  shear 

strain  in  the  material  axis  system,  one  would  use  this 
transformation  equation; 


cos^d  sin^d  2  sin  &  cos  B 

2  2 

sin  ^  cos  i9  -2  sin  ^  cos  9 

2  Z 

-sin  9  cos  9  sin  9  cos  9  cos  ^-sin  9 

With  6  =  +45  or  -45  ,  Eq .  (51)  is  realized. 

To  obtain  the  desired  stress - s trai n  data  sets  which 
were  entered  into  the  computer  program  of  this  study,  t,he 
stress  and  strain  values  derived  from  each  test  were  plot;ted 
on  a  single  graph  and  an  average  s  t  r  ess -S'tr  a  i  n  curve  was 
determined.  For  example,  to  determine  <y ^  vs.  six  curves 

were  plotted  of  the  axial  stress  (c/^  =  o'  ^  =  P  /  A)  versus 
the  maximum  strain  =  e  )  read  from  the  strain  rosette 

1  X 

for  the  0-degree  tension  specimens.  Ail  six  experimental 
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cr  MS.  £  curves  were  outout  on  one  plot,  and  a  curve 
that  was  based  on  least  squares  fit  was  plotted  through  the 
SIX  experimental  curves.  Values  of  a ^  vs.  were  then 

taken  from  this  curve  and  tabulated.  These  basic  property 
curves  are  presented  in  Chapter  V  with  other  test  results. 

Other  values  were  tabulated  as  well.  To  determine  the 
Poisson  ratios,  both  the  tension  and  compression  curves  of 
vs.  and  o ^  vs.  were  plotted  on  one  graph  and  values 

of  and  £^  were  chosen.  Dividing  £^  into  £^  yielded  the 
desired  results.  Additionally,  the  engineering  elastic 
constants  were  needed  to  give  the  computer  program  values  of 
basic  properties  from  which  to  start.  The  constants  were 
the  initial  values  of  basic  properties.  The  engineering 
elastic  constants  derived  from  these  tests  were  compared 
with  others’  findings  and  found  to  be  accurate.  The 
tabulated  values  of  nonlinear  material  properties,  along 
with  the  elastic  engineering  constants  were  entered  into 
data  files  for  the  analysis  portion  of  this  thesis. 

E .  Ultimate  Strength  Tests 

To  determine  laminate  s tress - s tr a i n  responses  to 
ultimate  tensile  load,  the  0T ,  90T .  SH,  Ql,  and  Q2  specimens 
were  tested  under  tension.  All  seven  specimens  of  each  of 
the  five  types  of  lay-ups  were  instrumented  with  a 
( CE A- 03 -  1 25UR- 3 50 )  strain  rosette  away  from  the  hole,  the 
far  field  gage,  to  determine  strain  to  failure  for  the 
laminate.  Four  of  the  seven  specimens  were  instrumented 
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with  a  'stacked'  rosette  ( WK -03 -062UR-350 )  at  the  side  of 
the  hole.  tSee  Figure  28  in  Section  B  of  this  chapter.) 

The  use  of  a  stacked  rosette  allows  one  to  determine  the 
strain  in  a  much  smaller  area.  Because  of  the  extremely 
large  stress  gradients  at  the  hole  (See  Figure  21  in  Section 
III  C) ,  use  of  a  rosette  with  three  grids  spread  apart,  as 
in  Figure  25.  would  produce  meaningless  results. 

Refering  to  the  finite  element  model  of  the  specimen 
geometry,  the  small  stacked  rosette  was  centered  on  the 
second  element  from  the  hole  on  the  transverse  centerline  of 
the  specimen.  The  rosette  approximately  covered  the  area  of 
that  element.  The  far  field  gage  was  placed  on  a  group  of 
elements  where  stresses  were  assumed  to  be  uniform. 

Figure  35  shows  this  strain  gage  placement  with  respect  to 
the  finite  element  mesh. 

For  each  of  the  five  types  of  laminates,  speci.mens  were 
subjected  to  tension  until  they  failed.  The  loads  and 
strains  throughout  the  test  and  at  failure  were  recorded. 
Results  from  the  small  rosette  at  the  hole,  as  well  as  from 
the  gage  at  far  field,  were  used  for  comparison  to  the 
computer  analysis.  The  tensile  stresses  at  failure  were 
averaged  to  obtain  average  failure  stresses.  These  average 
stresses  were  not  only  compared  to  analytical  findings,  but 
they  were  also  used  to  accomplish  the  final  phase  of 
testing . 
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A  videotape  was  made  of  all  the  tests  to  ultimate 


strength.  From  repeated  viewing,  one  could  see  and  hear 
occurrences  that  were  not  noticed  during  the  test.  Also, 
stills  (photographs)  were  taken  of  various  frames  which 
showed  the  initiation  and  progression  of  failure,  as  well  as 
deformations  and  the  ultimate  breaking  of  the  specimen. 

F .  Failure  Progress  ion  Tests 

To  determine  the  initiation  and  growth  of  ply  failure 
in  laminates  with  a  0.4-inch  diameter  hole,  0T ,  90T ,  SK,  Ql, 
and  Q2  specimens  were  tested  under  tension.  For  three 
percentages  of  the  average  failure  stresses,  three  specimens 
of  each  lay-up  were  tested.  The  specimens  were  placed  in 
the  grips  of  the  Instron  test  machine  and  loaded  to  a 
predetermined  load.  This  load  was  derived  by  multiplying 
the  predetermined  stress  percentage  by  the  cross-sectional 
area  of  each  specimen. 

For  each  average  failure  stress  determined  above,  three 
percentages  of  failure  load  were  chosen.  For  each  of  the 
five  laminates,  the  results  of  the  ultimate  strength  tests 
were  studied.  The  stress-strain  curves  were  examined,  and 
the  videotape  was  viewed  repeatedly.  From  this  post-failure 
analysis,  the  percentages  of  failure  stress  were  chosen  and 
are  listed  in  Table  6.  These  percentages  were  based  on 
occurrences  during  the  failure  process  that  warranted 
investigation.  It  was  hoped  that  testing  to  these 


Table  5.  Percentages  of  Average  Failure  Stress 


(to  which  each  laminate  was  tested) 


Specimen 

Percentage  *  1 

#  2 

#  3 

0T 

50 

80 

90 

90T 

80 

90 

95 

SH 

70 

85 

95 

Q1 

80 

90 

95 

Q2 

80 

90 

95 

percentages  would  render  visible  results  through  the 
technique  of  stereo  X-ray. 

G .  Enhanced  Stereo  X-Ray 

Once  laminates  had  been  loaded  to  a  percentage  of 
failure  stress,  they  were  taken  to  the  Air  Force  Materials- 
Laboratory  for  analysis.  There  the  growth  of  ply  failure 
was  photographed  using  the  techniques  of  stereo  x-ray.  To 
determine  which  plies  failed  in  a  laminate  under  tension 
around  the  hole,  each  set  of  three  specimens  was  subjected 
to  'opaque  penetrant  enhanced,  three-dimensional, 
x-radiaography ■  [32],  By  using  this  technique  one  can  see 

failure  through  the  entire  laminate  thickness  and,  ideally, 
pick  out  the  extent  to  which  each  ply  has' failed.  This 
visual  analysis  requires  some  education  and  acclamation  for 
the  person  who  has  never  performed  stereo  viewing  of 
composites . 
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At  the  facilities  of  the  Non- Des true t i ve  Evaluation 
Branch  of  the  Materials  Lab.  all  44  specimens  were 
photographed  using  the  techniques  of  Reference  [32]. 

The  specimens  were  penetrated  with  an  x-ray  opaque  fluid  and 
then  photographed  onto  x-ray  film  in  two  images.  The  right 
image  was  shot  first  while  a  thin  plate  of  lead  was  placed 
over  the  left  side  of  the  film.  Then  the  left  side  of  the 
the  x-ray  was  shot  with  the  lead  plate  on  the  right  side. 
This  set-up  IS  shown  in  Figure  36.  Once  the  film  was 
developed  and  stereo  x-rays  were  completed,  the  analysis  of 
the  images  began.  This  analysis  involved  using  a  stereo 
viewer  with  a  magnifying  lens  and  picking  out  the  damage 
progress  ion. 

With  the  experimentation  described  in  this  chapter, 
along  with  the  methods  for  evaluating  failure  in  specimens, 
results  were  obtained.  These  results  are  presented  next. 

It  was  hoped  that  this  combination  of  experiment  and 
analysis  would  give  a  complete  picture  of  how  APC-2 
laminates  fail  under  tension  around  a  hole,  and  whether  or 
not  the  computer  analysis  used  can  accurately  predict  this 
behav i or . 
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This  chapter  presents  and  discusses  the  analytical  and 
experimental  results  of  this  thesis.  The  computer  program 
embodying  the  concepts  of  Chapter  II  and  applied  in  Chapter 
III  was  used  to  determine  the  theoretical  responses  of 
[  (3  ],  [90  ].  and  [i45]  Gr/PEEK  laminates  subjected  to 

1  <5  1  <S  4S 

in-plane  uniaxial  loads.  These  responses  were  then  compared 
to  the  experimental  results  of  Chapter  IV.  The  experimental 
study  of  the  two  quas i - isotropic  laminates  follows  the 
ana  1 y 1 1 ca 1  - exper 1  men tal  study  of  the  unidirectional  and 
angle-ply  laminates.  But  first,  the  basic  property  results 
are  presented  showing  best  fit  curves  and  listing  the  data 
that  was  input  to  the  computer  program. 

A.  Results  o f  3as i c  Property  Tests 

The  results  of  each  test  conducted  during  this  phase  of 
the  experimentation  are  listed  in  Tables  7  through  11.  For 
each  table,  remarks  are  included  for  completeness.  Note 
that  the  average  failure  stresses  shown  in  the  tables 
represents  the  strength  o f  Gr / PEEK  under  simple  loading 
cond 1 t 1 ons , 
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Table  7.  Basic  Property  Tests  of  0-deg  Tension  Specimens 


Failure 

Failure 

Test 

X-Sect 

Load 

Stress 

Mumber 

Area 

(lbs  .  ) 

(  p  S  1  ) 

1101 

0 . 04983 

14 . 900  . 

298 , 700 . 

1002 

0 . 04988 

15 .700  . 

314 .900 . 

1003 

0 . 0498 

14 . 380  . 

288 . 300 . 

1004 

0 . 05013 

15 . 400  . 

307 . 200  . 

1005 

0.05107 

13 . 850  . 

271 . 200 . 

1006 

0 . 0498 

15 . 400  . 

309 , 200 . 

Average 

Failure  Stress  = 

293 , 300 . 

Table  3 . 

Basic  Property 

Tests  of  0-deg 

Compression 

Specimens 

Failure 

Fai lure 

Test 

X-Sect 

Load 

Stress 

Numbe 

r  Area 

( lbs  .  ) 

( ps  1  ) 

Notes 

550  1 

0 . 06083 

A 

5502 

0 . 05917 

9 . 439 . 

159 . 500  . 

5503 

0 . 05933 

10.701. 

180 .350 . 

5504 

0 . 05992 

9 . 850  . 

164  ,  400  . 

3 

5505 

0 . 0594  1 

1  1 . 360 . 

191  .  200 . 

Average  Fai 

lure  Stress  = 

173 . 900 . 

Notes  : 

A .  A  tab  to 

one  of 

the 

strain 

gages 

was  lost. 

and  the 

specimen  was 

scapped 

8  Slippage 

of  the 

grip 

on  the 

tab 

occurred . 

Test  was 

redone  Fir 

St  load 

was 

5,604  i 

bs  . 

Second  and 

ul 1 1  mate  1 oad 

was  ;  1  .  1  6  1  1 

bs  . 

Table  9.  Basic  Property  Tests  of  90-deg  Tension  Specimens 


Test 


X-Sect 


Failure 

Load 


Failure 

Stress 


Number 

Area 

200  1 

0 . 07884 

2002 

0 . 08108 

2003 

0 . 08008 

2004 

0 . 08208 

2005 

0 . 08325 

2006 

0 . 07992 

2007 

0.08116 

2008 

0 . 083 

2009 

0 . 08  192 

Average  Failure 


libs  .  ) 

( ps  i  ) 

1 . 200 . 

15 , 220  . 

1,183. 

14 , 590 . 

1,134. 

14  .  160 . 

1 . 209 . 

1 4 , 730 . 

1 ,210. 

14 . 530 . 

920 .  - 

990  . 

12  ,  200  . 

560  .  - 

1  . 092  . 

13  .330  . 

Stress  = 

14,110. 

Table  10,  Basic  Property  Tests  of  90*deg  Compression  Specimens 


Failure 


Failure 


Test 


X-Sect  Load 


Stress 


Number 

Area 

( lbs  .  ) 

( ps  1  ) 

Notes 

460  1 

4502 

0 . 06208 

0 . 06225 

2.113. 

33 , 940 . 

A 

4503 

0 . 062  17 

2,101. 

33 . 790 . 

4504 

0 , 06233 

2,144. 

34  .  400  . 

4505 

0 . 06004 

2 . 082 . 

34 , 680 . 

4506 

0 . 06008 

1 . 999 . 

33 . 270 . 

4507 

0 . 060  16 

1 .987 . 

33 . 030 . 

4508 

0 . 05933 

2 , 036 . 

34 , 320 . 

Average  Failure  Stress  =  33,920. 


Notes  : 


A.  Data  obtained  i;  this  test  was  not  useable. 
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Table 

1  1  .  Bas 1 

c  Property 

Tests  of 

Shear  Specimens 

Failure 

Duration 

Fai lure 

Test 

X-Sect 

Load 

of  Test 

Stress 

Number 

Area 

(lbs  .  ) 

(min. ) 

( ps  i  ) 

Notes 

300  1 

0 . 083  17 

4000  . 

22 . 0 

A 

3002 

0.08159 

4200  . 

18 . 3 

B 

3003 

0 . 08200 

4400  . 

18.3 

B 

3004 

0 . 08300 

4670  . 

38 . 3 

56 , 260  . 

3005 

0 . 08  400 

4800 

40 . 2 

57,140. 

3006 

0 . 08  400 

4630  . 

38 . 6 

55,120. 

Average  Failure 

Stress 

=  56,170. 

Notes  : 

A.  This  specimen  was  loaded  twice.  When  it  was 
realized  the  duration  of  the  test  was  much  longer  than  the 
other  laminates,  this  first  specimen  was  unloaded  to  check 
for  slippage  of  the  grips  or  the  tabs.  It  was  then  reloaded 
at  a  faster  rate  (0.1  inches  per  min.). 

E.  These  specimens  were  loaded  at  the  rate  of  0.  1 
inches  per  minute.  However,  it  was  realized  that  this  load 
rate  was  too  fast,  and  the  load  rate  was  switched  back  to 
0.05  inches  per  minute  for  the  rest  of  the  specimens. 

For  the  shear  i t45 1  specimens,  these  results  were 
Significant: 

' 1)  All  of  the  coupons  reached  strains  that  were 
beyond  the  linear  assumptions  of  this  study.  The  average 


strain  was  estimated  to  be  at  least  207.,  based  on  measuring 
the  tina*  gage  lengths  of  the  failed  specimens. 

■2’  Scissor ing  of  the  fibers  occured  in  each  specimen. 
Scissoring  is  a  phenomenon  where  the  fibers  in  an  angle-ply 
laminate  try  to  align  themselves  with  the  load.  The 
graphite  fibers  in  the  ductile  PEEK  matrix  accomplished  much 
scissoring  before  ultimate  failure. 

3c*  h  c :  t  h  e  above  a  n  c  ma lies  can  be  seen  in  the  railed 
shear  spe  cim,ens  of  Figure  3'’. 

The  original  stress-strain  curves  of  these  3he.ar  tests 
shewed  that  a  failure  point  needed  to  be  defined  before  the 
actual  failure  of  the  specimen.  On  the  original  curve, 
after  strains  reached  approximately  5.5  7.,  the  slope  of  the 
curve  began  to  flatten  out  and  continue  to  about  a  207. 
strain.  It  was  believed  that  both  the  large  strains 
indicative  of  thermoplastics  and  the  scissoring  effect  of 
the  angle-ply  laminate  were  responsible  for  this  apparent 
p.ast.city.  As  an  initial  estimate  of  reasonable  failure, 
the  data  derived  from  the  shear  tests  was  cut  off  at  a  value 
near  5.5  percent  strain.  This  value  of  strain  is  slightly 
above  the  approximate  bound  on  linear  assumptions  in 
?.  eference  [33,.  Strains  occurring  beyond  5  or  6  percent 
approach  phenomena  described  by  the  principles  of  large 
i . sp 1  ace men t  theory,  where  partial  second  derivatives 
begin  t  c  enter  the  a  n  a  1 y  s i s . 
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Basic  Property  Shear  Specimens 
Before  and  After  Testing. 


The  s tress - s tra 1 n  curves  of  each  of  the  basic  property 
tests  were  plotted  from  experimental  data  as  described  in 
Chapter  IV,  From  each  set  of  stress-strain  curves  a  best 
fit  curve  was  determined.  These  best  fit  curves  are  shown 
in  Figures  38  through  44. 

From  these  seven  curves,  corresponding  stress  and 
strain  values  (or  Poisson  ratios  and  strain  values)  were 
picked  off  and  tabulated.  These  values  were  the  basic 
property  input  for  the  computer  program  as  described  in 
Chapter  III.  These  bas i c  pr oper  1 1  es  o  f  APC -  2  Gr / PEEK  are 
listed  in  Tab  1 e  12. 

As  described  in  Chapter  IV,  initial  values  for  the 
engineering  properties  must  be  input  into  the  computer 
program  for  it  to  begin  its  analysis.  For  Young’s  modulus 
and  the  shear  modulus,  these  initial  values  were  simply  the 
initial  slopes  of  the  basic  property  curves  ( s tress / s trai n 
for  the  first  data  point) .  For  the  Poisson  ratios,  the 
initial  values  were  extrapolated  from  the  data  above.  Th.ese 
initial  values  are  given  as  engineering  elastic  cons  tan  ts 
for  APC -2  Gr/P£EK  in  Table  13. 

B.  S t r es s - S t r a  1 n  Responses  of  [0  ].  [90  ],  and  [P45] 

Laminates 

As  described  i  .n  Chapter  IV,  0T ,  90T  ,  and  SH  laminates 

with  holes  were  tested  under  tension  to  their  failure  load. 
Also  tested  to  ultimate  strength  were  the  quas i - i s o tr op i c 
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Table 

12.  Basic  Material 

(Stress-strain 

Properties  of  APC-2  Gr/PEEK 

data  from  basic  property  tests) 

0 

Tension 

A 

0 

Compress  ion 

Strain 

Stress 

Po  1  sson 

S  train 

Stress 

Poisson 

( in/ in) 

( ks  1 ) 

Ratio 

( 1 n/ 1 n ) 

( ks  1  ) 

Ratio 

0.001 

19.5 

0.31 

0 . 0008 

15.0 

0 .35 

0 . 002 

39 . 0 

0.315 

0.0016 

30 . 0 

0 . 356 

0 . 003 

58 . 5 

0 . 3167 

0 . 0024 

44 . 5 

0 . 359 

0 . 004 

78  .  1 

0.3165 

0 . 0032 

58 . 5 

0.3615 

0 . 005 

97 . 8 

0.315 

0 . 004 

72 . 5 

0 . 3635 

0 , 006 

118.0 

0.312 

0 . 0048 

36 . 0 

0 . 365 

0 . 00'’ 

138.5 

0 . 3085 

0 . 0056 

99 . 0 

0 . 3662 

0 . 006 

159.0 

0 . 305 

0 . 0064 

111.5 

0 . 3673 

0 . 009 

180.0 

0 . 3025 

0 . 0072 

124.0 

0 . 3684 

0.01 

202 . 0 

0 . 3 

0 . 008 

136.0 

0 . 3695 

0.011 

224 . 5 

0 . 298 

0 . 0088 

148.0 

0 . 3705 

0.012 

247 . 5 

0 . 296 

0 . 0096 

159.0 

0.3715 

0.013 

271.0 

0 . 2943 

0.0104 

170.0 

0 . 3725 

0.014 
*  0.0141 

295 . 0 

297 . 393 

0 . 293 

0 . 2929 

*  0.01101 

178.0 

<0 . 3733 

90  Tension 

90^  Compression 

Shear 

Strain 

Stress 

Strain 

Stress 

Strain 

Stress 

( 1 n/ 1 n ) 

( ks  1 ) 

( in/ i n ) 

( ks  1 ) 

(  1  n  /  1  n  ) 

(ksi  ) 

0.001 

1  .  55 

0 . 002 

3 . 20 

0 . 004 

3 . 090 

0 . 002 

3.10 

0 . 004 

6 . 40 

0 . 008 

5  .  200 

0 . 003 

4 . 60 

0 . 006 

9 . 60 

0.012 

6 . 535 

0.004 

6 . 05 

0 . 008 

12.75 

0.016 

7 . 440 

0 . 005 

7 . 45 

0.01 

15.50 

0 . 02 

8 . 080 

0 . 006 

8 . 85 

0.012 

18.00 

0 . 024 

8.610 

0 . 007 

10.20 

0.014 

20 . 25 

0.028 

9 . 020 

0 . 008 

11.50 

0.016 

22 . 20 

0 . 032 

9 . 330 

0 . 009 

12.55 

0.018 

24 . 05 

0 . 036 

9 . 578 

0.01 

13.30 

0 . 02 

25 .73 

0 . 04 

9.810 

0.011 

13.90 

0 . 022 

27 . 30 

0 . 044 

10.030 

0.012 

14.30 

0 . 024 

28 . 70 

0.048 

10 . 230 

*  0.0129 

14 . 586 

0 . 026 

29 . 90 

*  0.056 

10 , 540 

0 . 028 

30 . 90 

0 . 03 

3  1.80 

0.032 

32 . 50 

0 . 034 

33.10 

*  0.0352 

33 . 40 

*  ■  denotes  U 1 

timate  Stress 

-Strain  Val 

ues  . 

1  29 


F 


Tab-e  .2.  Engineerinjg  Elastic  Constants  for  APC-2  Gr/PEEK 
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specimens,  but  their  results  will  be  presented  in  a  separate 
section.  For  all  tests,  data  was  gathered  in  various  forms 
with  the  goal  of  determining  stress-strain  response  to 
ultimate  strength  of  the  laminates.  The  results  from  this 
phase  of  e.xpe  r  i  men  ta  1 1  on  are  presented  in  this  section.  As 
an  overall  view  of  the  results,  all  5  types  of  laminates  are 
shown  as  failed  specimens  in  Figure  45.  The  results  of  the 
quas 1  -  1  so trop 1 c  study  are  discussed  in  Section  D  of  this 
chapter.  As  described  in  Chapter  III,  the  0T,  90T,  and  £K 
laminates  were  also  analyzed  using  a  fully  nonlinear 
ply-by-ply  plane  stress  finite  element  technique.  These 
analytical  studies  of  the  laminates  are  presented  in  tandom 
with  e  X  p  e  r 1  me  ntal  results. 

1.  -0  i  Cn 1 d 1 r ec t 1 ona i  Specimens.  The  ?eneral 

results  of  6  tests  to  ultimate  strength  are  listed  in 
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Table  14. 

U 1 1 1 ma te  Strength  Tests  of  [0  ] 

1  a 

Spec  1  mens 

Failure 

Fa 1  lure 

Samp  1 e 

X-Sect 

Test 

Load 

Stress 

Number 

Area 

Number 

(lbs.) 

(  p  S  1  ' 

A1 -0T-O-0  1 

0 . 04988 

600  1 

14 , 900  . 

159,100. 

A1 -0T-O-02 

0 . 04988 

6002 

15 , 707  . 

165 , 400  . 

A1 -0T-O-03 

0 . 0498 

6003 

14 . 380 . 

153 , 350 . 

A1 -0T-O-04 

0 . 050  13 

6004 

1 5 . 400  . 

136 , 550 . 

A1 -0T-O-05 

0.05  107 

6008 

13 , 850  . 

133,450. 

A1 -0T-O-06 

0.0498 

6012 

15 , 400  . 

128,250. 

Average  Fail 

ure  Stress  = 

146,000. 

Table  14,  For  these  specimens,  a  repeatable  and  significant 
occurrence  took  place.  At  40*4  of  each  specimen's  failure 
load,  splitting  of  the  laminate  occurred.  It  began  at  the 
sides  of  the  hole  and  suddenly  extended,  parallel  to  the 
fibers,  through  the  tab.  Failure  of  each  speci.men  was 
dominated  by  this  phenomenon.  A  single  specimen  is  shown  in 
Figure  46.  Notice  that  the  top  segment  slipped  out  of  the 
left  tab  at  failure,  and  the  bottom  segment  slipped  out  of 
the  right  tab.  The  middle  segments  were  simply  unloaded. 

'^hen  the  splitting  occurred,  nearly  an  instantaneous 
drop  in  measured  load  occurred.  This  can  be  seen  in 
Figure  47  on  the  1 oad - d i s p 1 acemen t  curve.  Each  of  the 
d 1  scon 1 1 nui 1 1 es  tat  approximately  6000  lbs.;  represents  one 
side  of  the  s p e c i me n  splitting. 
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F 1 gure  46 .  Failed 


Spec  1  men 
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In  analyzing  this  phenomenon,  one  could  consider  the 
sudden  splitting  as  t.he  ultimate  failure  of  the  laminane. 
Thus,  40%  of  failure  load  from  Table  14  could  be .  cons  id'ered 
as  the  actual  ultimate  load  on  the  specimen.  Therefore',- 
after  40%  of  the  ultimate  load,  a  different  problem  eocists. 
The  problem  is  one  of  two  0.4-inch-wide  thin  plates  under,  a 
uniaxial  tensile  load,  with  an  unloaded  0.4-inch  wide 
segment  in  the  center.  As  a  comparison  to  the  unnotched 
0-degree  basic  property  speci me n ,  we  could  calculate  the 
failure  stress  of  one  of  the  seg  m^n t  s .  Ho  we  ver.  These"  out  eT- 
segments ,  with  no  stress  risers  in  their  geometry,  never 
failed  across  their  fibers.  Instead,  they  slipped  out  of-  _ 
their  tabs.  Thus,  direct  comparison  to  ultimate  strength  in 
the  0-degree  direction  was  not  possible. 

When  the  splitting  occurred,  the  strain  gages  at  both'., 
the  hole  and  far  field  read  discontinuous  jumps  in  strain.'. 
For  the  gages  at  the  hole,  three  of  the  four  gages  continued 
to  read  strain  after  the  splitting,  and  these  three  readings 
were  similar  past  40%  of  the  failure  stress.  The  fourch 
gage  most  likely  failed  due  to  the  shock  of  the  splitting.  ; 
These  three  sets  of  readings  were  averaged,  and  the  average 
experimental  curves  are  shown  in  Figure  48.  Again,  note  the 
; 'imps  in  the  readi.ngs.  The  analytical  s  t  r  e  s  s  -  s  t  r  a  i  .n  curve 
IS  also  plotted  agai.nst  the  laboratory  results.  The 
3 t r e s 3  -  3 t r a  1 n  response  was  modeled  accurately  up  to  the 
occurrence  of  splitting,  then  to  a  good  appro xi, nation  bevond 
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•hat  load.  The  analytical  discontinuities  that  occurred 
just  before  the  element  near  the  hole  failed  were  a  result 
of  the  element's  strain  response  reactin;?  to  other  elements 
failing  next  to  it.  The  current  version  of  the  program 
contained  only  a  gradual  unloading  scheme  (see  Section  II 
Q' ,  and  these  strain  responses  were  a  result  of  this  scheme. 
This  unloading  scheme  produced  some  anomalies  between 
experimental  and  analytical  results  for  both  types  of 
unidirectional  laminates,  and  these  will  be  discussed  in 
subsequent  paragraphs.  The  general  trend  was  accurate, 
however,  and  the  model  predicted  failure  of  that  area 
exactly  within  the  range  of  experimental  values. 

After  the  splitting,  most  far  field  gages  continued  to 
read  strain,  but  no  readings  were  repeatable,  as  they  were 
before  the  splitting.  For  six  specimens,  the  far  field 
3 t r ess - s t r a  1 n  response  is  shown  in  Figure  40.  Note  the 
discontinuities.  These  results  are  reasonable  since  the 
gages  were  placed  in  the  center  of  the  specimen,  which 
became  unloaded  after  40%  of  the  ultimate  load.  Also  shown 
on  this  plot  IS  a  superimposed  analytical  response.  Mote 
that  the  computer  program  accurately  measured  the 
3 t r es s - s tr a  1 n  response  until  the  splitting  occurred.  The 
splitting  never  reached  the  element  measuring  the  far  field 
response  because  of  the  gradual  unloading  scheme  used  in  the 
program.  Therefore,  the  element  continued  to  output  stress 
and  strain  values. 


Note  tha 


or  all  the  stress-strain  plots  of  this 


section,  the  axes  for  the  curves  are  the  nominal  stress  on 
the  laminate  (load/area)  versus  the  axial  and  transverse 
strains  in  the  element. 

2.  ^Unidirectional  Laminates  .  Seven  specimens 

were  tested  to  failure,  and  the  general  results  of  these 
tests  are  listed  in  Table  15. 

For  these  specimens,  a  predictable  failure  mode 
occurred.  The  specimens  simply  failed  at  the  most  narrow 
part  of  their  cross -sec t i onal  area.  This  type  of  failure  is 
shown  in  the  photograph  of  Figure  50. 

For  the  strain  gage  near  the  hole,  the  average 
experimental  s tress -s trai n  response  is  shown  compared  to  the 
a.nalytical  response  in  Figure  51.  Again,  the  model 


Table  15. 

Ultimate  Strength  Tests 

Of  [90  ] 

1  <5 

S  p  e  c  1  me  n  s 

Samp  1 e 

X-Sect 

Test 

Failure 

Load 

Fail ur e 
Stress 

.Number 

Area 

Number 

( lbs  .  ) 

(  ps  1  ) 

B2-90T-H-02 

0 . 098  15 

7002 

765  . 

7794  . 

B2-90T-H-03 

0 . 098  15 

7003 

850  . 

8660  . 

B2-90T-H-04 

0 . 098  15 

7004 

780  . 

7947  . 

32-90T-H-05 

0 . 09840 

7005 

308  . 

3211. 

B3-90T-H-08 

0 . 09720 

7008 

780  . 

3  0  2  5. 

34-90T-H- 12 

0 . 09720 

70  12 

7  12. 

7325  . 

34 -90T-H- 16 

0 . 098646 

70  16 

697  . 

7066  . 

Average  Failure  Stress 

=  7860. 
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acci^rai^ely  predicted  stress-strain  behavior,  and  it  exactly 
predicted  the  u 1 1 1 ma te  failure  of  the  area  next  to  the  hole. 
Its  prediction  of  7,551  psi  at  failure  fell  within  the  range 
of  experimental  vaules  of  failure  stress.  Analytically,  the 
strain  response,  ,  follows  closely  the  basic  property 
curve  for  90-degree  tension.  The  actual  area  failed  with  a 
stiffen  response.  The  reason  for  this  variance  probably 
relates  to  the  modeling  of  .gradual  failure  versus  the 
reality  of  sudden  failure.  This  concept  will  be  further 
discussed  in  the  next  section. 

The  s tress - s trai n  response  is  shown  in  Figure  52  for 
the  far  field  gage.  The  analytical  stress-strain  response 
modeled  the  experimental  response  exactly  as  shown. 

However,  beyond  the  experimental  failure  stress,  the  model 
still  showed  some  stiffness,  due  to  the  unloading  scheme. 

And  after  40  load  increments,  which  corresponded  to  a  stress 
of  8840  psi,  the  model  had  not  yet  failed  completely  across 
the  specimen,  and  the  computer  run  was  stopped. 

3,  Angle-Ply  Laminates  .  Seven  angle-ply 

laminates  were  tested  to  their  ultimate  strength,  and  the 
general  test  results  are  listed  in  Table  16.  The  failure  of 
these  specimens  was  characterized  by  hig,h  strain  to  failure, 
especially  in  the  region  of  the  hole.  This  high  strain  is 
apparent  in  the  failed  SH  specimen  shown  in  Figure  53.  Mote 
that  si.gnif  leant  necking  occurred  in  the  region  of  the 
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Table  io.  Ultimate  Strength  Tests  of  [-451  Soecimens 

4S  ■ 


Fa  1  lure 

Fa  1  lure 

Samp  1 e 

X-Sect 

Test 

Load 

Stress 

Number 

Area 

Number 

( lbs  .  ) 

i  p  S  1  ) 

C2-SH-O-0  1 

0 . 09752 

300  1 

2 , 628  . 

26  ,  948  . 

C2-SH-O-02 

0 . 09865 

8002 

2 .669 . 

27 , 055 . 

C2-SH-O-03 

0 . 09952 

8003 

2 . 630  . 

26 , 427  , 

C2-SH-O-04 

0 . 09993 

8004 

2  ,  423  . 

24  ,  247  . 

C2-SH-O-08 

0 . 0972 

8008 

2  ,352 

24,198. 

::2-3H-0-  1  2 

0 . 09856 

80  1  2 

2 , 400  . 

24,351. 

C.3-3H-0-  16 

0 . 09985 

30  16 

2,650. 

26.540, 

■i 

mm 

=  25.630. 

hole.  The  holes  itself  also  appears  as  an  ellipse  in  the 
specimen  of  Fi;?ure  53,  due  to  the  high  deformations. 

For  the  gage  at  the  hole,  the  very  high  strain  which 
occurred  in  these  specimens  started  early  in  the  loading. 

As  shown  in  the  next  section,  a  distinct  failure  region  was 
growing  at  the  sides  of  the  holes.  When  that  region  met  the 
gage,  the  gage  failed.  The  stress  at  failure  varied  for 
each  gage,  but  in  ail  cases  it  was  significantly  below  the 
failure  stress  of  the  entire  laminate.  Taking  the  average 
failure  stress  corresponding  to  the  strain  gage  at  the  hole, 
the  failure  of  the  region  over  which  the  gage  was  placed  was 
about  45 T.  of  the  average  failure  stress  of  the  entire 
specimen.  The  s t ress -s t r a i n  curves  for  the  area  at  the  hole 
were  averaged,  and  this  average  curve  is  plotted  in 
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riiuro  54  as  shear  stress  on  the  specimen  versus  shear 
strain  at  the  hole. 

Ana  1 y 1 1  cal  1 y ,  the  s tress -s trai n  response  and  ultimate 
failure  of  the  element  matched  the  actual  response  very 
well.  This  analytical  response  is  also  plotted  in  Figure 
54.  Note  that  the  first  analytical  s t r ess -s tr a i n  values 
were  vt'tpu'^  almost  at  the  u.lti.ma.te  failure  stress  ::r  the 
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the  hole  was  suddenly  very  large.  If  a  smaller  load 
increment  was  applied,  the  analytical  s t r ess -s t ra i n  curve 
wou.f  have  better  matched  the  experimental  response. 

Note  also  in  Figure  54  that  the  analytical  ultimate 
failure  stress  predicted  very  closely  the  failure  of  the 
area  near  the  hole.  Thus,  cutting  off  the  shear  basic 
property  data  at  values  corresponding  to  approximately  5.5% 
strain  yielded  accurate  results. 

The  3 tr es s - 3 tr a  1 n  curves  of  nominal  shear  stress  versus 
the  far-field  shear  strain  (as  defined  in  Section  IV  E)  is 
shown  in  Figure  55.  Note  that  repeatable  data  was  obtained 
except  for  one  errant  strain  reading.  The  analytical 
stress-strain  behavior  matched  this  experimental  response 
very  well,  also  shown  in  Figure  55;  and  this  result  met  the 
4  0  a  I  of  this  test. 
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Overall,  the  computer  analysis  provided  good  results  in 
modeling  the  s tr e s s - s tr a i n  behavior  of  the  0T ,  90T  ,  and  3H 
laminates.  Conclusions  will  be  drawn  for  the  results  of 
this  section  in  Chapter  VI. 

C  .  Results  o  f  the  Failure  Progress  ion  Study  for  ^  ^  ^  J  • 

[90  ],  and  [i45]  Laminates 

1  <S  -  4S  - 

The  results  of  studying  the  initiation  and  progression 
of  failure  in  the  laminates  we  re  derived  from: 

,1)  A  videotape  post-failure  analysis, 

\2':  A  stereo  x-ray  post-failure  analysis,  and 

'3)  A  progressive-failure  finite  element  technique. 
Mote  that  the  specimens  were  subjected  to  x-ray  only  after 
xoading.  They  were  not  x-rayed  before  loading  to  check  for 
fabrication  flaws. 

1  .  Fai lure  Progression  i  n  C  0  _  ^ ]  Unidirectional 
Specimens .  As  discussed  in  the  previous  section,  splitting 
occurred  for  these  specimens  at  40/i  of  their  ultimate 
strength.  This  phenomenon  thus  dominated  the  experimental 
stress-strain  behavior  of  this  laminate. 

The  splitting  and  eventual  failure  was  documented  by 
the  video  camera,  and  4  stills  were  taken  for  a  represent.a- 
tive  specimen.  Figure  56  shows  the  tensile  coupon  for  nest 
*  600C  at  the  beginning  of  the  loading.  For  this  test,  one 
volt  corresponded  to  4,082  pounds.  The  splitting  occurred 
at  a  reading  of  about  1.60  volts.  After  further  loading, 
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■  h  ■?  s  "  1  :  '  IS  3  1  1 h  n  1  V  p  e  r  c  e  ri  b  1  ■=  at  3.50  •/  ?  1  t  s  1  4  .  7 
lbs.'  :n  Figure  57.  The  ulti  ma  te  load  on  the  spec:  me  r.  :  s 
she  wn  in  Figure  53  at  4.05  volts  (16,500  lbs.),  and  the  next 
frame  of  the  videotape  provided  the  failed  laminate,  as 
shown  in  Figure  59.  implying  sudden  failure.  The  '-0.02' 
indicates  a  sudden  zero  loading  on  the  specimen. 

The  results  of  all  tests  to  percentages  of  the  average 
failure  stress  are  shown  in  Table  17.  (The  average  failure 
stress  was  146,000;  see  Table  14.)  A  stereo  x-ray  was 
taker,  of  all  3  specimens.  As  expected,  the  only  failure 
region  that  is  visible  is  the  splitting  along  the  fiber 
direction.  The  effect  of  stereo  viewing  provided  the 


Table  17.  Progressive  Failure  Tests  of  [0  ]  Specimens 

1  a 


Samp  1  e 
Number 

X-Sect 

Area 

Test 

Number 

♦/ 

/♦ 

Load 

(kips) 

Actual 

♦/ 

A2 -0T-H-05 

0  .  10  114 

6005 

50% 

7.417 

""3 . 334 

50 . 2 

A2-0T-H-06 

0 . 0996 

6006 

50 

6 . 997 

70.251 

48  .  1 

A2-0T-H-07 

0 . 09993 

6007 

50 

7  .  598 

76.033 

5  2.1 

A3-0T-'d-09 

0 . 09504 

6009 

30% 

11.133 

117.665 

30.6 

A3-0T-H-  10 

0 . 096  16 

60  10 

80 

11.228 

1  Id  .  "^64 

30 . 0 

A3-0T-H  -  1  1 

0 . 09832 

60  1  1 

80 

1  2 . 925 

131.458 

90 . 0 

A3-0T-H-  13 

0 . 09728 

60  13 

90% 

12.  126 

124 . 650 

♦  # 

A3-0T-H- 1 4 

0 . 09632 

60  1  4 

90 

12 . 044 

125.042 

35 . 6 

This  13  the  percentage  where  the  test  machine  was 
attempted  to  be  stopped. 


This  specimen  broke  unexpectedly. 
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]  Specimen  at  Maximum  Load  of  16,500  lbs 


]  Specimen:  Sudden  Failure;  Instant  Unloading 


n 


expected  result  that  the  splitting  occurred  perpendicular  to 
the  front  face  of  the  specimen  and  through  the  thickness.  A 
typical  specimen  stressed  almost  to  failure  is  shown  in 
Figure  60.  The  stereo  x-rays  for  the  failed  specimen  of  test 
*  6013  are  shown  in  Figure  61. 

For  an  analytical  comparison  to  growth  of  failure, 

Figure  62  shows  a  similar  progression  to  the  experimental 
results.  As  shown,  the  failure  began  at  the  element  that 
was  .most  highly  stressed  and  proceeded  parallel  to  the  fiber 
direction.  However,  the  analytical  result  shows  a  failure 
region  developing,  instead  of  a  single  line  of  elements 
failing.  The  reason  for  this  difference  is  in  the  unloading 
scheme  of  the  computer  program.  The  scheme  in  the  current 
version  of  program  provides  only  for  gradual  unloading, 
although  the  specimen  modeled  was  a  unidirectional  element. 
'See  Section  II  G.)  Thus,  as  the  elements  failed,  they  were 
still  modeled  to  have  stiffness.  (The  program  is  currently 
being  developed  to  handle  instantaneous  unloading.)  The 
results  of  Figure  62,  however,  do  show  an  accurate  trend  of 
failure  progression. 
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va.ue  of  55%  of  ultimate  load  would  decrease  to  a  value  near 
40%.  which  13  our  experimental  result.  Thus,  the  onset  of 
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failure  in  the  0-degree  specimen  could  be  predicted  in  this 
ma nner.  From  this  analysis,  the  splitting  of  the  entire 
laminate  could  be  defined  by  the  failure  of  the  first 
e  1  e  me  n  t  . 

Thus,  analytically,  this  problem  has  two  bounds:  (1; 

The  definition  of  sudden  splitting  with  failure  of  the  first 
element,  which  models  the  actual  failure,  and  (2)  the 
gradual  failure  of  the  elements,  which  does  not  model 
reality  but  is  an  academic  exercise.  No  interlaminar 
stresses  exist  in  the  0'degree  laminate,  and  this  gradual 
failure  scheme,  or  strain  softening  concept,  does  not  apply. 

a 

2 .  Fail ure  Progression  i n  90  Unidirectional  Specimens . 
As  implied  in  the  previous  section,  the  failure  of  these 
specimens  was  dominated  by  matrix  failure  beginning  at  the 
hole  and  extending  to  the  edge  of  the  laminate.  However, 
the  progress  ion  of  this  failure  was  barely  noticeable  in  the 
post-failure  analyses. 

From  the  videotape,  three  loading  'scenes'  of  a  typical 
test  are  shown  in  Figures  63  through  65.  In  Figure  63,  a 
typical  90-degree  specimen  is  shown  at  the  beginning  of 
loading.  For  these  tests,  one  volt  corresponded  to  a  load 
of  408  pounds.  In  Figure  64,  where  the  specimen  is  loaded 
to  nearly  100*4  of  its  failure  load,  no  failure  region  is 
shown.  Then,  after  advancing  one  frame  of  the  videotape, 
where  one  frame  equals  1/60  sec)  the  video  picture  shown  in 
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Specimen, 


Figure  65  rendered  a  failed  specimen,  again  implying  sudden 
: a  1  lure . 

The  results  of  each  test  to  a  percentage  of  the  average 
failure  stress  (7860  pounds)  are  presented  in  Table  18. 

These  tests  and  the  stereo  x-rays  provided  further  proof  of 
sudden  failure.  A  specimen  loaded  to  approximately  907.  of 
the  average  failure  load  showed  little  sign  of  being 
stressed.  In  Figure  66  the  white  regions  shown  could 
represent  the  beginning  of  a  failure  area,  or  they  could 
simply  be  a  hole  defect  caused  by  the  drill.  The  specimen 
shown  in  Figure  57  has  failed  through  the  thickness  but  not 
exactly  at  the  hole.  This  unexpected  failure  could  have 
been  initiated  by  flaws  or  notches  from  the  fabrication  of 
the  hole. 


Table  18 

Progressive  Fai 

lure 

Tests  of 

[90^^"  Specimens 

Sample 

X  -  S  e  c  t 

Test 

Load 

Stress 

Actual 

Mumbe  r 

Area 

Mumber 

7/ 

(lbs.) 

( ps  1 ) 

•/ 

B2-30T-H-0  1 

0 . 09843 

700  1 

307. 

639  . 

6489  . 

32.6 

33-90T-H-0e 

0 . 09t52 

7006 

30 

625  . 

5408  . 

3  1.5 

33-90T-H-07 

0.09744 

7007 

30 

632  . 

5486  . 

32.5 

33 -90T-H-09 

0 . 09744 

■’009 

907. 

697  , 

'"153. 

9  1.0 

53-90T-H-  10 

0 . 096 

70  10 

90 

694  , 

7229  . 

92 . 0 

33 -90T-H- 1 : 

0 . 09624 

70  1  1 

90 

695  . 

-f  n  0  0 

A,  <,  - 

9  1.9 

34 -90T-H- 13 

0 . 09752 

70  13 

9  57. 

737  . 

■"557  . 

96  .  2 

34 -90T-H- : 4 

0.09865 

70  1  4 

95 

746  . 

7562  . 

♦  ♦ 

34 -90T-H- 1 5 

0 . 0984 

70  15 

95 

735  . 

7470  . 

95.0 
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90jg]  Specimen  Stressed 
e  Stress . 


67.  Stereo  X-Ray  of  Failed  [90 


Analytically,  the  progression  of  failure  modeled  by  the 
computer  program  is  shown  in  Figure  68.  Again  a  region  of 
failure  developed  instead  of  the  failure  of  a  row  of 
ele.ments  across  the  specimen  and  parallel  to  the  fibers.  As 
in  the  0-degree  tests,  the  computer  program  used  a  gradual 
unloading  scheme  on  this  model.  Thus,  only  the  trend  of 
failure  progression  is  accurate. 

As  in  the  previous  subsection,  failure  in  the 

analytical  model  could  be  defined  by  failure  of  the  first 

element.  In  this  case,  failure  occurs  in  the  element  at  the 

hole  at  89’/.  of  the  average  ultimate  stress.  Extrapolating 

to  the  hole,  the  percentage  is  further  reduced  to 

approximately  36%.  This  prediction  of  sudden  failure  is 

lower  than  the  actual  result.  However,  the  fibers  in  the 

•» 

90-degree  specimens  may  not  have  been  exactly  at  90  ,  thus 

raising  the  ultimate  stress  value.  The  stress  concentration 
factor  drops  off  significantly  if  the  fibers  are  not  at  90 
(See  Figure  9  in  Section  II  D.) ,  and  this  may  account  for 
the  difference  in  results. 

3.  Failure  Progress  ion  i  n  Angle-Ply  Specimens  . 

As  mentioned  in  the  last  section,  a  distinct  failure  region 
grew  at  the  hole  during  the  loading  of  t45-degree  specimens. 
Physically,  strain  response  was  characterized  by  a  necking 
oi  the  specimen  at  the  hole  and  the  growth  of  a  triangular 
region  on  either  side  of  the  hole.  It  see  me  d  that  wh  e  n 
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‘hese  'triangles'  reached  the  midpoint  between  the  hole  and 
the  edge,  the  specimen  failed. 

This  progression  o.f  failure  is  shown  in  several 
photographs  from  the  videotape.  For  a  reference.  Figure  69 
shows  the  specimen  at  the  beginning  of  loading.  Note  here 
that  1  volt  equals  1020  pounds.  Note  also  that  the 
displacement  of  the  specimen  can  be  measured  qualitatively 
by  loo 'King  at  the  wr  iting  on  the  speci  me  n  at  the  left  of  the 
picture.  The  crosshead  is  moving  left  and  pulling  with  it 
the  left  end  of  the  specimen.  The  writing  on  the  specimen 
IS  slowly  moving  out  of  view. 

The  next  five  photographs  show  a  senes  of  failure 
states.  Figure  70  shows  the  initial  region  of  failure  and 
the  initial  necking  of  the  specimen.  Figure  71  shows  a 
distinct  triangle,  and  Figure  72  shows  the  development  of 
reflective  cracking  toward  the  edge  of  the  specimen.  This 
seems  to  be  a  three-dimensional  effect  caused  by  the 
interlaminar  stresses  at  the  edges  of  the  specimen.  As  will 
be  shown,  the  computer  program  did  not  model  this 
three-dimensional  phenomenon.  Continuing  the  loading, 

Figure  73  represents  the  maximum  load  (2,660  lbs.)  on 

th 1 sspec 1  men .  Finally,  Figure  7  4  shows  the  specimen  failure 

after  the  load  dropped  off  to  2670  lbs. 

The  results  of  the  tests  to  percentages  of  average 
failure  stress  (25,680  lbs.)  are  listed  in  Table  19.  The 
stereo  x-rays  from  each  of  these  tests  yield  results  si. mil. a r 
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Specimen  at  2540  lbs.  (957.  of  Ultimate). 


7 


Table  19. 

Progress 

ive  Failure  Tests  of  [ 

t45 ]  Spec  1 

4S  ^ 

mens 

Samp  1 e 

X-Sect 

Test 

Load 

Stress  Actual 

Number 

Area 

Number 

7.* 

(kips) 

( ks  1  ) 

*/ 

C2-SH-H-05 

0 . 09848 

8005 

70% 

1  .  773 

1 8 . 004 

70  .  1 

C2-SH-H-06 

0 . 09873 

8006 

70 

1  .  750 

17 . 725 

69 . 0 

C2~SH-H-07 

0 . 09848 

8007 

70 

1  .  755 

17.821 

69 . 4 

C2-SH-H-09 

0 . 095  12 

8009 

85% 

2 . 076 

21.825 

35 . 0 

C2-SH-H- 10 

0 . 096 

8010 

85 

2 . 095 

2  1  .  823 

85.0 

C2-SH-H- 1 1 

0 . 09632 

30  1  1 

85 

2.102 

21.823 

85  .  0 

C3-SH-H- 13 

0 . 09744 

80  13 

95% 

2.377 

24 . 394 

94 . 3 

C3-SH-H- 1 4 

0 . 09728 

80  14 

95 

2 . 373 

24 . 394 

94 . 3 

C3-SH-H- 1 5 

0 . 097  1  2 

8015 

95 

2 . 369 

24 . 393 

95.0 

to  those  of  the  videotape,  but  with  some  new  insights. 

Figure  75  shows  an  x-ray  of  a  specimen  loaded  to  707, 
of  Its  probable  failure  load.  The  beginning  of  necking  is 
evident,  and  the  hole  appears  slightly  elliptical.  Also, 
the  beginning  of  the  triangular  region  of  failure  is 
slightly  visible.  Figure  76  shows  an  x-ray  of  a  specimen 
loaded  to  857.  of  its  failure  load,  and  the  failure  region  is 
more  evident.  Figure  77  shows  a  specimen  close  to  failure. 

The  stereo  effect  of  viewing  these  pairs  of  x-rays 
showed  that  the  tips  of  the  triangle  exist  in  the  middle  of 
the  specimen  lay-up.  Close  to  the  hole,  failure  is 
occurring  through  the  thickness,  but  it  t.ipers  off  to  the 
center  of  the  iaminas.  Also,  this  failure  region  does  not 
appear  to  be  that  of  de 1  ami nat i on ,  but  instead  each  line 
shown  13  series  of  short  parallel  cracks  along  the  fibers. 
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Figure  75.  Stereo  X-Ray  of  [±45]  Specimen  Loaded  to  70% 
of  its  Ultimate  Load. 
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Figure  77.  Stereo  X-Ray  of  [±45]^^  Specimen  Loaded  to  957. 
of  its  Ultimate  Load. 


Further  evident  in  these  x-rays  were  minute  cracks 
extending  along  the  fibers  from  the  hole  to  the  edge  i 1  the 
specimen.  These  occurred  in  various  layers,  but  they 
apparently  did  not  degrade  the  laminate.  Also  evident  in 
these  x-rays  is  that  the  edges  appear  to  be  failing  where 
the  necking  is  occurring. 

For  the  computer  model  of  the  failure  process,  a 


31  m.  liar  region 

of  failure 

was  shown . 

This  analytical 

progression  of 

failure  is 

shown  in  F 1 

g  u  r  e  7  8. 

Mote  t  h  a 

triangular  regi 

0  n  s  s  h  0  wn 

in  the  video 

3  t/  1  1  3 

and  in  the 

stereo  x-rays  are  not  reflected  in  the  analytical  failure 
progression.  A  possible  explanation  for  the  V-shape  of  the 
analysis  is  that  in  the  experimental  specimen,  the 
* hree-d 1  mens lonal  edge  effects,  caused  by  interlaminar 
stresses,  restrict  the  V-shape  from  occurring;  and  a 
triangular  region  is  formed.  Thus,  the  computer  model 
provides  an  accurate  account  of  the  failure  process  in  a 
two-dimensional  format,  but  it  cannot  measure 
three-dimensional  effects. 

D .  Results  o  f  the  Quasi-isotropic  Study 

As  described  in  Chapter  IV,  the  two  types  of  quasi- 

isotropic  laminates,  [ 0 /  + 4 5 / 90 / - 45 ]  and  [0/±45/90j  .  were 

zs  zs 

tested  for  experimental  comparison  of  s t r e s s - s t r a  1 n 
responses.  Also,  an  experi me ntal  investigation  of  the 
initiation  and  progression  of  ply  failure  was  conducted. 
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A;§ain,  the  purpose  of  this  test  was  to  compare  experimental 
behavior  of  the  laminate  and  to  draw  conclusions  about  their 
potential  for  delamination  based  on  comparisons  to  similar 
Gr/Ep  specimens  in  Reference  C23]. 

1.  Results  of  Tests  to  Ultimate  S  tren^j  th .  The  results 
of  all  tests  to  failure  of  the  [ 0/ +  45/90/ -45  ]  (Ql)  and 
[0/X45/90J  (Q2)  tensile  coupons  are  listed  respectively  in 

Tables  20  and  21.  Note  that  the  average  failure  stresses 
are  identical,  considering  the  range  of  all  failure 
stresses.  A  failed  specimen  of  each  laminate  is  shown  in 
Figures  79  and  80.  Note  here  that  the  failure  surfaces  are 
also  very  similar.  The  only  difference  seems  to  be  that  the 
Ql  specimens  (designated  Dl)  show  a  more  ±45-degree 
dominated  failure;  and  alternatively  the  Q2  (El)  specimens 
show  a  more  90-degree  dominated  failure. 

This  difference  in  failure  surfaces  proved  to  be 
negligible  when  the  stress-strain  curves  of  the  two 
quas 1 - 1  so trop 1 c  specimens  were  compared.  Figures  81  and  82 
show  the  far  field  response  to  tensile  loads  for  the  Ql  and 
Q2  specimens,  respectively.  Note  that  if  these  curves  were 
mapped  together  there  would  be  no  discernable  difference  in 
the  responses . 

Figures  83  and  84  show  the  experimental  s tress - s trai n 
response  at  the  holes  of  the  Ql  and  Q2  specimens.  Note 
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Tab  1 e  20 .  Fa i 

lure  Tests 

of  [0/+45/90/-45:„ 

2S 

Specimens 

Samp  1 e 

X-Sect 

Test 

Fai lure 
Load 

Fai lure 
Stress 

Number 

Area 

Number 

( lbs  .  ) 

( ps  1  ) 

Dl-Q-H-0 1 

0 . 09383 

9001 

3 , 986  . 

42 ,480 . 

Dl-Q-H-02 

0 . 09255 

9002 

4 , 000  . 

43 . 220  . 

Dl-Q-H-03 

0 . 09472 

9003 

4 ,240  . 

44 .760 . 

Dl-Q-H-04 

0 . 09504 

9004 

4  ,  160 

43 . 770 . 

Dl-Q-H-08 

0 . 09352 

9008 

3 , 964  . 

42 . 390  . 

Dl-Q-H- 12 

0 . 09376 

9012 

4 , 084 . 

43 . 560 . 

D2-Q-H- 16 

0 . 0939  1 

9016 

4 , 292  . 

45 . 700 . 

Average  Failure  Stress  =  43,700. 


Table  21. 

Failure  Tests 

of  [0/1 

45/90]  Specimens 

Sample 

X-Sect 

Test 

Fai lure 
Load 

Fai lure 
Stress 

Number 

Area 

Number 

(lbs  .  ) 

(psi) 

El-Q-H-01 

0 . 09383 

120  1 

3,714. 

39 , 580 . 

Ei-Q-H-02 

0 . 09496 

1202 

4 , 265  . 

44,910. 

El-Q-H-03 

0 . 09584 

1203 

4,391. 

45 ,820 . 

E  1  -Q-H-04 

0 . 09632 

1204 

4 . 330  . 

44 ,950  . 

El -Q-H-08 

0 . 09352 

1208 

3 . 936  . 

42 , 090  . 

El-g-H- 12 

0 . 09496 

12  12 

4 .020  . 

42 . 330 . 

E2-Q-H- 16 

0 . 09391 

1216 

4,154. 

44 , 230 . 

Average  Failure  Stress  =  43,420. 


(By  assuming  that  test  *1201  is  a  bad  test,  the  average 
failure  stress  would  be  44,050  psi.) 


Figure  79.  Failed  [ 0/45 /90/ -45 ] 2g  Specimen 


Figure  80.  Failed  [0/±45/90}  g  Specimen 


C0Uf>3NS  0-0- 


(ZE^POkIe  at  FaIZ  FIEU?  fc» 

^p6ciMef4^ . 


COUPON 


again  that  the  responses  of  each  laminate  are  nearly 
identical . 

2,  Resul ts  o  f  Fai lure  Progress  ion  Study .  For  the 
quasi-isotropic  laminates  tested  to  ultimate  strength,  the 
videotape  post-failure  analysis  was  conducted.  Also,  a 
stereo  x-ray  post-failure  analysis  was  conducted  on 
specimens  subjected  to  percentages  of  their  average  failure 
stresses  . 

Typical  photographs  were  taken  off  the  videotape.  For 
the  Q1  specimen,  a  series  of  failure  states  is  shown  in 
Figures  85  through  88.  Likewise,  a  group  of  photographs  is 
shown  in  Figures  89  through  92  that  represents  progressive 
failure  of  the  Q2  specimen.  For  both  tests,  the  voltmeter 
read  1  volt  for  every  2,041  pounds.  Notice  that  failure  for 
both  specimens  does  not  begin  to  show  through  the  specimen 
until  late  in  the  loading. 

Q1  and  Q2  laminates  were  tested  to  percentages  of 
their  ultimate  loads,  according  to  Chapter  IV,  and  the 
results  of  these  tests  are  listed  in  Tables  22  and  23. 

Each  specimen  tested  was  subjected  to  stereo  x-ray  and 
typical  results  are  shown  in  Figures  93  through  95  for  the 
Q1  specimen  and  in  Figures  96  through  98  for  the  Q2 
specimen.  For  the  Q 1  specimen,  at  80%  of  the  average 
failure  load,  only  a  few  45-degree  cracks  were  evident,  and 
these  were  mostly  in  the  middle  plane.  See  Figure  93. 
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Figure  86.  e0/45/90/-45 Specimen  at  4020  lbs.  (95%  of 

Ultimate  Load)  . 
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Figure  87.  [ 0 / 45/90 / -45 ] Specimen  at  4200  lbs.  (99%  of 

U1 timate  Load) . 
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Samp  1 e 
Number 

X-Sect 

Area 

Test 

Number 

7. 

Load 

(kips) 

Stress 
( ks  1 ) 

Actual 

% 

Dl-Q-H-05 

0 . 09368 

9005 

80% 

3.317 

35 . 408 

81.0 

D1 -Q-H-06 

0. 09512 

9006 

80 

3 . 368 

35 . 408 

8  1.0 

Dl-Q-H-07 

0 . 09616 

9007 

80 

3 . 374 

35 . 087 

80 . 3 

Dl-Q-H-09 

0 . 09263 

9009 

90% 

3.651 

39.415 

90 . 2 

Dl-Q-H- 10 

0 . 0936 

90  10 

90 

3 . 702 

39.551 

90 . 5 

D 1 -Q-H- 1 1 

0 . 09383 

9011 

90 

3 . 697 

39.401 

90 . 2 

D2-Q-H- 13 

0 . 09015 

90  13 

9  5% 

3 . 729 

4  1  .  364 

94 . 7 

D2-Q-H- 14 

0.09143 

90  14 

95 

3 . 809 

4  1  .  660 

95 . 3 

D2-Q-H- 15 

0 . 09232 

90  15 

95 

3 . 838 

4  1  .  573 

95  .  1 

Table  23.  Progressive  Failure  Tests  of  Q2  Specimens 


Samp  1 e 
Number 

X-Sect 

Area 

Test 

Number 

% 

Load 

(kips) 

Stress 
( ks  1 ) 

Ac  tua 1 
% 

El-Q-H-05 

0 . 096 

1205 

80% 

3.317 

34 . 552 

79 . 6 

El-Q-H-06 

0 . 09624 

1206 

80 

3 . 368 

34 . 996 

80 . 6 

El -Q-H-07 

0 . 09616 

1207 

80 

3 . 374 

35 . 087 

80 . 8 

El-Q-H-09 

0 . 09383 

1209 

90% 

3,651 

38.911 

89 . 6 

El -Q-H- 10 

0 . 09368 

1210 

90 

3.702 

39.518 

91.0 

El-Q-H- 1 1 

0 . 09616 

121  1 

90 

3 . 69'^ 

38 . 447 

88 . 5 

E2-Q-H- 13 

0.09143 

1213 

95% 

3.729 

40 . 785 

93 . 9 

E  2  -  Q  -  H  -  1  4 

0.09135 

12  14 

95 

3 . 809 

4  1  . 697 

96 . 0 

E2-Q-H- 15 

0 . 09232 

1215 

95 

3 . 838 

4  1  .  573 

95 . 7 
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-Ray  of  [0/±45/90]2g  Specimen 
late  Failure  Load. 


Further  loading  to  90%  produced  a  slight  failure  region  as 
shown  in  Figure  94.  Upon  magnification,  the  entire  hole 
showed  cracks  all  the  around  which  resembled  the  rays  of  the 
sun  in  90  and  45  degree  directions  (with  some  0-degree 
cracking  as  well) .  Still  further  loading  to  95%  produced 
the  result  shown  in  Figure  95.  This  region  appeared  to  be 
dominated  by  short  cracks  in  the  matrix  in  the  ±45-degree 
and  90-degree  directions  through  the  thickness.  Another 
interesting  note  is  that  whenever  a  crack  extended  out  from 
the  hole,  it  was  criss-crossed  by  short  cracks  of  other 
1  ami nas . 

For  the  Q2  specimen  loaded  to  80%  of  its  failure  load, 
only  slight  cracks  were  evident  all  the  way  around  the  hole, 
along  with  a  few  single  long  cracks  in  the  90-degree  and  / 
d45-degree  directions.  A  specimen  loaded  to  90%  of  the 
average  failure  load  began  to  show  significant  cracking  in 
the  90-degree  direction  on  one  side  only,  and  it  appeared  in 
the  back  of  the  specimen.  A  specimen  further  loaded  to  95% 
shows  developed  failure  regions  on  both  sides  of  the  hole 
and  through  the  thickness  for  all  laminas. 

As  a  final  note  for  the  x-rays  of  t.he  quas  i  -  i  so  tr  op  i  c 
study  (as  well  as  for  those  of  the  other  laminates)  failure 
of  the  composite  for  most  specimens  seemed  to  start  at  the 
bottom  of  the  specimen.  The  specimens  were  fabricated, 
drilled,  and  gaged  consistently:  therefore  there  .must  have 

been  consistent  damage  to  the  hole  to  produce  consistent 
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failure  i 
and  other 
d 1 scussed 
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itiation  at  the  bottom  of  the  hole.  More  on 
points  of  the  failure  progression  studies  wi 
next  in  Chapter  VI. 


this 
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Cone  1  us  ions 


VI  . 

Based  on  the  results  of  this  thesis  and  on  co mp a r i s o n s 
to  past  studies,  many  conclusions  can  be  drawn.  Much  was 
learned  during  the  course  of  this  thesis,  and  not  all  of  it 
was  directly  related  to  the  objectives  first  listed  in 
Chapter  I.  However,  only  conclusions  relating  to  those 
objectives  will  be  discussed  in  this  chapter,  and  the 
appendices  that  follow  will  present  other  information  gained 
in  this  study. 

A.  Cone  1  us  1 ons  on  the  Per i vat  ion  o  f  Bas i c  ?r ooer  1 1  e  s 

The  following  conclusions  were  noteworthy  for  this 
phase  of  experimentation: 

1.  The  material  response  of  Gr/PEEK  is  highly 
nonlinear,  especially  in  shear. 

2.  The  strains  measured  in  the  shear  specimen  far 
exceeded  the  limits  of  linear  strain  theory.  Large 
displacements  took  place.  But  these  displacements  were  a 
combination  of  the  ductility  of  the  matrix  and  the 
scissoring  of  the  fibers.  As  the  fibers  turned  from 
45°toward  0*  the  specimen  gained  strength.  Therefore,  a 
point  on  the  original  experimental  shear  stress-strain  curve 
between  the  small  strains  and  the  scissoring  was  assumed  to 
equal  failure  of  the  laminate.  The  point  that  was  chosen 
yielded  very  good  results. 
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3.  Ideally,  a  convergence  study  should  be  done  to 
correlate  the  cut-off  point  of  shear  s tress -s tra in  basic 
property  values  to  the  ultimate  strength  of  the  ±45-degree 
tension  specimen  with  the  hole.  This  procedure  would  use  the 
experimental  results  of  the  ±45-degree  test  specimen  with 

the  hole  and  relate  the  analysis  with  a  given  range  of 
cut-off  points  in  the  data. 

4.  Standard  strain  gages,  which  were  accurate  up  to 
approximately  4%,  were  used  for  this  thesis.  For  Gr/PEEK 
basic  property  specimens,  "high  elongation'  gages  should  be 
used  to  measure  strain  on  test  specimens.  For  the  shear 
specimens,  the  gages  failed  before  the  specimen.  And  for  a 
high  temperature  study,  given  the  nature  of  thermoplastics, 
much  more  strain  will  develop  in  all  basic  property 
specimens  during  testing. 

5.  Comparisons  to  Gr/Ep  [19]  were  made  and  the  basic 
properties  of  Gr/PEEK  were  found  to  be  similar  in  the 
0-degree  direction,  about  2  times  more  ductile  in  90-degree 
tension,  similar  in  90-degree  compression,  and  more  than  2 
times  more  ductile  in  shear. 

6.  The  basic  properties  of  0-degree  tension  were 
dominated  by  the  fibers,  as  expected.  The  basic  property 
curve  was  concave  up.  This  ‘wool  thread  effect  is  a 
property  of  graphite  fibers,  and  it  was  shown  in  previous 
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studies  as  well.  (The  wool  thread  analogy  is  referenced  to 
a  remark  by  Sandhu,  in  a  conversation  with  him.) 

7.  As  a  final  conclusion  on  the  properties  of 
Gr/PEEK,  to  accurately  predict  an  elastic  response  of 
Gr/PEEK,  a  nonlinear  form  of  the  stress-strain  constitutive 
relations  must  be  used. 

f 

a .  Conclusions  on  the  Behav i or  o f  Gr/PEEK  Unidirectional 
and  An;^le-Pl  V  Lami nates  . 

For  the  [0  ],  [90  ].  and  Ci45]  specimens,  several 

i  <S  1  c> 

conclusions  can  be  drawn; 

1.  Experimentally,  the  stress-strain  responses  of  the 

[0^^]  test  specimens  were  dominated  by  the  discontinuity  of 
splitting  along  the  fibers.  These  specimens  failed  in  two 
distinct  parts:  (a)  sudden  splitting  longitudinally  from 

the  hole  to  the  ends  of  the  specimen  and  (b)  sudden  failure 
apparently  involving  tab  failure  of  the  two  outside  segments 
of  the  specimen. 

2.  Computer  analysis  of  specimens  matched  the 

experimental  s tr ess -s t r a i n  response  up  to  the  splitting  of 
the  specimen.  The  initiation  of  failure  was  accurately 
predicted.  The  progression  of  failure,  however,  was 
overestimated  due  to  the  use  of  a  gradui»l  unloading  scheme 
which  attempts  to  model  interlaminar  stresses.  Since  no 
interlaminar  stresses  exist  in  a  unidirectional  laminate. 
Therefore,  to  model  reality,  the  failing  elements  should 
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have  been  .i.iven  zero  stiffness  instantly-  The  trend  of 
failure  was  demonstrated,  but  the  ultimate  failure  of  the 
model  would  have  occurred  quickly  with  an  instantaneous 
unloading  scheme.  (This  program  is  still  developmental; 
the  instantaneous  unloading  option  will  be  available  soon  for 
f uture  work . ) 

3.  Experimentally,  the  stress-strain  response  of  the 

[90  1  specimens  was  almost  linear,  and  it  culminated  with 

1  O 

catastrophic  failure  in  the  transverse  direction. 

4.  Computer  analysis  of  the  specimens 

accurately  modeled  the  stress-strain  response.  Initiation 
of  failure  was  slightly  underestimated,  but  small  deviations 
in  the  90-degree  orientation  may  have  affected  the 
correlation.  Progression  of  failure  was  overestimated  as  in 
the  0-degree  specimens  due  to  the  gradual  unloading  scheme. 
However,  the  trend  of  failure  was  demonstrated. 

5.  Experimentally,  the  failure  of  is  gradual, 

and  It  IS  one  that  can  be  observed  as  failure  progresses. 

It  IS  characterized  by  high  strains  to  stress.  Shear 
stresses  contribute  to  almost  all  of  the  failure  in  the 
specimen , 

6.  The  stress-strain  response  of  the 
specimens  was  modeled  accurately.  Both  initiation  and 
progression  of  failure  in  the  speci.men  was  predicted  very 
well.  The  point  of  cut-off  for  the  basic  properties  proved 
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to  be  accurate.  However,  three-dimensional  effects  that 
were  present  in  the  experiment  were  not  modeled  by  the 
analysis.  The  use  of  the  gradual  unloading  scheme  works 
well  with  this  angle-ply  laminate. 

C .  Cone  1  us  1 ons  on  the  Quas i-isotropic  S tudy 

Based  on  the  results  presented  in  Chapter  V,  the 

t,endency  of  the  [0/±45/90  ]^^  Qr/PEEK  specimens  to  delaminate 

13  not  as  great  as  that  of  Gr/Ep.  In  Reference  [23],  the 

delamination  moment  coefficient  vDMC)  for  this  lay-up  was 

critical  for  Qr/Ep,  due  to  the  existance  of  two  90-degree 

plies  in  the  middle  of  the  lay-up.  The  DMC  for  the 

[ 0 / 4 5 / 90 / - 45  ]  for  Gr/Ep  was  far  from  critical  due  to  the 

zs 

45-degree  plies  in  the  center.  But  considering  that  both 
lay-ups  of  Gr/PEEK  produced  the  same  stress-strain  behavior, 
the  tougher  PEEK  material  is  more  resistant  to  delamination. 

D .  General  Cone  1  us  ions 

1.  A  good  method  for  the  prediction  of  the  onset  of 
failure  in  composite  laminates  is  Sandhu’s  fully  nonlinear- 
elastic  finite  element  program  with  its  Strain  Energy 
Failure  Criterion.  The  prediction  of  the  onset  of  failure 
proved  to  be  accurate  in  this  study. 

2.  For  prediction  of  growth  of  failure  this  program 
yields  an  accurate  trend. 
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3.  An  instantaneous  unloading  scheme  should  be 
considered  for  unidirectional  specimens. 

4.  A  gradual  unloading  scheme  works  well  for  multi¬ 
directional  laminates,  but  the  program  does  not  measure 
three-dimensional  effects. 

5.  In  gener.al,  the  use  of  progressive-ply-failure 

finite  element  analysis  is  a  viable  procedure  for  predicting 
the  failure  characteristics  of  composite  laminates  with 
stress  CO  ncen  t  r  .a  1 1  ons  .  It  has  been  shown  that  the 

incremental  constitutive  relations,  in  conjunction  wi ch  tne 
failure  criterion  proposed  by  Reference  CS!,  provide 
suitable  analytic  means  for  predicting  the  behavior  of 
unidirectional  and  angle-ply  laminates  under  general  stress 
states  using  only  the  lamina  property  data  obtained  under 
si.mple  load  conditions.  This  technique  has  the  potential 
for  improved  accuracy,  and  could  be  extended  to  the  study  of 
multidirectional  laminates  with  some  additional  work. 

6.  Experimentally,  the  failure  progression  due  to 
oensile  loading  for  all  laminates  was  often  aided  by  extra 
stress  concentrators  from  the  drilling  of  the  hole. 
Fabrication  damage  should  have  first  been  checked  by 
subjecting  unstressed  specimens  to  snereo  x-ray  for 
comparison  to  the  failure  progression  results.  Future 
projects  should  consider  iihis  procedure. 
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In  study  in;?  the  failure  behavior  of 


t  e  o  . 


visual  pcst-failure  techniques  are  important.  With  the 
videotape,  photojlraphs  from  the  videotape,  and  with  the 
stereo  x-rays,  failure  phenomena  were  more  closely 
correlated  than  could  have  been  obtained  with  only  the 
results  from  strain  gages. 


H!  .  Closing 

This  thesis  has  accomplished  the  ob;ective3  that  were 
listed  in  Chapter  I.  It  has  completed  a  comprehensive 
experimentation  program  that  may  help  guide  future 
experimentation  in  related  areas  such  as  high  temperature 
applications.  Also,  this  study  has  shown  that  the  computer 
program  provides  good  results  but  requires  optimization  to 
include  various  unloading  schemes. 

The  author  hopes  that  the  information  presented  i.n  thi 
thesis  may  contribute  a  small  increment  to  the  advancement 
of  materials  technology  and  the  development  of  future 
romposite  material  applications. 
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Append i x  A 

Proijress  1  on  o f  Study 

This  appendix  describes  the  logistics  effort  involved 
in  accomplishing  an  AFIT  thesis  sponsored  by  the  Flight 
Dynamics  Laboratory.  Much  prior  planning  was  required  for 
this  thesis  to  be  accomplished.  Work  began  with  the 
ordering  of  materials  from  manufacturers.  In  addition  to 
the  composite  panels  used  to  fabricate  the  specimens,  the 
supply  of  strain  gages  and  tab  material  also  required 
advanced  planning.  For  this  study,  the  time  between  placing 
the  order  and  having  specimens  ready  to  instrument  was  about 
SIX  months.  For  a  typical  AFIT  student  whose  thesis  program 
only  lasts  nine  months,  this  lead  time  can  severly  affect 
the  amount  of  work  that  can  be  reasonably  accomplished  in  an 
exper  1  ,men ta  1  thesis.  Thus,  the  following  sequence  of  events 
IS  chronicled  for  the  future  AFIT  student  who  wishes  to  do  a 
study  using  experimental  facilities  at  Wr i gh t - Pa t terson  AF3 . 

Once  the  prospectus  was  accepted,  the  material  that  was 
chosen  (Qr/PEEK)  had  to  be  manufactured  in  the  desired 
lay-ups.  In  April  1937,  there  were  two  potential  choices. 

The  Gr/PEEK  material  could  have  been  ordered  from  a  reput-able 
manufacturer,  or  it  could  have  been  manufactured  at  the 
•Structures  Test  Facility  (Bldg.  65)  of  the  Flight  Dynamics 
Laboratory.  The  press  with  which  Gr/PEEK  is  manuf acutered 
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was  not  due  to  arrive  to  Bld^.  65  until  June  1987,  and  the 
personnel  who  would  make  the  material  had  yet  to  be  trained 
Therefore,  the  choice  to  order  from  the  manufacturer  was 
made,  and  an  order  was  placed  in  May  with  a  promised  due 
date  of  late  August.  The  company  which  supplied  the 
material  was  Imperial  Chemical  Industries  (ICI) , 
headquartered  in  England.  The  material  was  ordered  from  th 
Fiberite  Corporation  in  Orange,  California,  which  is  a 
subsidiary  of  ICI.  The  order  form  was  obtained  from  Mr. 
Mick  Yardich  of  the  Aero/Astro  Department  of  AFIT;  and  it 
included  the  sizes,  quantities,  and  lay-ups  of  the  material 
desired.  It  was  submitted  to  AFIT/ENY,  which  ordered  the 
material  through  Base  Supply. 

To  determine  the  quantities  of  panels  needed,  a  test 
plan  was  written  in  which  the  number  of  specimens  to  be 
tested  was  sketched  within  the  sizes  of  available  panels. 
This  test  plan  became  the  guide  to  follow  for  all  groups 
involved  with  this  experimentation. 

This  test  plan  wa s  a c c o mp lished  with  the  help  of  Dr. 
Palazotto  of  AFIT/ENY  and  Dr.  Sandhu  of  .AFWAD/'FIBCA . 
.Addressed  to  AFWAL/FI3C,  the  test  plan  described 
experimentation  which  involved  five  different  groups 
within  FIB,  the  Structures  Division  of  the  Flight  Dynamics 
Library.  While  the  panels  were  being  manufactured,  the  tes 
plan  was  finalized,  and  each  group  involved  met  to  discuss 
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the  logistics  of  the  program.  An  outline  of  the  test  plan 
f  1 1  1  0  ws  . 

The  test  plan  was  in  letter  format,  and  the  first  sect 
listed  general  program  information  about  the  project.  The 
first  page  contained  the  following  text: 


.-ROM;  AFIT/ENY  (Lt  R.  John  Martin,  54731)  29  SE?  8 

''V3JECT:  Test  Plan  for  Master's  Degree  Thesis 

TO:  AF'AJAL/FIBC  i  Mr  .  G.R.  Holderby) 


Program  Information 


a  . 

Organization 

b  . 

Project  Number 

■O  . 

Project  Title 

d  . 

Security  C 1  ass i f i ca 1 1 on 

e  . 

Project  Engineer 

t 

Project  Advisor 

g  . 

Project  Sponsor 

h  . 

Instrumentation 

Engineer 

J  • 

Test  Engl neer 

k  . 

USAF  Fabrication 

Engl neer 

1  . 

Test  Loca t ion 

AFIT/ENY 

24010366 

"Study  of  Ply  Failure  in 
Thermoplastic  Composites 
Unci  ass i f i ed 
Lt  Martin 

Dr.  Anthony  Palazotto, 

AFIT/ENY 

Dr .  R . S .  Sandhu , 

AFWAL/FIBCA 

Mr .  J .  Pappas 

Mr.  Stalnaker/  Mr.  Bates 

Lt  Fe 1 ker 

AFWAL/Structures  Test 
Faci  1 1 ty  ,  Bldg  .  65  , 

Area  B,  WPAFB 


Next  came  a  description  of  the  test  plan,  which 
included  many  diagrams  and  figures  explaining  the  testing 
that  was  to  take  place.  The  general  test  plan  began  bv 

t 

stating  objectives.  The  remaini.ng  portion  was  devoted  to 
the  technical  description  of  the  project.  The  following 
subjects  were  addressed  by  this  technical  description: 


205 


Materia 


3  . 

4  . 

5  . 

6  . 


Stacking 
Pane  1 3 
Specimens 
Testing 
Report 


equence  of 


Specimens 


October  lOa".  t  wo  mo  n  t  h  s  after  t  h  e  v  had  been  p  r  i  m :  s  e  i  The-.' 

were  taken  to  the  Non  - Des true 1 1 ve  Zvaluatirn  Bran  on  frr 
inspection  by  C-scan.  After  it  wa  s  d  e  t  e  r  m  i  n  e  d  t  h  a  --  ear  n 
panel  contained  no  significant  flaws,  the  panels  were  taken 
to  Bldg.  65  to  the  Composites  Group  which  fabricated  the 
3  amo 1 es , 


Panels  were  cut  into  subpanels,  which  were  then  tabbed. 
Once  tabbing  was  complete,  the  subpanels  were  cut  into 
3  p  e  c 1  me  ns  according  to  the  test  plan.  Speci me  ns  were 
drilled  as  necessary  to  create  a  0.4-inch  hol.e  at  their 
centers.  The  process  of  fabricating  the  specimens  from 
the  manufacturer's  panels  took  a  total  of  a  month  and  a 
half.  As  specimens  were  completed,  they  were  taken  from  the 
Composites  Group  to  the  Instrumentation  Group  for  the 
application  of  strain  gages. 

The  Instrumentation  Group  accomplished  the  tedious  task 
of  "gaging"  each  specimen.  The  strain  gages  used  were  three 
types  of  Ml cro-Measurements  Precision  Strain  Gages: 
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1  CEA-V33- i  25'JR-350 

2,  CEA-03-062UR-350 

3.  WK-03-060WR-350 

The  CEA  gages  were  st,rain  gage  roset-tes  which  contained 
three  separate  grids  as  shown  in  Figure  25  of  Chapter  IV. 

The  WK  gages  were  stacked  rosettes,  i.e.,  the  three  grids 
were  stacked  on  one  another  to  achieve  more  pinpointed 
measurements  of  strain.  The  ‘03'  is  a  measure  of  the 
coefficient  of  thermal  expansion  of  the  strain  gage  materia* 
which  must  be  compatible  with  that  of  the  material  being 
gaged.  The  third  set  of  alpha-numeric  set  of  characters 
tells  the  size  and  type  of  gage.  For  example,  ’060'  is  one 
dimension  of  the  grids,  i.e.,  0.060  inches.  '  WR  ‘  indicates 

a  stacked  rosette.  The  ’350*  indicates  the  resistance  in 
ohms  at  room  temperature  for  which  the  gage  is  designed. 

Test  equipment  used  in  this  thesis  was  set  up  for  350  ohms. 
Any  other  resistance  in  a  gage  would  have  caused  a  major 
reconfiguration  of  the  equipment. 

Applying  strain  gages  requires  a  skill  level  attained 
only  through  training  and  experience.  The  technicians  in  the 
Instrumentation  Group  followed  a  multi-step  process  and 
demonstrated  excellent  workmanship.  Strain  gaging  the 
specimens  took  about  four  weeks  to  accomplish,  but  this  time 
period  was  partly  concurrent  with  the  testing  of  specimens. 

From  the  instrumentation  group,  the  specimens  were 
taken  downstairs  in  Bldg.  65  to  be  tested.  Total  testing 
t*me  also  took  about  four  weeks.  The  technicians  and 
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I  r.<i^r  3  ^cco  mp  lishin^  vhis  t-estiin^  we  re  e  x  t.  r  e  me  *  v  r.  e  .  r  i  - 
and  provided  en’ihusiastic  support  for  this  project. 


Data  from  the  tests  we  re  extracted  by  the  co  mp  u  t  e  r  s  .  r. 
31d?.  55.  The  engineers  <atherinig  this  data  were  also  very 
supportive  of  this  thesis.  From  this  data  came  much  of  the 
experimental  results  used  in  this  thesis. 

As  a  final  note,  the  author  feels  that  much  was  learned 
just  by  participating  in  all  the  phases  of  the 
experimentation  process.  Getting  a  hands-on  education  net 
only  helped  get  the  work  done,  but  it  allowed  for  better 
’.;nder  s  tand  1  ng  of  the  results  when  it  came  to  writing  the 
thesis  . 


Since  ihermoplaistic  romposines  have  a  semi  cry  s  na  i  ne 
matrix,  they  require  different  processing  procedures  than 
the  therraoset  matrix  composites.  The  .manufacturing  of 
-aminates  involves  the  stacking  of  layers  of  prepeg.  p^^acing 
into  a  press,  heating  to  above  the  melting  temperature  c: 
the  matrix  under  low  pressure,  applying  a  c o n s c  .  i da 1 1  :  n 
pressure,  and  cooling  rapidly  to  control  the  degree  of 
crystallinity  to  ensure  consistent  properties  [34!  . 

A  prepeg  is  a  thin  layer  (0.00525  inches-  of  the 
composite  material,  .A  laminated  panel  is  .made  up  of  several 
plies  of  prepeg  tape.  A  single  lamina  is  produced  by 
cutting  pieces  of  the  prepeg  in  the  desired  orientation  see 
F.gure  9  9  [34].),  aligning  and  butting  the  pieces  together 
and  spot  welding  each  seam  (see  Figure  100  [34 [^ 

Once  the  desired  number  of  laminae  have  been  .made,  a 


'prelam  stack'  of  eight  plies  is  prepared  bv  tacking 
together  each  ply  in  the  desired  stacking  sequence.  A 
laminate  panel  is  made  up  of  one  or  more  of  these  prel,am 
stacks  [34], 


The  prelaminated  panel  is  then  placed  into  a  mold 


cavity.  The  mold  cav 

1 1  y  13 

a  picture 

f  r  ame 

of  s 

tai  n  *  e 

steel  with  thickness 

0  .  10-0 

.15  inches 

less 

than 

t,he  :  : 

consolidated  panel. 

See  F 1 

g  ur  e  10  1, 

This 

f  r  ame 

1  s 

h  *5  d 
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S’.irr  o’^nded  bv  lavers  o:  hi^h  *  empera* ’jre  al'aminum  fo:.. 
3*ainles3  steel  ^lazinig  plates,  and  the  press  platen.  The 
glazing  plate  and  the  aluminum  foil  are  pretreated  with  a 
mo  l.d -release  agent  13  4], 

The  press  platens  are  electrically  heated  and  are 
equipped  with  temperature  control.  The  pressures  required 
for  contact  and  consolidation  are  achieved  by  computer- 
controlled  hydraulic  pistons  connected  to  the  top  press 
platen. 

To  achieve  the  proper  cooling  rate,  one  a  1 1 e r n a o  i  ■/ e  i 
to  use  a  two -stage  process  where  the  mold  stack  is  removed 
from  the  hot  press  after  consolidation  and  placed  in  a  col 
press  for  cooling.  Another  alternative  is  to  keep  the  mol 
stack  in  the  press  and  cool  the  press  platens  internally 
win.-,  an  air- water  mixture  coolant  system  with  the 
accropriate  relation  bet  we  en  air  and  wa  ter  pressure  to 
crciuce  the  cooling  rate  required. 

The  1  a.mi  na  1 1  on  /  c  ons  o  1  1  da  t  ed  process  can  be  summarized 
in  the  following  steps: 

Preheat  the  press  platens  to  the  melt  temper atur 
3  Place  the  mold  stack  between  the  platens  and 
accly  a  contact  pressure  of  T0  psi.  Heat  the  stac.x  ; 
^  m.inutes  per  each  prelam  stack  with  a  maxim.um  of  30 


3'  Apply  a  conso 1 : daP i on  pressure  or  200  psi  for 
m  p  r.  u  t  e  3  . 

o 

(4)  Cool  the  panel  at  a  minimum  of  50  F  per  minut 
achieve  the  proper  decree  of  crystallinity.  The 
cooling  pressure  should  be  300  psi. 

'51  Remove  the  mold  stack  from  the  press  after  5 
minutes . 

.5)  Carefully  remove  the  panel  from  the  mold. 
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